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Some Factors Affecting Large Transport 
Aeroplanes with Turboprop Engines 


The Thirteenth Wright Brothers Lecture* 


A. E. RUSSELLT 
The Bristol Aeroplane Company, Lld., England 


FOREWORD 


To deliver a lecture in commemoration of the Wright Brothers 

before the Institute of Aeronautical Sciences is a great honour of 
Which I am deeply aware. This honour one feels is due not 
@lely to the technical situation but to some more subtle link 
between the Institute and the Royal Aeronautical Society and 
our two countries. 
» Iam grateful to your Council for the privilege of being included 
Mtheir list of scientists and engineers from the United States and 
from England who, on this date, annually recall the first con- 
Molled and sustained flight and present some aspect of aero 
Mautics of topical interest. 


INTRODUCTION 


T {HE MAJORITY OF LECTURERS in this series have been 
scientists, responsible for the direction and control 


Mi research, who have indicated the path leading to- 
Wards improved performance of aeroplanes or power 
Plants. They have presented their work in the knowl- 
idge that they serve to educate other scientists working 
it the same field, directing bodies for the function of 
their policy and engineer in their application to 
design. 

As a experience I have 
fo report embraces a field. the 
last few years The Bristol Aeroplane Company, Ltd., 
has devoted a considerable proportion of its effort 


the 
smaller 


design engineer, 


During 


Howards the design and development of large trans- 
Port landplanes with turboprop engines. With this 
a& my subject I am aware that various papers have 
been published recently dealing with various aspects of 
the same type of aeroplane and power plant. Their 
Nope covers the field of broad conjecture, deductions 
* Presented before the Institute of Aeronautical Sciences in the 
US. Chamber of Commerce Auditorium, Washington, D.C., 
December 17, 1949. 

t Chief Designer, Aircraft Division 


along paths of established aerodynamic and thermo- 
dynamic principles, and reports on the performance of 
experimental and development aeroplanes and power 
plants. For a working basis one appears to be on fairly 
safe ground if one accepts the conclusion that, for 
transport aeroplanes, the crossover point in the choice 
between turbojet and turboprop engines is at medium 
ranges, with the turboprop giving greater economy at 
long range. 

Probably of all the world’s records the one for long 
distance bears least relationship to feasible transport 
operation. To achieve it, all normal standards have 
to be ignored; the take-off is long and hazardous, the 
aeroplane is overloaded to a point where the factors 
of safety for strength and performance standards rely 
on the cooperation of fate, and not less trust is placed 
in the survival of the motors. The sole restraint is on 
the patience in awaiting a favourable weather report. 

In civil operation one does not find the airworthiness 
authorities more sympathetic when economic difficul- 
ties arise and the pay load inevitably depends on the 
difference between the gross weight of a fully certified 
aeroplane and the sum of the unloaded weight and the 
fuel weight. 

The selection of turboprop engines offers one solution 
to the problem of long range, and the size of aeroplane 
in which they are installed must be adjusted to suit 
the load to be carried. With established technical data 
this is just a case for simple arithmetic. 

The most obvious change indicated by general investi 
gations into the effect of using turboprops is a rise in the 
economic cruising speed and altitude. Some have 
derived satisfaction from the that 
is consistent with earlier steps, for it can be shown that 
in the decade 1930-1940 the standard of performance 


of transport aeroplanes gave cruising speeds of the 


fact this increase 
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order of 200 m.p.h. at 10,000 ft.; in the decade 1940— 
1950, 300 m.p.h. at 20,000 ft.; and now in the next 
stage, 400 m.p.h. at 30,000 ft. This seems to clinch the 
issue. 

The review I wish to make is from the engineer's 
angle, for at these higher cruising speeds and altitudes, 
together with expanding requirements for range and the 
increasing size of aeroplanes, new technical problems 
have been introduced and old ones have been amplified. 
The role of the designer of transport aeroplanes is to 
find sound solutions to these problems while still satis- 
fying the needs of the operator by providing him with 
efficient load-carrying capacity and without departing 
from those standards of safety established and defined 
by past experience. 

The main section of the paper deals with certain 
technical aspects relevant to the design of large aero- 
planes. To give perspective to this part, two sub- 
sidiary sections are included. 

The first is a brief review of the characteristics of the 
turboprop engine immediately available for the type 
of aeroplane under consideration. 

The second is an analysis of the primary aerody- 
namic characteristics determining operating efficiency 
and, hence, defining the geometry of the aero- 
plane. 


(1) TurRBopRop ENGINE PERFORMANCE 


The turboprop engine possesses, in common with the 
turbojet engine, the fundamental advantages of rotating 
as opposed to reciprocating machinery—e.g., freedom 
from low-frequency vibration and noise, mechanical 
simplicity making for ease of maintenance and reliabil- 
ity, and the absence of any immediate limitation on the 
output of a single unit. The first two of these features 
have so far been amply demonstrated in service. 


In these respects the turbojet engine appears to be 
extremely attractive, offering as it does many advan- 
tages from the point of view of the air frame—e.g., 
grouped power units and lower installation drag. 
However, the saving thus attained is offset by the 
relatively poor propulsive efficiency of jet units at the 
speeds that are appropriate when long range and high 
pay load are the prerequisites. Also a handicap of the 
turbojet engine for civil aeroplanes is its comparatively 
low take-off thrust. If cruising conditions alone are 
made the criteria, a good case can be made for the 
turbojet over a considerable range of speed and distance, 
since the saving in weight is sufficient to compensate for 
the greater fuel load. However, when take-off per- 
formance is a major criterion, assuming no recourse to 
assisted take-off by special means, then a larger unit is 
required, and the weight-saving associated with jet 
propulsion is not so evident. The turboprop, on the 
other hand, can inherently provide ample take-off 
thrust, and the power-plant weight may, in fact, be of 


the same order as that for a turbojet for equal take-of 
performance. . 

The comparison between the turboprop and th: 
piston engine is more subtle, since each utilises the sam, 
final propulsive means, and the turbine, when fitty 
with a heat exchanger or using a high compressiq 
ratio, has a thermal efficiency of the same order as thy 
of the piston engine. 

The differences here mainly concern the operating 
characteristics of the two engines and are essential} 
linked with the characteristics of the aeroplane unde 
consideration. For example: 

(a) The compression ratio of the turbine engine 
not constant but increases with reduced temperature 
This leads to an improvement in efficiency with jp. 
crease in altitude. 

(b) The compression ratio and peak temperatur 
are functions of r.p.m.; this results in a deterioration 
in efficiency at part load. 

(c) The ram compression ratio, besides increasing 
the air mass flow, also enables the working fluid to e. 
pand over a higher ratio with consequent increase ¢ 
power and improved economy. 

The three effects combine to make the economic 
operation of propeller gas turbine essentially restricted 
to high altitudes, speeds, and powers. As the neces. 
sity for operating the engine at a high percentage ¢ 
maximum power increases the economical latitude and 
speed of the aircraft, these conditions can be fulfilled 
simultaneously without detriment to the aircraft 


performance. 
Factors Affecting the Performance of Turboprop Engine 


Having established the basis of the turboprop en 
gine, it remains to discuss the methods available for 
improving its overall performance as an aeroplane 
power plant. As usual in power plant design, once 
the fundamental principle has been established and 
executed efficiently further economy and performance 
can only be bought at the expense of increased com: 
plexity and weight. Thus, a balance between thes 
must be arranged to result in the minimum power 
plant weight and specific fuel consumption, with due 
consideration of other factors that may affect the 
overall aircraft performance. A fundamental re 
quirement for high efficiency is to reduce the rate o 
energy rejection in the exhaust, and it is evident that 
two alternative methods are available. Exhaust heat 
energy may be recuperated by means of a heat ex 
changer and fed back into the cycle, or the exhaust 
temperature may be reduced by increasing the ex 
pansion ratio. Another method is to raise the peak 
temperature, but this is not such a fundamental vari- 
able since it is not applicable to the ideal simple cycle. 
These devices to improve efficiency will be further 
considered, together with their effect on engine 


weight. 
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Heat Exchangers and Other Modifications to the Simple 

Cycle 

The use of a heat exchanger of moderate thermal 
ratio (about 40 per cent) raises the efficiency of the 
practical cycle to the same order as that of the piston 
engine without recourse to a compression ratio greater 
than 5:1. However, it immediately necessitates duct- 
ing the compressed air before combustion into juxta- 
position with the exhaust gases, with inevitable in- 
crease in bulk and weight of the installation and the 
creation of parasitic energy losses. A further objec- 
tion is that exchanger matrices are extremely sus- 
ceptible to high gas temperatures, which at present 
prevents the heat exchanger engine from taking ad- 
vantage of increase in the latter. Despite this, how- 
ever, there may well be a case for this type of engine in 
very long-range aircraft, especially when a satisfac- 
tory regenerative matrix has been developed capable 
of thermal ratios of up to 90 per cent. 

Intercooling and reheat are well-known devices for 
improving the efficiency of a gas-turbine cycle, but for 
aeroplane applications their complexity is likely far to 
outweigh the small decrease in fuel consumption. 


Development of the Simple Cycle 


It would seem that the immediate programme of 
development for turboprop engines should concentrate 
on improvement of component efficiencies in order to 
exploit fully the potentialities of the simple cycle when 
optimum values are prescribed for the main parameters. 
This approach will result in a light, simple power unit 
of low frontal area and attractive fuel economy. 

It is convenient at this stage to separate the thermal 
cycle from the means of propulsion, and consider the 
compressor turbine combination as a gas generator 
producing ‘“‘gas horsepower.’’ The best means of 
dividing the available power between propeller and jet 
may be examined subsequently. 

It should be noted that this procedure is satisfactory 
for turboprops, since the propulsive efficiency is nearly 
independent of the thermal efficiency; but this is not so 
for turbojet engines, and such a procedure may lead to 
erroneous conclusions. 

The specific output and thermal efficiency of the 
gas generator are functions of four main variables. 
These are compressor temperature rise (or compression 
ratio), peak temperature, compressor efficiency, and 
turbine efficiency; other parameters have fairly well- 
defined values. Fig. 1 shows the effect of the first 
two variables for representative values of component 
efficiency, where thermal efficiency is plotted against 
gas horsepower for a gas-generator unit. It is evident 
that it will pay to use as high a gas temperature as 
possible for maximum power output, and this is also 
advantageous from the point of view of efficiency, 
especially at high compression ratios. The peak tem- 
perature assumes this importance in the practical 
cycle because of the fact that the unit would require a 
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power input if no fuel were supplied and the peak 
temperature were equal to the compressor delivery 
temperature. A certain fuel input is necessary to 
overcome the sundry losses before the efficiency attains 
a positive value; this quantity becomes of less im- 
portance as the peak temperature is raised. At ex- 
tremely high values, we approach the ideal cycle condi- 
tion where the efficiency is dependent solely on com- 
pression ratio. To exploit the advantages of high 
temperature, we are dependent on advances in mate- 
rials, or, alternatively, a system of blade cooling must 
be developed. The latter expedient immediately 
introduces a compromise, since substantial losses are 
associated with the provision of a cooling flow. The 
increased power output resulting in low engine weight 
is also compromised to a certain extent, and it would 
seem that low specific weight is more likely to be 
immediately obtained by aerodynamic and structural 
refinement of the components. 

For the purpose of arriving at the best type of turbo- 
prop for the immediate future, Fig. 2 has been pre- 
pared for a fixed value of maximum temperature—1.e., 
827°C., which represents the best that can be achieved 
with current materials without recourse to turbine blade 
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cooling. A variation in compressor efficiency is in- 


corporated to allow for the fall-off likely to occur at 
high ratios because of the unfavourable effect of reheat 


in the stages and the difficulty of matching them satis- 
turbine 


The 
efficiency is retained as a variable together with com- 


factorily over the operating range. 


pressor temperature rise; the former is not neces- 
sarily dependent on expansion ratio, since the beneficial 
effect of reheat in a turbine tends to compensate for 
It will be seen that turbine 


efficiency has a marked effect both on the absolute 


the difficulty in matching. 


values of efficiency and output and on the optimum 
compressor temperature rise. 

Turning to the optimum value of the latter param 
take-off altitude 
cruise conditions may influence the value chosen. How 


eter, a compromise between and 
ever, for a turbine efficiency of 86 per cent, which is 
reasonably attainable, it is seen that the optimum 
temperature rise for efficiency is in the region of 300°C. 
(9:1 compression ratio), while the value for best power 
output per unit air flow is just over 200°C. (about 
5'/2:1 compression ratio). 

Clearly, the question of weight now arises, and the 
high compression engine will be penalised on this score 


SCIENCES FEBRUARY, 1950 
because of both the lower specific output and th 
weight of extra compressor stages and the change froy 
light alloys to steel in certain parts of the engine stry. 
ture to withstand the high temperatures. 

Before attempting to evaluate this compromise, jt js 
necessary to convert the power and efficiency figures oj 
Fig. 2 into thrust horsepower and specific fuel cop. 
sumption. The optimum division of work betwee; 
propeller and jet may be determined for any on 
operating condition, but for a specific aeroplane ap. 
plication this also is subject to a compromise with take. 
off requirements and is governed by extraneous limita. 
tions on propeller diameter, speed, etc. In this respect 
the high cruising r.p.m. of turboprops may be an 
embarrassment, but the use of an independent power 
turbine enables the problem of take-off thrust to be 
overcome by overspeeding while retaining optimum 
r.p.m. for cruise. 

Assuming the use of this device, the equivalent 
powers and fuel consumption have been estimated for 
the cruise condition of Fig. 2 and a turbine efficiency 
of S86 per cent. Together with a curve of estimated 
weight against compressor temperature rise, this en- 
ables Fig. 3 to be drawn. 

This is a plot of engine plus fuel weight versus com- 
pression ratio for a range of durations, all at 30,000-ft 


390-m.p.h. cruise conditions. It can be seen that the 
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LARGE TRANSPORTS 
rinimum engine plus fuel weight is achieved at a com- 
pression ratio of about 8:1 for durations above 10 hours. 
Up to durations of 5 hours, the importance of engine 
weight predominates, and compression ratio is not 
critical; low compression ratio engines are then opti- 
mum from the point of view of simplicity and initial 


cost. 


The foregoing is still not the complete story on the 
selection of the turbine-engine parameters. The varia- 
tion of power and consumption with altitude and speed 
is aflected by the choice of compression ratio and peak 
temperature. For instance, a high-compression engine 
is more susceptible to ambient temperature than one 
of low-compression ratio, so the high-compression 
engine will show more improvement with altitude and 
more under tropical conditions. High 


temperatures, on the other hand, reduce the effect of 


power loss 


ambient temperature slightly. 


There is as yet not sufficient evidence on the sub- 
ject of scale effect to decide on tke optimum size of 
unit. Fundamentally, the and length of 
the blade-carrying elements will vary as the square 


diameter 


root of the power for a given design, so the specific 
weight would be expected to increase with size. How- 
ever, manufacturing limitations tend to make small 
engines heavier, and over the present range of propeller 
turbines the specific weight of similar designs appears 
to be roughly constant. It does seem probable, how- 
ever, that component efficiency may be expected to 
increase slightly with engine size. In certain cases it is 
beneficial to couple two small engines together 
for the Brabazon Mk. 2 the coupled layout is more 
compact and more suitable for buried installation than 


In this particular case also, the 


e.g., 


one large engine. 
choice of engine compression ratio is clearly near the 
thermodynamic optimum, since the range required 
lays all the emphasis on specific fuel consumption. 


(II) AEREROPLANE PERFORMANCE 


When reciprocating engines are employed and a 
normal relation is maintained between maximum con 
tinuous boost and maximum economic boost, the power 
available at low altitude (5,000 ft.) at the former 
condition is severely limited; this arises from the fact 
that it is desirable to arrange that the full throttle 
height at cruising r.p.m. is at least as great as the de 
sired cruising altitude, usually in the region of 20,000 
ft. 

The take-off power is also limited by the basic char 
acteristics of this type of engine, though to a lesser ex 
tent than the maximum continuous power. 

These conditions will often result in take-off per- 
and airworthiness performance standards 
having an influence on the choice of wind loading and 


formance 


aspect ratio; the necessary compromise resulting from 
consideration of those items of the performance may 


WITH T 


URBOPROP ENGINES 71 
have an adverse effect on the pay load/range perform- 
ance. 

For aeroplanes with reciprocating engines, the pay 
load/range characteristics improve with increasing 
wing loading, and such increases have been made, 
through increase of weight, during the development of 
all successful aeroplanes, since means have been found 
to circumvent critical conditions. 

A different situation obtains when turboprop engines 
are employed; this arises from the fact that the power 
output of a turbine engine increases continuously as the 
altitude is reduced from cruising altitude to sea level, 
and it is unlikely that aeroplane performance at the 
lower levels will be important to the extent of deter- 
mining gross weight or wing design. For turboprop 
engined aeroplanes, attention will therefore be directed 
towards attaining the best possible efficiency under 
cruising conditions. 

When turboprops are used, high cruising altitudes are 
essential; there is no other way to obtain an acceptable 
pay load range performance; any other advantages 
that may accrue are incidental, and any penalties in- 
volved are small in relation to the gain in economy and 
must be accepted. Also, it is essential to operate the 
engine at the maximum permitted continuous r.p.m. 
and power in order to achieve the minimum specific 
consumption, and it follows that at the lower alti- 
tudes, where the minimum power required is less than 
the power available, the equilibrium forward speed is 
far above the economic speed. The only way to re- 
duce the power available so that an economic condi- 
tion is obtained is to operate at the highest altitude 
consistent with the maintenance of sufficient reserve 
climb performance. This altitude increases as weight 
is reduced by the consumption of fuel and is rather 
nearer to the ceiling than to the altitude for maximum 
speed. 

Thus it is necessary to cruise at a speed closer to the 
most economic speed than when reciprocating engines 
used. It follows immediately that the induced 
drag—and, hence, the wing span 
more dominant role when turboprop engines are used. 


are 
plays an essentially 


Cruise Procedure 

When the cruising r.p.m. of the engines are main- 
tained constant, the performance of an aeroplane at dif- 
ferent weights and over a wide range of altitudes may 
be presented in the form shown in Fig. 4, the merit of 
which is that the effect of alternative flight plans on 
speed, specific range, cruising altitude, and reserve 
club performance is immediately apparent. This form 
of presentation was devised by the B.O.A.C. and was 
first given by Dykes in his paper at the Second I.A.5.- 
R.Ae.S. Conference.° 

Table 1 compares the result of following various 
flight plans, since the weight of the aeroplane falls be- 
cause of the consumption of fuel. The initial weight 
is not important in this comparison and has been taken 
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TABLE | 

Flight Mean Cruising Still-Air Range, Fuel for Difference Maximum 
Plan Speed, M.p.h. Miles 5,000 Miles from Plan (3) Altitude, Ft 
(1) 358 4,900 0. 386W +0.048W 33,000 
2) 319 5,900 0.328W —0.010W 39,000 
(3) 342 5,700 0.338W qe 39,000 
(4) 348 5,100 0.372W +0.034W 30,000 


as 300,000 Ibs. The still-air range (column 3) assumes 
that the fuel available amounts to 37.5 per cent of the 
The flight plans examined are: (1) 
(2) speed for 


initial weight. 
maximum cruising speed throughout; 
maximum range throughout (nearly zero potential 
rate of climb); (3) altitude such that the potential 
rate of climb is 200 ft. per min.; (4) as in (3), except 
that the maximum altitude is limited to 30,000 ft. 

[An additional plan would be to fly at constant E.A.S. 
throughout a flight; depending on the speed chosen, 
the result lies somewhere between Plans (2) and (3).] 

It is seen that with one exception the differences be- 
tween the mean speeds are relatively small; on the 
other hand, the available range with a given quantity 
of fuel or the weight of fuel required to fly a given dis- 
tance varies in a most significant way; for, if the pay 
load of a fully equipped civil aeroplane be taken as 
0.10W for Plan (3), then for Plan (1) it would be only 
one-half that amount. 

The reduction in pay load or range due to adopting 
Plan (3), instead of Plan (2) (which gives the greatest 
range), is seen to be small, whilst giving a useful in- 
crease in speed. Plan (2) would probably not be ac- 
ceptable because of the small potential rate-of-climb 
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available, whereas a rate-of-climb of 200 ft. per min, 
available if speed is reduced, should be sufficient ty 
maintain a consistent flight path. The conclusion js 
that the best plan is No. 3, flight at constant potential 
rate of climb. It may be noted that this corresponds 
also to almost constant lift coefficient. 

If towards the end of a flight it is necessary to re. 
strict the altitude to, say, 38,000 ft., for reasons of 
limitation of cabin differential or due to the perform. 
ance of the air-conditioning system, the penalty js 
small. 

Plan (3) will be used for subsequent investigation oj 
the effects of other variables on long-range operations 
including the effect of engine failure. 


Initial Weight 


Fig. + can also be used to show the effect on range 
of variation of the initial weight for the particular 
aeroplane to which it relates. Starting at a series oj 
different weights, the fuel consumed in flying different 
distances can be determined, A plot of the results 
would show that, for a wide variation of the initial 
weight, covering all likely values, the fuel consumed for 
a given distance is a constant percentage of the initial 
weight. 

This result is also indicated by the Breguet range 
formula, which has as its basis flight at constant life 
coefficient. 

It is clear from Fig. 4 that the mean cruising speed is 
sensibly independent of initial weight. 

Thus, for a fixed geometry of aeroplane, engine 
power, and set range, the potential pay load is approxt- 
mately a constant percentage of the gross weight, such 
tendencies as there are being favourable with an in- 
crease of gross weight in spite of the lower operating 
altitudes. 


Effect of Temperature 


If an aeroplane is required for world-wide service, 
it is necessary to take account of the effect of tempera- 
ture on both the power output of the engines and on 
the power requirements of the aeroplane. 

The influence of ambient temperature on the per- 
formance of a typical aeroplane with all engines operat- 
ing is shown in Fig. 5. The three full lines show specific 
range in ICAN Standard conditions and in conditions 
approximating to temperate maximum conditions 
(15° and 26°C., respectively, higher than ICAN Stand- 
ard). The dashed lines refer to cruising altitudes and 
cruising speeds are marked on the full lines. 
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LARGE TRANSPORTS 

In this example the cruising altitude comes down by 
about 3,000 ft., and the range and speed by approxi- 
mately 5 per cent, for each 10°C. increase above ICAN 
Standard. For temperatures lower than ICAN Stand- 
ard, the performance is, of course, improved to a corre- 
sponding extent. 

Thus it appears that the ambient conditions in flight 
have a greater influence on the performance of aero- 
planes with high-compression turboprop engines than 
has been the case with reciprocating-engined aero- 
planes. 

A partial solution to the problem this presents is to 
permit an increase in the cruising r.p.m. of the turbines 
under high ambient temperature conditions; the engine 
manufacturer and those responsible for certifying ap- 
proval of his product are then set a difficult task in 
formulating a representative type test. The stress 
conditions in a turbine engine are complex and depend 
on blade temperature, gas loading, and r.p.m. 

At the other end of the scale, it may be necessary to 
restrict the r.p.m. in substandard temperatures. 

Investigations show that the relation between poten- 
tial rate-of-climb and cruising lift coefficient is little 
affected by change of atmosphere, and, as constant 
potential rate-of-climb implies almost constant lift 
coefficient, the selected flight plan can be followed by 
the crew in any atmosphere from the indications of an 
incidence meter; too large an incidence would imply 
too great an altitude and vice versa. 
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Engine Failure 


Certain performance requirements in a high am- 
bient temperature after the failure of an engine may 
be an important design criterion. This can arise if the 
air-line operator demands that a certain minimum alti 
tude must be maintainable, a necessary precaution 
during certain climatic periods or conditions, such as the 
monsoon, when the phrase “flying above the weather”’ 
has more than usual popular appeal. Such require- 
ments are of some importance in the initial selection of 
the engine characteristics, for a turboprop engine de- 
signed as an optimum for ICAN Standard conditions 
will not necessarily be the best design if operation 
through the tropics is involved. 

Fig. 6 shows the performance of the aeroplane to 
which Fig. 5 relates after the complete failure of a 
power unit, and Table 2 brings together typical values 
from both figures to show the effect of power-unit 


failure; the effect of ambient temperature is also 


TABLE 2 
ICAN 
Stand ICAN ICAN 
Atmosphere ard +15°C. +26°C 
Mean speed (m.p.h.) 350 332 313 
All Engines Mean altitude (ft.) 36,000 32,000 28,000 
Operating | Mean specific range 
(relative) 100% 95.5% 90.5% 
Mean speed (m.p.h.) 307 287 27 
One - quarter | Mean altitude (ft.) 30,000 25,000 20,000 
Engines In- < nn eniiiice 
atten ) uviea specinec range 
operative 4 ~~ . - f 
I (relative) 95.5% 90.5% 85.5% 
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brought out clearly. The values given relate to the able effect, for a reduction in induced drag leads to bot} 


- : , ; : Airw¢ 
performance during the period of consumption of the — higher altitudes and higher cruising speeds. Table 4, j olve 
e e ° . ° ° e . ° s 7 Vv 
last half of the fuel carried. It will be seen that the which the aspect ratio is varied at constant wing are It | 
effect of power unit failure on range is the same as the | summarises the results. os 
. ° saeert 3 . {rom 
effect of a rise of 15°C. in ambient temperature but that adit 
Ee : , Addi 
the power unit failure results in a much greater loss of TABLE 4 mage 
y 
speed. Mean Mean there: 
Aspect Speed, Miles Fuel for Range of Cruising a 
Wing Design Ratio M.p.h per Lb. 5,000 Miles Altitude, Ft for 1 
8 340 0.048 0.38817 23,000-37,000 possi 
Having found that the turboprop aeroplane must be 10 345 0.048 0.348W 26,000-39,000 gine | 
mae 12 350 0.054 0.3091 30,000-41.090 
aT > . solyy ¢ ce > . ~ et > , , 
operated closely to a prescribed cruising procedure, chutt 
the effect on operating efficiency of varying the wing , 
: ; ae Power Loading t 
area and the aspect ratio can be examined. The actual _ ,; 
. . : ° sc ° ‘ : yay 
selection of these values will, of course, depend on a While turboprop engines are in their present stage of “ 
t i ae ne § 
balance between aerodynamic gains and structural development, the power produced and the specifi . . 
; . ; : ; ee impo 
penalties. fuel consumption achieved at the time of initial ser 
oe ee ee nae oo names weig] 
Thus a limit is imposed on span by the amount of ice may differ from the values used in the aeroplan 
structural stiffness than can profitably be provided to design calculations. II 
inhibit catastrophic aeroelastic effects; this is a func A change in turboprop engine power may be con 
tion of the proportions of the wing box, which are in sidered in two stages, the effect at a constant gros Dyna 
turn determined by the aspect ratio. At the relevant weight and the effect of the gross fixed weight can b 0 
speeds there are major technical obstacles to designing creased. 
; aaa ; pr eee . a gross 
a wing of very high aspect ratio without incurring an At a fixed initial weight, the primary effect on rang ; 
‘ ae ‘ . : : . ou 
uneconomic weight penalty. and pay load of an increase or decrease in engine power “— ; 
; . ene a ae ae boas. See fund: 
As primary technical considerations do not impose a "> P rhaps surprisingly, nil. Phe truth of this is best sider 
limit to wing loading, the effect of this variable on demonstrated by the following simple analysis: iti 
pay load range efliciency can be examined at constant —__ At constant Cz, Cp, is constant and o1” is propor aye 
aspect ratio without invalidating the conclusions. tional to IV, or satie 
The effect of a change of wing area on a given aero Va (W/o) stres 
plane is to displace the selected cruise procedure line stati 


: int Hence D,IV and d.hp. (or e.b.hp. if propeller efficiency is 
both horizontally and vertically, the cruising speed “a of tl 
rae Pad Sart. ‘ : bia constan 

rising and the specific range falling with increase of wing 


loading. Table 3 gives a summary of the results of c(W*/e Norn 
varying the wing area at constant aspect ratio. ; ores ‘, 
; Hence, assuming that the specific fuel consumption 1s Ci 
constant, lbs. perhour. «(IW ¢°). plific 
TABLE 3 Then specific range (= J Ibs. per hour) « (Il ¢ since 

Initial Mean Mean (0° WW’) or « IV’~', so that, for constant specific fuel lor 
Loading Moh Rene a a consumption and the selected flight plan, the specific usa 
15 298 0.0500 0.3331 29,000-42,000 range is a function of weight only. In fact the specifi — 
55 341 0.0480 0.34811 26,000-41,000 fuel consumption falls slightly with increasing alti arbi 

65 349 0.0463 0.3621 24,000-40,000 me : , a ‘ 2 

os ans 00445 0 37811 29’ 0100-39009 tude, and, if an increase in engine power at a given alti at 
tude can be effected at constant specific consumption, As 1 

the provision of greater power will result in a small bec 
It is seen that a reduction of wing loading raises the improvement to the specific range. If, on the other aero 
cruising altitude and, hence, the economy, but the hand, the provision of greater power is accompanied by reac 
amount of the improvement is comparatively dis- 4 small general reduction of specific consumption, th dyn 
appointing. The reason is that the increase in drag yalye of the extra power .in respect of range will be lor 
coefficient has a large detrimental effect. As will be yyore marked. Where long range is of paramount app! 
shown later when allowance is made for the different —jnterest, a moderate reduction in the specific consump T 
wing structure weights at the several loadings, it is tion is of more value than a large increase in power. that 
apparent that there is no optimum value of the initial It is not so easy to generalise on the value of an asst 
wing loading in relation to pay load range performance. — jnerease in power in improving range or pay load via an iner 
Over the range of loadings tabulated, the sum of the  jyerease in all up weight, since the magnitude of the obt: 
weights of the fuel required and of the wing structure is — increase, the characteristics of the particular aeroplan¢ of a 
substantially constant. concerned, the question of the provision of sufficient nat 
Purely from performance considerations, an increase structural strength at the increased weight and with reir 


of aspect ratio can be expected to have a more favour- the new distribution of weight, and a whole series 0! The 
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LARGE TRANSPORTS 


Airworthiness and Operational requirements are in- 
volved. 

It is clear that the circumstances are quite different 
from those obtaining when piston engines are employed. 
Additional power to permit an increase of gross weight 
by overcoming some critical operating condition and 
thereafter throttling back to a more favourable power 
for fuel consumption in cruising flight is no longer 
possible. Also, in the example quoted where one en- 
gine was inoperative, fuel is not saved by deliberately 
shutting down an engine or-group of engines. 

It appears that for any long-range aeroplane the 
pay load will be substantially a constant percentage of 
the gross weight and that the only features of prime 
importance in a gas-turbine engine are the installed 
weight and the specific fuel consumption. 


(III) ANALYSIS OF STRUCTURAL CHARACTERISTICS 


Dynamic Effects 

On being faced with the design of an aeroplane of 
gross weight far beyond our previous experience, some 
thought was given to the need for including certain 
fundamental effects that had hitherto not been con 
sidered in normal methods of design. With greater 
size, greater masses, and greater structural deflections, 
it was apparent that the dynamic response under the 
action of rapidly applied external forces could produce 
stresses different from those estimated by the usual 
static analysis. For this investigation a knowledge 
of the behaviour of structure in vibration was required. 


Normal Modes of Vibration 

Calculations involving dynamic response are sim 
plified by the use of the normal modes of vibration, 
since elastic and inertia couplings are absent and, 
for the purpose required, aerodynamic coupling is 
usually so small that it can be ignored. It should be 
remembered however that the use of these modes is 
arbitrary in the sense that the elastic and inertia terms 
are combined irrespective of the aerodynamic terms. 
As the speed is increased the inertia-elastic coupling 
becomes progressively of less importance than the 
aerodynainic-elastic coupling, so that a point may be 
reached at which it is preferable to use a series of aero- 
dynamic-elastic (i.e., divergence) modes. However, 
for conventional subsonic aircraft the normal mode 
approach is likely to be the more accurate. 

The most important property of a normal mode is 
that, so long as structural damping is absent (a valid 
assumption for aircraft), its shape is fixed for any given 
inertia and elastic system. So-called resonance modes, 
obtained by applying an oscillating force, are composed 
of a number of normal modes, one of which predomi- 
nates. The contributions of the lesser modes must be 
removed before a true normal mode can be determined. 
The three methods of determining the normal modes 


With TF 


ENGINES 


URBOPROP 


of vibration of an aircraft structure—namely, by cal- 


culation, model tests, and full-scale tests—will be 
disctissed. 
Calculation of Normal Modes 

The inertia and elastic forces that determine the 


shape and frequency of a normal mode can be related 
by a simplified form of Lagrange’s equation® for each 
coordinate. 

Since an aircraft structure is continuous, the number 
of generalised coordinates and, hence, the number of 
equations are infinite. An approximation to the elastic 
and inertia properties of a structure must be made by 
considering a limited number of discrete masses con- 
nected by flexible but massless rods. The number of de- 
grees of freedom (or the number of masses) chosen 
depends on the complexity of the structure and the 
number of modes to be determined. If only the first 
or fundamental mode of a wing is required, it is normally 
sufficient to consider six degrees of freedom (i.e., three 
masses per wing), but if four to six normal modes are 
needed, about 20 degrees of freedom must be consid- 
ered. Since the work involved in the calculation is 
roughly proportional to the square of the number of 
degrees of freedom, the computing time for six modes 
is about ten times that required to determine the 
fundamental mode only. 

Having idealised the aircraft structure into the re- 
quired number of masses and stiffnesses, the equations 
of motion are solved by an iterative process. A mode 
of vibration is assumed, and the inertia forces corre- 
sponding to this mode of vibration are determined. 
These inertia forces when applied to the elastic system 
will give a deflected mode in general different from that 
The process is repeated until the 
At each stage 


assumed originally. 
assumed and derived modes coincide. 
the modes are made ‘‘free-free’’—1.e., the centre of grav- 
ity of the aeroplane is given a displacement such that 
the total energy of the system is zero. For modes 
above the first, the contributions of all modes lower 
than the one considered are removed during iteration 
by applying the orthogonality condition. 

This method has been in use for some years at The 
Bristol Aeroplane Company, Ltd. The computing 
time has been shortened by the use of influence coef- 
ficients for the displacements due to inertia forces; 
more recently, it has been further shortened by the 
use of matrices and, although the matrices involved 
may be of the 20th order, the method leads to quicker 
and more accurate computing. The development of 
electronic calculating machines may eventually reduce 
this computing time even further and enable more 
detailed calculations to be undertaken. 

Whereas in calculating the first few modes of vibra- 
tion it is sufficient to consider the flexural and torsional 
stiffnesses of the wing and the flexural stiffness of the 
body in connection with the masses, if modes above the 
third are required, the effects of shear deflection and 
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shear lag must be taken into account. Ina particular quent tendency towards overriding aeroelastic criteri, 


ure il 
application, the inclusion of shear deflection in the it may well be necessary to devote this time to ob. ere fe 
calculation reduced the frequency of the second un- taining reliable values of the elastic properties of thy Som 
coupled flexural mode by about 5 per cent, the third by _ structure at an early stage. sure a 
13 per cent, the fourth by 21 per cent, and the fifth by A vibration model must be supported in such a way § esult: 
30 per cent. The corresponding reduction in the fre- that the “free-free’’ modes are obtained, and this js | with |: 


quency of the coupled modes was about 15 per cent usually done by suspending the model on long rubbe J ynreli: 
on the fifth mode. The shape of the modes was also cords such that the rigid body frequency of the mode § jn tha 


affected but not to such an extent as the frequencies. jg small compared with its lowest fundamental fre. J) ness 0 

In the case of fuselages where shear stiffness is a quency. A ratio of less than one-third is necessary J methc 
greater proportion of the total stiffness, the reduction if there is to be no interference of the modes regardles J spring 
of frequency is even greater. Again, in a particular of the position of the points of support. giving 
application’ the frequency of the first mode was re- As with the calculations so with the model, the ac. § port t 
duced by 8 per cent, the second uncoupled flexural curacy of the results clearly depends on reproducing The 
mode by 34 per cent, the third by 49 per cent, and the _ the relevant stiffness characteristics and, again, if only an in 


fourth by 60 per cent. Fortunately, the effect of shear the first three modes are required, the model can be J thoug 
deflection can be included in the calculations without relatively crude. If accurate knowledge of higher |) stant, 


much difficulty by modifying the matrix of influence modes is required, there seems to be no alternative to J quenc 
coefficients. building a fairly complete scale model of shell form } to th 
The result of including the effect of shear lag across such that shear deflections and shear lags are repro. alterr 
the wing box’ and around the body periphery is again a duced. When a high degree of geometric simularity is J field, 
reduction of frequency, the reduction being of the maintained, there is, of course, no need to indulge in J a nun 
same order as that due to shear deflection. However, complicated estimates of structural stiffness. ing, | 
shear lag seems to have little effect on the shape of the This type of model would be more difficult and Excit 
modes, and it is probably sufficient to allow for its costly to make owing to the limitations of the minimum } the ¢ 
effect by a final adjustment to the frequencies. possible thicknesses of important members. Increasing | that : 
The introduction of shear deflection and shear lag, the length, mass, or frequency ratios or reducing the J» &! 
together with a greater number of degrees of freedom, modulus of model to full scale does not change the speed 
including perhaps engine nacelle and tail-plane flexibili- order of the overall problem, which becomes more dii- there 
ties, makes the calculation of the higher modes of vibra- ficult the larger the aircraft. at th 
tion a lengthy process. The sixth mode is probably Thus, while such a detailed model may be justified An 
the upper limit to which calculation can normally be for certain research investigations, it is not suited to | lect 
applied. For the higher modes there is a distinct ad- routine design analysis. For this purpose it is neces- veloc 
vantage in the use of a vibration model. sary to rely on experience to modify the results ob- since 
tained from a simple model or to carry out calcula- tiona 
The Vibration Model tions that include the factors appropriate to the modes that 
By building a model that reproduces the elastic and required. fund: 
inertia properties of the aircraft to some suitable scale the § 
and by applying an exciting force to this model, it is Resonance Tests on Complete Aircraft that 
possible to measure the frequencies and the shapes of It has been the practice for some time to carry out vm 
the normal modes in much the same way as the modes ground resonance tests on prototype aircraft prior to “a . 
are derived from resonance tests on complete aircraft. their first flight.? This serves as an experimental check vo 
Most vibration models used so far have been of qa on their flutter and vibration characteristics. This is - 
simple design. The flexural and torsional stiffnesses 4 late stage to discover, for the first time, that some onl 
have been simulated by suitably shaped bars to which intrusive mode asserts itself in embarrassing company. fine 
concentrated masses are attached at a few stations. Even the results of ground tests cannot be accepted com 
The number of discrete masses can be somewhat higher Without reserve, for the stresses induced by the vibra- dies 
than the number used in calculations, with consequent tions are often insufficient to promote the buckled lead 
improvement in accuracy. However, the reliability of State of the skin that may be appropriate to flight in 
the model results do depend on the assumption that the conditions. However, I am not aware of an instance i 
elastic properties of the aeroplane can be estimated ac- Where this discrepancy has been a just reason for optt stro! 
curately. This is a matter of some difficulty, especially | ™ism in critical cases. whe: 
if the structure includes cutouts or other discontinuities. As with the model, the aircraft must be supported tile 
Also, in the case of sweptback wings where the old in such a way that there is no interference with the cont 
methods may no longer suffice, the time required to ap- lowest frequency modes.° This condition becomes pha 
ply the more comprehensive theories is considerable. increasingly difficult to meet as the size of the aircraft ae 


However, with increasing performance and the conse- increases. The simple expedient of reducing the pres- ireq 
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sure in the tyres of the landing wheels is satisfactory 
only for small aeroplanes. 

Some aircraft have been supported on very low-pres- 
sure air bags with varying degrees of success. Good 
results have been obtained in the higher modes, but 
with large aeroplanes the lower modes have been quite 
unreliable. The bags so far developed have a defect 
in that they have either inadequate lateral shear stiff- 
ness or too great a vertical stiffness. A more promising 
method is to use high-pressure low-friction liquid 
springs with a large reservoir. These are capable of 
giving very low frequencies and can be used to sup- 
port the aircraft at the jacking points. 

The most simpie way of providing excitation is by 
an inertia vibrator producing sinusoidal forces, al- 
though, unless the speed of rotation is maintained con- 
stant, harmonies may be present at very low fre- 
quencies. Moving coil exciters, in which a coil attached 
to the aeroplane is supplied with variable frequency 
alternating current and vibrates in a strong magnetic 
field, have the important advantages that by employing 
anumber of exciters, coupled electrically to give phas- 
ing, the separation of particular modes is facilitated. 
Excitation at low frequencies is still a problem with 
the electric exciter, which has the further disadvantage 
that an external rigid mounting is necessary. Lateral 
excitation is thus more difficult. An accurate steady 
speed control is essential for either type of exciter, since 
there is a tendency for the speed to seek and remain 
at the frequency of a powerful mode. 

Amplitudes of vibration are usually measured by 
electrically integrating velocity or acceleration. The 
velocities pickup is more accurate at low frequencies, 
since the acceleration pickup response, being propor- 
tional to frequency, tends to exaggerate any harmonics 
that may be presented in the motion and suppress the 
fundamental. The acceleration pickup, however, has 
the same advantage as the inertia type exciter—.e., 
that it requires no external rigid support and it may be 
easily fixed to the aircraft in any position. Coupled 
with an efficient band-pass filter to remove harmonics, 
the acceleration pickup shows distinct advantages over 
other types. The accurate measurement of phases is 
the most necessary (and the most difficult) part, since 
the motion will be a combination of several modes.°® 
One method is to use a double-beam oscillograph for 
comparing the wave traces from the pickups and ex- 
citer. Another is to use a single-beam oscillograph 
feeding in the two-wave traces in perpendicular direc- 
tions to produce a Lissajous figure from which the phase 
can be deduced. Both methods become inaccurate if 
strong harmonics are present in the wave traces, as 
when there is motion in a plane perpendicular to the 
axis of the pickup. The most accurate is to record 
continuously the vibration on a film and measure the 
phases. The use of a multichannel amplifier and re- 
corder enables an accurate recording of the amplitude, 
frequency, and phase of the vibration of a dozen or 


more pickups to be made simultaneously. Analysis of 
this film on the drawing board is less liable to errors 
than direct measurements off a flickering oscillograph 
display. As so frequently happens, the best method 
takes the longest time. 


Flight Resonance Tests 


The fundamental symmetric and antisymmetric nor- 
mal modes can be found by disturbing the aircraft in 
flight (e.g., by applying sudden elevator or aileron move- 
ments) and recording the subsequent oscillation. The 
presence of aerodynamic stiffness, damping, and iner- 
tia will modify the normal mode, and a correction must 
be applied to the results to produce the true ‘‘vacuum”’ 
modes. Higher modes require an inertia exciter and, 
being more heavily damped, may be difficult to deter- 
mine. In any case, the use of an exciter is dangerous 
if there is any possibility of flutter at the exciting fre- 
quency. The method will thus usually be restricted 
to the first two modes where it has distinct advantages 
(in particular the absence of support difficulties and the 
more normal loading of the structure) over the ground 


resonance test. 


Dynamic Response 

In a structure, if the time of application or variation 
of the external loads is of the same order as the natural 
periods of vibration, the effects of flexibility can be 
important. In the case of a simple mass system an 
appropriate rate of loading will produce a maximum 
spring force twice that which would occur if the load 
were applied slowly. Also, for a harmonic force ap- 
plied to the same system, the general solution to the 
equation of motion provides a condition by which, 
where the force equals a natural frequency of the 
system, the forces may be limited only by damping in 
the system and by which, where the time of applica- 
tion of a transient impressed force is of the same order 
as a natural period of vibration of the system, the forces 
may be increased by up to 100 per cent. Thus, for 
calculating the loads in a flexible aeroplane, knowledge 
is required of the natural frequencies and modes (or the 
corresponding deflected shapes) in free vibrations. The 
description of each mode will depend on the relative 
importance of flexure and twist in its component parts. 

If wing and body flexure and wing torsion only are to 
be considered, the displacement of any point on the 
aircraft can be expressed in terms of two quantities: 
(1) the vertical displacements, along a fore-and-aft 
centre line or along a lateral axis, corresponding to unit 
displacement in that mode and (2) the rotations of 
wing sections about the lateral axis for the same unit 
displacement. The shear, bending moment, and torque 
in the wings for unit amplitude of vibration of the 
aircraft in any mode can then be calculated. 

In a free vibration, the “inertia forces’ at any time 
are equal and opposite to the ‘‘elastic forces.’’ There- 
fore, given the displacement of the reference section in 
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any mode, the elastic effects may be found as those due 
to the inertia forces arising in an equivalent vibration 
of amplitude equal to the displacement. These inertia 
forces can be used in conjunction with displacement 
at the reference section to give the equivalent loading 
on the wing or body at any time during or after a given 
set of forces. 

By considering the displacements defined by the 
mode considered, the inertia forces can be integrated 
inwards from the wing tip to produce the inertia shear, 
bending moment, and torque. These unit values may 
be tabulated or presented graphically. 

The response of any given mode (i.e., one degree of 
freedom) to an applied external force is found in a 
similar way to the response of a single mass-spring 
system. The basic relations are identical, since the 
displacements at all points in the structure in a given 
mode are fixed by the deflection at the reference section. 
The strain energy is equal to the work that would be 
done by slowly applied forces that produced the actual 
deflection at the reference section; it is also equal to the 
maximum kinetic energy of free vibration (of the same 
amplitude) in that mode. 

Thus, to find the ‘static response” of the mth mode 
to slow application of a force P, we write 


n= P C,/E, 


where 6n is the amplitude of the mode at the refer- 
ence section, (PC,) is the generalised value of the 
force, and F, is the generalised elastic stiffness of the 
mode. 

If the force P is applied at some place and the struc- 
ture is constrained to deflect with unit amplitude at the 
reference section in the mth mode, then the displace- 
ment at the point of application and in the direction 
of the force is defined as the “influence coefficient”’ 
C, of that force in the nth mode. 

If the forces are applied slowly, the deflection in each 
mode gradually builds up to a constant value. When 
these values are reached, all parts have equal accelera- 
tions, and the deflected shape is given by the sum of the 
deflections in all possible modes. But for normal 
structures where the masses and stiffnesses are 
easily separated, there is an infinite number of modes. 
While in such a case it is not possible to calculate the 
deflected shape by using these modes, the effect of 
rapidly applied or varying forces on the deflected 
shape and on the internal loading can be found by 
considering different modes separately. The ‘‘static’’ 
deflected shape of a structure under uniform accelera- 
tion may be simply found as the elastic deflection under 


not 


a balancing system of applied and inertia forces. 

With a spring mass system, if the variation of force 
with time is known, the momentum 
and the velocity changes to the mass can be integrated 
so that displacements due to successive impulses can be 
found. The free vibrations set up produce displace- 
ments in excess of the static deflection and the amount 


increments of 
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of this excess is called the overswing. The static for 
that would be required to produce this deflection ; 
salled the ‘equivalent static force,’’ and the differeng, 
between this and the impressed force at any time 
called the ‘‘true dynamic loading.” 


The loading applied to the structure (i.e., inerti; 
shear, bending moment, and torque) can be found }y 
multiplying the equivalent values as calculated for uni 
amplitude of vibration by the total deflection referre 
to the datum. Separating this deflection into its ty 
parts, the static deflection for all modes of vibration js 
identified with the inertia loading due to the acceler 
tion of the centre of gravity produced by the external 
forces at the time. The overswing for all modes give 
the total dynamic loading. 
effect becomes small for modes of high frequency. 


In practice, the dynami 


It may be found that the maximum nett values oj 
shear and pitching moment on the aircraft are not zero, 
in which case there are errors in the position of the 
reference datum. These nett values may be used ty 
produce linear and angular accelerations such that the 
final resultant inertia forces and moments are made 
zero and the “‘free-free’’ condition is satisfied. 


Applied Forces 


‘Applied unsteady forces to an aeroplane in flight 
may arise from either an intentional change in the 
Taking 
the case of an aeroplane flying into an uprising current 


flight path or from atmospheric disturbances. 


of air, the buildup of aerodynamic loads is complex 
During the time the aircraft is accelerating vertically, 
the vertical velocity thus acquired reduces the é 
fective gust velocity. A certain amount of pitching 
may also occur, depending on the stability characteris- 
tics of the aeroplane, though this is not a major effect 
perhaps as regards the tail plane. 

The Wagner effect is'' also important. A_ sudden 
change of incidence does not immediately produce the 
lift appropriate to steady-flow conditions; there is a 
time lag during which the system of vortices is built up. 
The buildup is asymptotic, and a measure of the rate 
is such that the wing must traverse a distance of about 
three wing chords to build up SO per cent of the full 
lift increment. By slowing down the application oi 
the gust load, this lag acts as an increase of gust gradi- 
ent distance by allowing more time for vertical velocity 
to accrue. 

A theoretical solution giving the lag for a sudden 
change of incidence in two-dimensional flow has been 
worked out by von Karman and Sears'* and for the 
lift in a sharp-edged gust by Sears and Sparks.' 
Further work by Jones''! deals with wings of finite 


aspect ratio. 


The Gust Case 


In a gust the three main factors 
response, vertical motion, and Wagner effect 
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complicating the problem on their own, suffer as well 
irom considerable coupling. A short list of these 
couplings is as follows: 

(1) The change of incidence due to vertical motion 
of the aeroplane does not produce an alleviating air 
load immediately. 

(2) The flexural response of the wing is aerody- 
samically damped, and this is subject to lag. 

(3) Dynamic response in the torsional modes causes 
wing twist, which modifies the lift. 

(4) The 
transient changes in lift modify the lift on the tail 


free vortex sheets associated with the 
plane. 

These and many other interactions must be consid 
ered when a rigorous analysis of the behavior of an 
aeroplane in a gust is undertaken. 

In view of the complicated nature of the aeroplane 
response, it is not surprising that reliable evidence on 
gust shapes and velocities is lacking. Two approaches 
to the design methods can be employed. 

In the first case, let us suppose that a certain test aero 
plane has accumulated sufficient I’-g records to estab 
lish statistical data. It is then possible to specify 
that at a given speed a certain value of the normal ac 
celeration can be used as a design value to ensure 
safety. 

Now, if another aeroplane were being designed 
provided that the Wagner lag, aeroplane response, and 
dynamic effects were the same as those of the reference 
aeroplane it could be designed to an arbitrary sharp 
edged gust and the lag and response effects could be 
neglected. This approach is not adequate in all cases, 
since it can be demonstrated that the value oi the wing 
loading has an important effect on vertical response 
However, a simple correction can be made by analysing 
the response with different wing loadings in a sharp 
edged gust, and the alleviating factors thus obtained 
can be correlated to the approved standard. This 
procedure does not allow for differences in Wagner 
lag or in dynamic response, which may be important. 

In the second case, the analysis is more detailed. 
Returning to the test aeroplane, if the I’-g recorder is so 
placed that structural accelerations due to dynamic re 
sponse are negligible, a calculation can be made for a 
specified gust in which aeroplane response and aero 
dynamic lags are allowed for. The question is then, 
what should the gust shape and value be? The evi 
lence published so far is not convincing. 

British design practice is to assume a vertical velocity 
of gust of 50 ft. per sec. E.A.S. at cruising speed. A 
range of gradient distances is considered, and the maxi 
mum velocity is sustained. The assumed gust shape 
must be used in the response checks of the reference 
aircraft. If a different shape of gust is taken, a dif 
lerent velocity must be used to produce the same ac 
Such negligible 


differences to the vertical accelerations to be used in 


celeration variations will produce 


design, 
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The design cruising speed for a civil aeroplane is cur- 
rently fixed at a somewhat arbitrary figure, which is 
taken as constant E.A.S. with altitude, the limiting 
altitude being that at which this speed is equal to the 
maximum potential cruising performance of the aero 
thereafter it falls with increasing altitude be- 


plane; 
Furthermore, the vertical 


cause Oi engine power limits. 
gust velocity (E.A.S.) is reduced above 25,000 ft. in the 
current ICAO requirements. Consequently, the criti- 
cal altitude for gust stressing is seldom higher than this 
figure. For constant E.A.S. cruising and gust speeds, 
the intensity of gust loading increases with altitude be 
cause of the increase in compressibility factors. In the 
case of an aircraft having I’, equal to 300 m.p.h. 
E.A.S.) up to 25,000 ft., the increase of gust intensity 
due to compressibility may be 20 per cent above the 
sea-level value. 

Statistically, 
metrically more often than not. 


an aeroplane will enter a gust asym 
The colli 


sion produces the highest accelerations. 


head-on 

If a rolling 
motion is set up by an oblique intercept, although the 
downward moving wing aggravates the conditions 
when it enters the gust, the overall effects are com 


pensated by the effectively greater gradient distance. 


Application in Design - 

These methods of detailed gust analysis, though labo 
rious to apply, give the designer a clearer indication 
of the characteristics of a particular case and may offer 
scope for alleviating some of the worst effects arising 
from gust loads on flexible high-speed aeroplanes. 

Such calculations the and landing cases 


have been carried out by The Bristol Aeroplane Com 


for gust 
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PEAK ACCELERATIONS IN WINGS DETERMINED BY 
MODAL DYNAMIC RESPONSE CALCULATIONS 





50 FT/SEC. UPGUST WITH GUST ALLEVIATOR. 
UNFACTORED. 
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PEAK ACCELERATIONS IN WINGS DETERMINED BY 
MODAL DYNAMIC RESPONSE CALCULATIONS. 





12 FT/SEC. SYMMETRIC LANDING, 0-0-8. 
UNFACTORED. 

Fic. 9. 
pany, Ltd., on all types since 1944. The following gen- 
eral effects appear from the interactions between dy- 
namic response and damping: 

Fig. 7 shows the response of an inertia-elastic system 
having one degree of freedom to an applied load whose 
The 
The static 


time history resembles the assumed gust shape. 
system is first assumed to be undamped. 
response is that which would accrue if the system had no 
mass; it represents the asymptotic condition where 
The maximum 
dynamic response is the peak value of the actual re- 


the natural frequency tends to infinity. 
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sponse, which does not necessarily occur at the sam, 
instant as the maximum applied load. It is clear that 
for loads whose time of application is less than th 


le 


natural period of vibration of the system, the amplity, 
can be up to 100 per cent greater than would be foun; 
by neglecting dynamic effects. 

In fact, the aeroplane modes excited in a gust yj 
have aerodynamic damping, and this damping is mog 
important for the higher modes. The results so moj. 
fied are also shown in Fig. 7. The indication is that a 
important simplifying assumption can be made; th 
dynamic overswing of modes having 7°/) greater than? 
can be neglected. 

The procedure then is to treat the problem as if jt 
were static and then add the difference between the 
dynamic and static responses of the lower modes 
vibration of the aircraft. 

Figs. 8 and 9 show the maximum transient acceler 
tions in the wing of a particular aircraft due to its dy- 
namic response to gust and landing cases. It is note. 
worthy that the transient accelerations can be as much 
as Six or seven times the steady values; this large mag 
nification is due to the interactions of the large number 
of degrees of freedom. 

At any given altitude the higher cruising speeds as. 
sociated with turbine engines will increase the severity 
of the aerodynamic loads due to the gust. There is, 
however, an important alleviating effect, in that the 
aerodynamic damping for wing flexure increases as 
forward speed increases, and this reduces the total 
response in the fundamental mode. In fact, at high 
speeds it may make the system ‘‘dead-beat’’ or better 

This again is similar to a simple spring mass system 
in which the load rises rapidly to a maximum and then 
falls off in intensity. If sufficient damping of the 
system be present, then the maximum amplitude and 
therefore the stress induced in the system can be les 
than would be induced by the same maximum load 
statically applied; started otherwise, the maximum 
stress induced can be less than that calculated assuming 
the system rigid. 

For an aeroplane the damping depends on the dy- 
namic head, the compressibility factor, and the ¢- 
fective change of incidence due to vertical velocity 0 
vibration. It can be found that damping will decrease 
in importance with increasing altitude and _ increase 
rapidly with increasing speed. 

These points are illustrated in Fig. 10, which shows 
the results of calculations carried out on an aeroplane 
of 230-ft. span and 300,000-lb. gross weight. The 
results were obtained by step-by-step methods taking 
into account Wagner effect, freedom of the aeroplane 
to move vertically and to pitch, aerodynamic damping, 
and stiffness but neglecting effects. 
The conditions are for a gust of 50 ft. per sec. E.A.S. 


“virtual inertia” 


vertical velocity and a gradient distance of 100 ft. 
The points of interest are: (1) The general shape 


of the curves for the aeroplane are similar to those ol 4 
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(2) Since each of the curves is pro- 


simple system. 
portional to the applied loads of its own case, the more 


rapid fall-away after reaching the maximum illus- 
trates the greater importance of damping at higher 
speeds. (3) At 370 m.p.h. the flexible aeroplane 
actually experiences less stress in the gust than would 
be the case if the gust were considered a static problem. 

Hence, the highest speed associated with a given gust 
intensity will not necessarily give the worst design con- 
ditions, and it would appear that, as cruising speeds for 
large aeroplanes rise, the importance of the “‘gust at 
cruising speed” will diminish. 

It is interesting to investigate the general trends of 
the effects of varying size and speed on aeroplane re- 
sponse to gust conditions. Simple physical arguments 
can be used to demonstrate the effect of the variables. 

To give a picture of the quantitative effects involved, 
a family of similar aeroplanes differing only in size 
have been considered. These machines are repre- 
sentative of current practice, with spans from 70 to 
230 ft., and have the same plan form, mass distribution, 
wing loading, and strength in terms of ultimate normal 
acceleration. All are assumed of the same construction 
utilising the same maximum 
These assumptions enable the vibration characteristics 
of all members of the family to be expressed in terms 
of span only, and it may be mentioned that only the 


permissible stresses. 


fundamental mode of wing flexure is of primary impor- 
tance in this connection. It is reasonable to suppose 
that the results obtained for this family may be ap- 
plied without major error to most orthodox present- 
day designs. Calculations have been made which as- 
sume the currently required gradient distance of 100 
ft., invariable with altitude, and allowance has been 
made for the combined effect of aeroplane motion and 
aerodynamic lag in the buildup of the air loads upon the 
time history of the applied loads. To a close approxi- 
mation, this allowance amounts to increasing the gradi- 
ent distance by 20 per cent and assuming a sine varia- 
tion of loading with time in place of the linear. 

The effect of speed to a first approximation can be 
found by first obtaining the undamped velocity of 
vibration and then computing the equivalent damping 
forces. As the damping forces are proportional to the 
square of the speed, they increase in relation to the 
applied force as the speed rises. 

Fig. 11 shows a typical plot of overstress factor 
against speed at sea level, the dotted curve showing 
the value of the factor if zero damping is assumed. 
This curve shows clearly the importance of damping 
at all operating speeds and particularly at the higher 
speeds. 

These curves, unfortunately, show that the maximum 
overstress factor occurs at around the present-day 
values of the cruising speed. Similar plots for the 
Whole family at various altitudes have been made and 
used to construct Fig. 12, which shows the maximum 
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EFFECT OF FLEXIBILITY ON WING ROOT BENDING 
STRESS DUE TO PASSAGE THROUGH A GUST. 
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VARIATION OF OVERSTRESS FACTOR WITH SPEED 


I70FT. SPAN IOOFT. GRADIENT DISTANCE 
AT SEA LEVEL. 





OVERSTRESS FACTOR 
- 














overstress factors as a function of span, with altitude 
as parameter. 

The most striking thing about these curves is the 
rapid curving over at small values of the span, fol- 
lowed by a roughly constant factor for all values 
of the span from 100 to 250 ft., with, apparently, a 
drop in favour of very large spans indeed. There is, 
of course, a simple physical explanation as to why, 
contrary to earlier expectation, the very large span 
aeroplanes do not, in fact, experience much greater 
overstress than the (comparatively) smaller machines. 
At a given altitude and forward speed, the aerodynamic 
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VARIATION OF OVERSTRESS FACTOR WITH SPAN 
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damping forces are proportional to the vertical velocity 
of deformation of the wing. Now, in flying through the 
same gust, the large machine will experience deflections 
greater than the small one roughly in proportion to 
span—t.e., a 250-ft. span machine will experience ap- 
proximately two and one-half times the deflection of 
the 100-ft. span. 
takes place, however, is clearly not much greater than 
for the smaller span, so that the rate of deflection and, 


The time in which this deflection 


consequently, the damping force also are about twice 
that of the smaller aeroplane. Thus it comes about 
that what the large machine loses due to its slow fre- 
quency of oscillation it gains due to superior damping 
in flexure. 

altitude. 
The maximum overstress factor, about 22!5 per cent 


This figure also indicates the effect of 


at sea level, has become no less than 50 per cent at an 


altitude of 36,000 ft. An increase of the gradient 
distance to 200 ft. in place of 100 ft. would reduce this 
figure to approximately 30 per cent. 

The general conclusions that can be drawn from this 
study of dynamic response to a gust constant in shape 
and velocity are as follows: 

(1) At any given altitude, the overswing factor in- 
creases with increase of speed until aerodynamic 

With further in- 


crease of speed, the overswing factor gets smaller until 


damping overrides the gust effect. 


the wings become ‘“‘dead-beat”’ or better. 

(2) The maximum value of overswing factor is ap- 
parently independent of the size of aeroplane. 

(3) The speed at which the overswing factor reaches 
its maximum value decreases with increase of size of 
aeroplane. 

(4) The overswing factors increase with increase of 
altitude due to decreased damping. 


SCIENCES—FEBRUARY, 1950 

As the true magnitude of these dynamic effects js 
controlled by the rate at which a gust builds up, it js 
most important that we should learn more about acta) 
gust structures. 
that the gradient is proportional to the square of th 


If, for example, it could be establishes 


true velocity of the gust, overswing would be inde. 
pendent of altitude. 

In any case it appears that for high cruising speeds 
the gust cases are actually less severe on large aero. 
planes than on small ones, and on this count there is no 
limit to size. 


Gust Alleviation 


The detection of gusts at long range by radar depends 
on cloud formation. While this system appeared prom- 
ising at one time, recent flight experience has shown 
that severe disturbances occur in clear air and that 
long-range gust warning is not therefore possible. We 
must assume therefore that the aircraft cannot alter 
course to avoid a gust or reduce speed before encounter- 
ing it. In certain circumstances it may therefore be 
desirable to find means whereby the applied loads on 
the wings in gusts can be reduced so as to permit a re- 
duction of structure weight. There are two methods 
of achieving this. One is by so arranging the response 
of the aeroplane that alleviation takes place automati- 
cally; the other is by detecting and measuring some 
quantity associated with the gust (e.g., gust velocity 
and using it to operate a servomechanism that supplies 
the alleviation. 

An example of the first type of alleviation is to cause 
the wing to twist nose down as the load due to an up- 
gust builds up, thus reducing incidence at the wing 
tip. This entails either a forward distribution of the 
flexural and inertia axes or a skew hinge at the root 
The flutter and structural problems of these solutions 
for straight wings lead one to seek more feasible alter- 
natives. 

As an example of the second type of alleviation, a 
symmetric up-gust, which is usually the most severe 
case, might be alleviated by nose-down pitching of the 
whole aircraft. Analysis, however, shows that the 
angular accelerations are greater than normal elevator 
power can provide. While this approach might be 
made to work and produce some favourable effect on 
the structure, it would no doubt be balanced by equal 
but opposite reactions of the passengers. 

By far the most attractive solution makes use of 
symmetric motion of the ailérons, either by direct me- 
chanical coupling or by a servomechanism. This sys- 
tem, while potentially capable of reducing the bending 
moment at the wing root by 50 per cent and of almost 
completely eliminating the bending on the outer wing, 
does not involve unacceptable large pitching or vertical 
motion of the aeroplane. 

Many types of gust detection and measurement have 
been suggested. One proposal makes use of wing 
strain. 
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TRANSPORTS 


LARGE 


The difference in extension between a free rod and the 
tension surface of the wing across the span is used to 
operate the ailerons.'© Two methods of operating the 
ailerons are possible; in one, the rod operates the 
ailerons directly. The system of using wing strain is 
mechanically simple but suffers from appreciable time 
lags. The peak dynamic response of the wing occurs 
after the maximum gust velocity is reached. This 
response activates the ailerons whose effect is subject 
to Wagner lag, and there is a still further delay before 
the change in lift produces the dynamic response in the 
alleviating sense. This means that in practice the 
aileron is out of phase with the gust by 0.1 sec. or longer 
for a large aeroplane. This greatly reduces the effec- 
tiveness of alleviation, and the implications as regards 
flutter must be closely watched. 

As an example of the flutter implications, consider 
the wing vibrating in its flexural fundamental mode. 
When the tip reaches its extreme “up”’ 
strain between the tension surface and the free rod will 
Owing to the lags 


position, the 


demand an “up” aileron angle. 
in the system, this lift due to aileron will not act until 
the wing is on its downward beat; thus, the downward 
lift due to aileron acts as negative flexural damping 
and may be sufficient to reduce the total damping to 
zero. 

The vertical acceleration of the aircraft is a good 
measure of the gust velocity, but alleviation systems 
based on this quantity suffer from even greater lags 
than those using wing strain and must be rejected on 
that account. 

A device ahead of the wing which measures directly 
the change of incidence due to a gust is probably the 
most promising method, operation of the ailerons 
being by some power or servo system. The essential 
need is to keep the time lags for the set of operations 
within the small margins permissible so that the allevia- 
tion is synchronised with the arrival of the gust. Ina 
large aeroplane by virtue of a long body extending for- 
ward of the wings some advanced warning is possible; 
in the Brabazon I, a distance SO ft. forward gives a 
time of 0.2 sec. between the arrival of the gust at the 
nose and at the wing. Also, the dynamic response pro- 
vides some smoothing and permits an additional lag of 
the order of 15 per cent of the wing fundamental 
period. 

Even so, with large power-operated controls, these 
times still leave the problem a difficult one; the ef- 
ficiency of alleviation falls off rapidly with increasing 
lags. Fig. 13 shows the distribution of lags in the re- 
sponse to the British variety of gust for the Brabazon I, 
the total lag being just acceptable. 

A device of this nature needs to be automatically 
compensated to deal with changes of incidence due to 
changes of flight speed. As static and manoeuvre mar- 
gins can also be adversely affected, the system should be 
insensitive to minor disturbances. 

Another precaution that must be observed is in the 
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placing of the measuring meter, for, if it has an ampli- 
tude in any mode of the same sign as the amplitude at 
mid-aileron, the damping in that mode will be reduced. 
Unlike the previous case of alleviation from wing 
strain, this effect is not primarily due to the lags in the 
system. If we consider the case of a detector placed 
just forward of the aileron, when the wing is passing 
through its zero position in the upward direction, the 
device will detect a fictitious down-gust and demand 
“down” movement of the ailerons. This will cause a 
lift on the wing acting in the sense of negative flexural 
damping. The effect is most marked in the case of the 
fundamental mode, where the total aerodynamic damp- 
ing can be reduced by as much as 50 per cent; this 
could easily produce flutter at low speed. 

The effect can be eliminated by placing the detector 
at a position where the amplitude of the lower modes 
relative to the mid-aileron are small and negative. 
The effect on the higher modes can be neglected. 

This elementary discussion of the lag and stability 
characteristics of gust alleviating systems is intended 
only to make clear some of the physical factors involved. 
The system as a whole forms a servo loop; the various 
contributory lags, both aerodynamic and mechanical, 
must be separately analyzed and the stability of the 
whole servo loop studied for a suitable range of fre 
quencies before any particular system can be claimed 
to be safe. In general, however, the systems with low 
lags and small feedback will present the less severe 
problems and are therefore to be preferred. 

Flutter 

The problem of flutter is one of the earliest associated 
with flying but has, until comparatively recent times, 
been solved merely by solving the problem of strength 
coupled with the tactical distribution of lead weights. 
We are now becoming proficient at solving the prob- 
lem of strength and are disturbed if our test structures 
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do not fail within 1 or 2 per cent of the design loads 
(however arbitrary these loads may be). At the same 
time, this steady refinement of design has resulted in a 
reduction of structure weight for given design loads. 
Refinement of structural design has reduced the stiff- 
ness of the structure, and this, coupled with a steady 
increase of cruising speeds, has brought the flutter 
problem into its present prominence. 

Flutter is by no means confined to large aeroplanes, 
although at one time it was considered that flutter of 
such aeroplanes would involve considerably more de- 
tailed analysis than for small aeroplanes. However, it 
now appears that it is not so much the size of the aero- 
plane which matters but the type of structure employed 
and, in fact, that any new aeroplane designed to 
modern strength theories will require a complete flutter 
analysis. 

Any part of the aircraft structure where a structural 
deflection results in a change of aerodynamic load on 
the part is liable to flutter. Thus arise the familiar 
problems of wing flexure-torsion, fuselage-tail plane, 
and fuselage-fin flutter. Added to these we have fur- 
ther degrees of freedom due to the motion of a dozen 
separate analyses which may have to be undertaken 
before the aeroplane flies—in fact, some considerable 
time before the aeroplane flies. 
likely to be designed by the flutter requirements, these 
analyses must all be completed in the early stages of the 


If the aeroplane is 


design. 

The criteria available are of limited scope; the wing 
flexure torsion flutter criterion!” '* and the spring and 
servo tab criterion’ *° are examples and cannot be ap- 
plied to unusual weight or structure distribution. Since 
most new aeroplanes are unusual in one or other of 
these respects, the criteria cannot be taken as more than 
rough guides. Thus the prospect of having to design 
an aeroplane on the basis of flutter calculations is alarm- 
ing. 

It seems improbable 
parameters they include, can be developed to cover the 
multitude of inertia and elastic systems present in 
The only answer seems 


that criteria, however many 


different types of aeroplanes. 
a complete calculation for each variation of mass, 
stiffness, or geometry. The major difficulty to be 
tackled in this field seems, therefore, to be excessive 
time required for each flutter calculation. 

Another approach to this problem can be made by 
introducing the flutter model either in a wind tunnel 
or in free flight. That there are many difficulties 
is immediately apparent when one considers the task of 
producing a model that not only has the requisite 
elastic and inertia properties but also must have the 
correct aerodynamic properties. Because of scale ef- 
fects it is likely that the model would of necessity be 
large, ruling out the wind tunnel for all except local 
flutter effects. Flutter model testing may therefore be 
limited to those cases where the aerodynamic deriva- 
tives cannot be determined with sufficient accuracy for 
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calculations. This is not to suggest that the experi. 
mental approach is useless—it is essential that our 
aerodynamic derivative theories be checked by direct 
experiment—but that this approach is unsuitable for 
ad hoc testing of particular aeroplane designs and js 
best suited to the verification of the basic flutter data. 

There remain two (very similar) approaches that we 
may make to this problem and which both show 
promise. These are the electric analoge and the electric 
computing machine. 

The electric analoge is a development of the well- 
known analogy between the equations determining the 
oscillation of a structure and the equations of the os- 
cillation of an electric circuit. In this, the inductance of 
the circuit replaces the inertia, the resistance replaces the 
damping, and the capacity replaces the elasticity of 
both structure and air. Thus a circuit may be built up 
which will oscillate at a certain frequency defining an 
equivalent flutter speed and mode. Once constructed to 
simulate any particular type of flutter problem, the 
answer would be given immediately and the effect of 
variation of structural and aerodynamic parameters 
easily found. 

The electronic computing machine is a device for 
applying a specified programme of operations to a given 
set of parameters in a short time. Thus the normal 
methods of flutter calculations are amenable to treat- 
ment by this machine, as, in fact, are all other calcula- 
tions that require a definite sequence of computations. 
Once again, the short time required for each complete 
calculation enables a greater variety of structural and 
aerodynamic parameters to be considered in the time 
available. 

Both devices suffer from a high initial cost, and it is 
in this respect that the electronic computing machine, 
with the advantage of its greater adaptability to all 
types of design calculation, may be effectively cheaper 
to run. 


Wing Flexure-Torsion Flutter 


Although it is impossible to obtain a really accurate 
value for the flutter speed for any general analysis rely- 
ing on a criterion and fairly broad assumptions, it 1s 
often instructive to observe the general trend as vari- 
ous structural and aerodynamic parameters are altered. 
The wing flexure-torsion flutter criterion!’ (i.e., rigid 
aileron) is a reasonable guide to the flutter speed and, 
although it neglects the effect of variations of mass dis- 
tribution (e.g., wing engines) and stiffness distribution 
between different aeroplanes, may be assumed true 
for this general investigation. 

The criterion may be stated as 


Vec = K, (1 — M*)* (mO/s,,,”) (1) 


where J\/ is the Mach number corresponding to the 


flutter speed, Vc, and where a minimum value ol 
(1 — M*)’* = 0.775 is applied for all values of Mf > 


0.8. 
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LARGE TRANSPORTS 

K, is a function of wing loading, taper ratio, and 
inertia axis position, although the variation with wing 
loading is extremely small and may be neglected. 

It is now required to determine the effect on mé 
and, hence, |’»c of various parameters such as wing 
loading (w) aspect ratio (A), spar box dimensions, etc., 
or typical wing structures. 

By assuming a certain inertia factor (m) and inertia 
relief factor (7,), it is possible to design a wing struc- 
ture and obtain general expressions for skin and spar 
web thicknesses in terms of the various parameters. 
From this thickness distribution a value of m may be 
determined. This analysis results in the following ex- 
pression for Vee: 


= 10°A, (1 — M*)" X 


Vee 
nw (1 


Vi fA LK £(1 4 =f (2) 
*? py? ” k; bt 


= Br} 


where A, and A; are constants being about 0.002 and 
0,008, respectively, while A, is about 65 (where ft.lIb. 
sec. units are employed). 

An idea may be obtained from this expression of the 
effect of various structural parameters on the flutter 
speed. The effect of aspect ratio (A) and of wing load- 
ing (w) is shown in Fig. 14, where it is evident that, 
although increase of aspect ratio reduces the flutter 
speed, this may be offset by an increase of wing loading. 
Thus it appears that, if high cruising speeds are de- 
sired for aeroplanes of the future, an increase of wing 
loading or a reduction of aspect ratio will be inevitable 
unless the criterion of wing design is allowed to change 
from strength to stiffness. It must be remembered 
that a wing designed for strength is inefficient as a tor- 
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DATA: 


sion member, the skin thickness decreasing from root 
to tip as the first or second power of the span, whereas 


on a wing designed for torsional stiffness only the 


thickness should remain 
wing tapers towards the tip. 


constant or increase if the 


If the stiffness becomes a 


design condition, it is therefore likely to modify our 
ideas of the optimum form of compression structure 


for the wing box and will probably result in a smaller 


ratio of stringer area (As) to skin area (bt) than has 


hi 


progressively reduced towards the wing tip. 


therto been considered desirable, this ratio being 
Fig. 15 


shows the effect of As/bt and of design compressive 
stress (fc) on the flutter speed, assuming, however, that 


As/bt is constant along the span. 


A method often proposed for increasing the torsional 


stiffness and the flutter speed of a wing is to increase 


the interspar chord (k:). 


This generally has an un- 


favourable effect on the bending efficiency of the box, 
but the effect on flutter speed is marked, as may be seen 


from Fig. 16. 


ae 
i.e 
ve 
values for the constants A,, Ko, and K3. 


It is of interest to note here that the other major 
roelastic effect that fixes the torsional stiffness— 
., aileron reversal—tresults in an equation for re- 
rsal speed similar to Eq. (2) above but with different 
The aileron 


reversal speed as a function of wing loading and aspect 
ratio is plotted on Fig. 17 and in all cases is higher than 
the flutter. 


The above analysis indicates therefore that the 


adoption of wing loading of about 80 Ibs. per sq.ft. 
as a method of reducing wing structure weight (based 
purely on bending considerations) is favourable from 


the aeroelastic point of view and may allow some in- 


crease of aspect ratio (provided, of course, this does 
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not seriously increase the structure weight required for 
bending strength). However, as has been mentioned 
previously, no account has been taken of the effect of 
It is 
impossible to consider these effects in detail in a general 
analysis such as this, although it may be said that en- 
gines placed forward of the flexural axis usually result 
There are possibilities 


engine and other masses on the flutter speeds. 


in a raising of the flutter speed. 
here of reducing the required torsional stiffness by a 
favourable distribution of masses (e.g., engines placed 
at the wing tip), but much more work is required be- 
fore such an effect may be included with confidence in 
the torsional stiffness criterion. 


Wing Analysis 

The following features that have a bearing on wing 
design have been considered: (1) the effect of a varia- 
tion of wing loading and aspect ratio on the perform- 
ance of an aeroplane and on the weight of fuel required; 
(2) an examination of the dynamic response of an aero- 
plane in the gust case and the consequent effect on the 
stresses and the structure; and (3) the provision of 
structural stiffness necessary to avoid wing flutter and 
aileron reversal. 

The task now is to determine the layout of a wing 
that will satisfy airworthiness requirements and permit 
the highest possible efficiency in operation. This in- 
volves the selection of wing loading and aspect ratio in- 
volving a balance between the weight of the structure 
satisfying the demands of strength and stiffness and 
the weight of the fuel appropriate to the particular 
performance characteristics. For this purpose some 
rule must be established whereby the effect of these 
variables on structure weight can be assessed. The 
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L950 


influence of size must also be included. Those who ; 
the past based their forecasts on the square cube |gy 
and established a limit to the maximum possible size ¢j 
aeroplane must by now have diminished confideng 
On the other hand, those who entirely ignore the effeg 
or rely on some mystic solution can introduce larg 
errors into their overall evaluation. 

Methods of weight analysis fall into four separat 
categories successively approaching an increased orde 
of accuracy. These stages are appropriate to: 

(1) A general survey for the purpose of suggesting 
probable trends and comparing the relative merits oj 
different forms of layout and types of power plant. 

(2) A particular design study of an aeroplane whos 
form and type of power plant have already been decided 
as a result of a general survey. In this case the de 
signer is using as a basis an approximate specification 
of user requirements, and his purpose is to determin 
the probable gross weight and the optimum combina 
tion of span, wing loading, and other characteris 
tics so as to obtain the best possible performance and 
economy of operation. 

(3) In the detail design stage when weight estimates 
can be made from production drawings. 

(4) In the manufactured state when actual weigh 
ings can be made of details, components, and the com 
plete aeroplane. 

For the first two stages one is forced to rely ona 
semiempirical cum statistical analysis, taking the results 
of known aeroplane to establish general coefficients 
Care is needed in the selection of the data for this pur- 
pose, since each aeroplane is subject to many variables 
and compromises, often known only to the designers 
so it is important to separate the various controlling 


factors. It is likely that one would find a variation 
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LARGE 


in the general terms of formula used for weight esti- 
mation as derived by different technicians, but no 
doubt each design team feels that its own methods have 
special merit. In novel circumstances the first esti- 
mate may be an inspired guess, but this should be fol- 
lowed at the first opportunity by a more fundamental 
approach. 

‘It is not my intention to describe in detail any par- 
ticular form of estimate, but, so that an argument can 
be pursued, it is desirable to outline the effects of 
relevant factors on wing structure weight. 


A wing can be split up into different parts each ful- 
filling a main function—viz., the compression struc- 
ture and ribs, the tension structure, the shear members, 


leading and trailing edges, flaps, and controls. 


Taking the compression structure and ribs, it has 


been demonstrated that the instability failing stress f 


of a typical panel is given by the equation 
f=KVPE/I 
where 


P = compressive end load per in. width 


lI 


panel pin centre length 
a constant depending on the geometry of the 


ll 


l 

K 
stiffeners 

The value of A is a measure of the radius of gyration 

per inch width of the panel and can be used to compare 

the relative efficiencies of different forms of stiffener. 


This formula gives a good indication of the layout of 
wing structure that should be adopted to obtain mini- 
mum weight. At a chordwise section at or near the 
wing root, the instability failing stress f should be 
equal to the yield stress of the best material available. 
Ifa form of wing skin stiffening having a high value of 
K is selected, the desired stress can be achieved by ar- 
ranging the wing geometry to give the required ratio of 
P/l to satisfy the equation. That is to say that, by a 
method of successive approximation the width and 
depth of the main structural box at the root is made to 
control the value of the compressive end load per inch 
resulting from the spanwise applied bending moment 
and the rib pitch fixes the panel pin centre length. 
Clearly, this optimum ratio of P//] should be main- 
tanned as far as possible outwards along the span. 
To assist in this aim the wing should have a straight 
taper. Also, resistance to the applied bending moments 
should be confined to a main structural box, and con- 
tributions from other secondary portions such as the 
leading and trailing edges should be kept as small as 
This could be achieved by providing a 
honstrain-carrying connection between the primary 
and secondary structures. The structure of the second- 
ary portions can then be made as light as possible, 
since the only forces it has to carry arise from the air 
These are transmitted directly 


possible. 


loads on its own cover. 
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to the main box. This is a conception that in practice 
is difficult to attain completely. 

The weight of the sum of the ribs and the compres- 
sion structure can now be estimated. Away from the 
wing root, the rib pitch plays an important part in 
maintaining an optimum value of the ratio P/I. 

If the weight of a typical rib is expressed in the 
form dt, where d = depth of rib, /, = equivalent thick 
ness of rib, and p = density of material—the weight per 
unit area of a typical portion of the compression struc- 
ture and ribs can be simplified into the form 


W,? = p [(dt,/l) + (1/K)V PL/E] 


The rib pitch can now be varied to make this weight 
a minimum by equating Ow,/0/ to zero so that 


W,! (min.) = p (dt,,P/E)'" [3/(4k?)]' 


P, the end load per in. width, can be expressed in 
terms of the load factor, wing load, wing span, and 


wing depth, so that 
W/W « (Nb? t,/w2)”* (1/K“’) 


It is now possible to see what saving can be had by an 
improvement in skin-stiffener design. The value of A 
appropriate to a well-designed Zee section is 0.5, 
whereas it is possible to achieve a value of 1.1 if a com- 
bination of inverted Y sections and minor intermedi- 
ates is used. It is probably optimistic to assume that 
this arrangement could be maintained all the way to the 
tip, since the section thicknesses would be limited by 
manufacturing processes; therefore, a mean value for 
K of 1.0 should be used. 
term containing A shows that the weight of the ribs 
and compression structure could be reduced by about 
13 per cent by a change from Z stiffeners to the more 


Inserting these values in the 


efficient combination. 

By comparing the component parts of suitable large 
aeroplanes, a coefficient can be found whereby weight 
estimates for other aeroplanes can be made. The 
part of the equation defining the compression structure 


plus rib weight is then 


W/W = 43 p (N,*t,/Ew?)'”? (1/K”) 


The remaining parts of the wing—i.e., the tension, 
shear, and secondary structures 
relationships with the load factor and wing geometry. 
Their weights can be expressed in a more simple form, 
and a series of coefficients can be found by equating 
them to the known weights of the similar components 


have more direct 


of existing aircraft. 

With a basis for estimating the weight of the struc- 
ture for different sizes of wings subjected to specified 
applied loads and taking into account the dynamic 
effects previously discussed, it is possible to draw up a 
general picture from which the weight of a particular 
wing can be assessed. 

The gust cases will usually provide the design strength 
criteria for the wings; typical load factors for various 
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wing loadings and cruising speeds are shown in Fig. 18. 
These curves apply at altitudes below 25,000 ft., since 
some relaxation in the gust indicated velocity is per- 
mitted at greater altitudes. As the performance of 
aeroplanes fitted with turboprop engines is such that 
the maximum indicated cruising speed increases with 
decreasing altitude, to avoid excessive load factors it 
may be necessary to limit the maximum permissible 
indicated speed below some specified altitude. This 
should not be a handicap in normal operation. 

The factors shown in Fig. 18 should not be used in 
the formula for weight estimates without making an 
allowance for overswing. Following the discussion 
where it was shown that, while a constant overswing 
stress factor of 1.2 may be taken at a speed of 375 m.p.h. 
at ground level for any wing span between 100 and 250 
{t., the factor increases with increasing altitude, the 
reason being the reduced effect of aerodynamic damp- 
ing. While this may be so, it is unlikely that such un- 
favourable conditions would be accepted in the design 
of an aeroplane when alleviation is possible by the 
introduction of some suitable device. 

With these assumptions the wing structure weights 
shown in Fig. 19 have been produced. Two sets of 
curves are mae the one for which the ultimate 
strength factor is 3.75 and the other for the appropriate 
gust factor at a true speed of 375 m.p.h. at 16,000 ft., 
with a representative overswing factor of 1.2, and a 
factor of safety of 1.5. The difference between these 
two sets of curves shows the effective saving to be had 
by efficient gust alleviation. 


Obviously, Bristol aeroplanes fall into their a 
propriate place among the curves as the coefficieny 
in their equations were derived from them. 


Optimum Wing Loading and Aspect Ratio 

At this stage it is convenient to combine the resyl 
of Fig. 19 with those of Tables 3 and 4 to form an opin 
ion of the influence of wing loading and aspect ratio py 
operating economy (see Table 5) 


TABLE 5 


Fuel for 
Wing Aspect 5,000 Miles Wing Structure Fuel + Win 


Loading Ratio Range Weight Weight 
45 10 0.333 W 0.132W 0.4651 
55 10 0.348 0.123 0.471 
65 10 0.362 (a) )0.117 0.479 
75 10 0.376 0.112 (0). 488 
45 10 0.333 0.160 (). 493 
55 10 0.348 J0.140 (). 488 
65 10 0.362 (b) )0.125 0). 487 
46 10 0.376 0.115 0.491 
55 8 0.388 0.111 0). 499 
55 10 0.348 (a) <0.123 0.471 
55 12 0.309 0.135 0.444 


(a) With gust alleviation 
(b) Without gust alleviation. 


The following conclusions can be drawn: 

(1) With no gust alleviation and with an aspect 
ratio of 10, the combined weight of fuel and wing struc. 
ture weight is sensibly constant over a wide range oi 
wing loadings, and, hence, there is no optimum value 
to give maximum pay load. This obviously encourages 
the adoption of as high a wing loading as possible so as 
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to produce the smallest possible size and the highest 
possible cruising speed. An upper limit may be set by 
any of the following reasons: 


A minimum acceptable altitude after an en- 


(a) 
gine has failed so as to be able to fly above 
rough weather. This is only likely to be 
critical in the tropics. 

(b) The need to be able to take-off from short 
runways at airports situated in the tropics 
and at considerable altitude. 

(c) A volume limitation to contain the fuel. 


With jet pipes passing through the wing, 
landing-gear stowage space, heated air pas- 
sages in the leading edge for anti-icing, and the 
rear portion occupied by flaps and_ their 
operating mechanisms, the wing area can, in 
the case of long-range aircraft, be determined 
by the fuel capacity. 


(2) When gust alleviation is incorporated, a greater 
pay load can be carried if a moderate wing loading is 
selected. The difference between 55 Ib. per sq.ft. and 
75 Ib. per sq.ft. wing loading is equivalent to a differ- 
ence of pay load equal to nearly 2 per cent of the gross 
weight. 

(3) The aspect ratio should be as high as possible. 
Of the aerodynamic variables, this characteristic shows 
the highest profit. With a high aspect ratio the aero- 
elastic criteria are more easily met if the wing loading 
is also high. Even at moderate wing loadings there 
seems to be no reason why values in excess of 12 or even 
15 should not be adopted. 


Remainder of Structure 

Fuselage.—An analysis of the effects of different 
forms of fuselage is complicated by many interrelated 
operating conditions. If an aeroplane is designed to 
operate on one specific route, only then can the fuse- 
lage size and passenger accommodation be arranged to 
provide facilities just in accordance with the available 
capacity pay load. 

If a long-range aeroplane is required to carry greater 
pay loads at shorter distances, then its efficiency at 
long range may be seriously impaired. 

Suppose, for example, that a civil transport aero- 
plane with a range of about 5,500 miles has the most 
simple form of single-deck circular fuselage, having a 
floor area just sufficient to accommodate the appro- 
priate number of passengers. Then, if a considerably 
greater number of passengers must be carried at lesser 
double-deck or ‘‘double-bubble’’ fuselage 
must be provided. Assuming that the gross weight 
remains the same, the following changes take place: 
(1) a small increase of drag due to greater surface area; 


ranges, a 


(2) the greater surface and floor areas, together with the 
additional windows and doors, increase the structure 
weight; (3) similarly, the weight of furnishings and 
soundproofing are increased; (4) heavier air-condition- 
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ing plant and ducting; (5) heavier electrical installa- 
tions. 

The sum of these changes may raise the weight of the 
complete furnished and equipped body from about 
14 to 17 Assuming, 
too, that the maximum pay load is increased from 9 
per cent to 15 per cent, the total increase of loaded body 
Now at constant gross 


IS per cent of the gross weight. 


weight is about 10 per cent. 
weight there must be an equivalent reduction of fuel 
weight, thus reducing the relieving loads on the wing 
and thereby increasing the bending moments. The 
consequent increase of the wing structure weight is 
about | per cent. The landing weight is also greater, 
so that those parts designed by landing loads will also 
suffer. 

On this reckoning, if this double-deck fuselage is 
adopted for a long-range aeroplane, the basic opera- 
tional weight may be 4-5 per cent greater than would be 
the case for a circular fuselage. After making allow- 
ances for the remaining structural parts, items of equip- 
ment, crew, and passenger services, the two arrange- 
ments are compared in Table 6. The fuel for range is 
taken from Fig. 4. 


TABLE 6 
Comparison of Single- and Double-Deck Fuselage Nett 
Load-Range 


Pay 


Double Deck 


Single Deck 
Volumetric pay load 9% 15% 
5,500-mile range 9% 4.5% 
1,800-mile range ~Nett pay load 9% 9.0% 
3,900-mile range 9% 15.0% 


In this particular example, at a fixed all-up weight of 
aeroplane, the arrangement, with the body designed 
specifically to match the pay load at 5,500-mile range, 
shows at this range an advantage in the ratio of 2:1 over 
the arrangement with the body designed to match the 
pay load available at 3,900-mile range. At this lesser 
range the advantage is reversed. 

The only general conclusion I wish to draw from this 
example is that a penalty is inevitable if one tries to 
compromise long-range utilisation by excess capacity 
appropriate to lesser range. 

When the principal application for which an aero- 
plane is intended permits a greater load to be carried, 
then the double deck may be the most economic ar- 


rangement. 


Runways and Landing Gear 


It has frequently been asserted that for large land- 
planes the cost and constructional problems in pro- 
viding runways of adequate length and strength would 
become prohibitive. This need not in fact be the case, 
since it can be shown that multiwheel landing gear will 
more than compensate the effect of the increased weight 
upon the runway by suitably spreading the load. 

Several methods of spreading the load upon the run- 


way are open to selection. The aeroplane may have a 
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number of widely spaced main units each carrying one 
or two tyres, or it may have few main units with groups 
of wheels. Track units are a third choice but are 
fundamentally unsuited to the prepared runways under 
consideration. 

Before the merits of various forms of landing gear 
are assessed, it will be necessary to consider the general 
problem of runway bearing capacity. Runways and 
perimeter tracks can be classified broadly as ‘‘flexible 
pavements’ and “rigid pavements.’’ In the former 
group the subsoil is compacted and covered with a 
surface of tarmacadam, asphalt, or similar material 
that protects the subsoil from moisture and from break- 
up in dry weather, and, if the subsoil is properly graded, 
this construction can withstand considerable ‘contact 
pressures without appreciable settlement. The “‘rigid 
pavements’ are essentially beam systems capable of 
distributing wheel loads over an area appreciably 
larger than the tyre contact area. Concrete or rein- 
forced concrete slabs upon consolidated subgrade are 
in this group. 

“Flexible pavements’’ suffer damage by indentation 
reaching a permanent critical deflection. ‘Rigid 
pavements” characteristically suffer slab corner break- 
ages and sometimes transverse failures due to succes- 
sive hoggings and saggings under rolling loads. Fail- 
ures are aggravated by water penetration and erosion 
or differential settlement of the subsoil. The character- 
istic behaviour of the two different groups of pavements 
has led to differing approaches for the assessment of 
the loading effects on runways. The one is related to 
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the effect of a load upon the deflection of a continuoys 
flexible pavement, and the other is dependent upon the 
stresses in the corner of a discontinuous slab pavemen 
laid upon prepared subsoil. 

From theory and experiments in soil mechanics 
numerous workers have shown that the capacity of g 
runway structure to sustain a load IW depends on the 
contact pressure p and that Wp? is a constant. Ry. 
perimental values for the indices, based upon tests 
on both “rigid” and ‘‘flexible’’ pavements, have pro. 
vided the relationship of 


W°-5 49-44 = constant 


This agrees closely with the theoretical conclusions of 
Rodin’s deflection theory,*' which gives 


W°* p** = a constant 


As a pair of wheels and tyres are moved apart, the 
carrying capacity of the runway as limited by the 
critical deflection increases, since, by the method of 
superposition, the deflection beneath one wheel due to 
both can be shown to diminish. Average dual wheels 
as now in use are spaced about two tyre widths apart 
and permit an increase of about 50 per cent over the 
load carried on a single wheel with the same tyre pres- 
sure.*! When the dual wheels are about 4 tyre widths 
apart, the gain is 75 per cent and is about as great a 
gain as is economical with dual wheels. 

If wheels are spaced in tandem pairs instead of side- 
by-side, the reduction in severity of runway loading is 
comparable with that obtained by side by side. The 
use of larger numbers of wheels with tyre contacts 
spaced both laterally and fore and aft will in most cases 
give a satisfactory reduction in runway loading. A 
suitable arrangement is the four-wheeled bogie, which 
can be about two and one-half times less severe on the 
runway than an alternative single wheel. This saving 
and most of the foregoing points are illustrated by Fig. 
20. In this diagram the curved lines give constant 
runway loading and are drawn for the strengths of run- 
ways defined by the ICAO Runway Strength Classes, 
which are written in terms of unique combinations of 
wheel loads and tyre pressures.*? The reduction in 
tyre contact pressure necessary for the carriage of very 
heavy loads on a single wheel is seen to be prohibitive. 
The benefit derived from the use of dual wheels is 
shown by the reduction from single-wheel load to 
“equivalent single-wheel load’’ shown plotted at the 
same contact pressure for typical aircraft as Nos. 13, 
14, 15, 16, and 17. The load on one main wheel group 
of the Brabazon I is about 135,000 Ibs. This is re- 
duced to an equivalent of 100,000 Ibs. on the runway 
by the use of dual wheels with tyres at 20-in. centres. 

When the four-wheeled eight-tyred bogie, which 
replaces the dual twin-tyred wheels, is fitted, the 
equivalent effect is reduced to that of a load of only 
50,000 Ibs. at a contact pressure of 90 Ibs. per sq.in. 
In a similar way a smaller aeroplane, having a static 
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LARGE TRANSPORTS 
load per wheel group of about 56,000 Ibs., has an equiva- 
lent single-wheel load of about 36,000 Ibs. at 90 Ibs. per 
sq.in. with dual wheels (side by side) and an equiva- 
lent single-wheel load of only 22,000 Ibs. when a four- 
wheeled bogie is employed at a contact pressure of 110 
Ibs. per sq.in. 

The lines of constant runway strength enable ready 
comparisons to be drawn. It will be seen that the 
four-wheeled bogie at 110 Ib. per sq.in. tyre pressure 
requires a lower ICAO class of runway than the same 
aircraft with dual wheels at 90 Ibs. per sq.in. 

In the development of Rodin’s deflection theory of 
failure in reference 21, the natural outcome was to use 
Cwp as a Runway Loading Index, where C is a correc- 
tion factor for multiple wheels and tyre contact shapes, 
and on this basis an aircraft can be compared with an- 
other aircraft or with a runway strength even when the 
contact The ICAO strength 


classes have equivalent indices as follow in Table 7. 


pressure is different. 


TABLE 7 


Jearing Capacity for 


Isolated Wheel Load of Index 


Tyre Pressure Load 


Class (Lbs.) (Lbs. per Sq.In.) (Lb. per In.) 
| 100,000 120 3,460 
2 75,000 100 2,740 
3 60,000 100 2,450 
4 $5,000 100 2,120 
5 30,000 85 1,600 
6 15,000 70 1,025 
7 5,000 35 394 


If large aeroplanes are equipped with load-spreading 
landing gear, it can be expected that runways greater 
than ICAO Class 3 strength need not be provided for 
aeroplanes up to 300,000—400,000-Ib. gross weight. 

At least two factors conspire to restrict tyre pres- 
sures and sizes. High landing and take-off speeds will 
result in braking and tyre difficulties if the relative size 
of wheels is reduced as performance and weight in- 
crease. Also, a trend towards higher tyre pressure 
and surface speeds will increase the liability to tyre 
damage from stones and other objects and will demand 
an uncomfortably high standard of runway surface 
cleanliness. This situation is a deterrent in the use of 
very high pressures, although tyres at inflation pres- 
sures of 1SO Ibs. per sq.in. may be found to be satis- 
factory after suitable development. 

The most suitable form of landing gear for larger 
aeroplanes is one having groups of wheels. That this 
trend in development does not require an increase in the 
proportion of weight is supported by the evidence of 
Fig. 21, in which the percentage of the gross weight 
devoted to landing gear has been plotted for some aero- 
planes. It can be observed that no increase in per- 
centage weight need be expected with growth in size. 

One further aspect of the landing-gear arrangement 
deserves comment. Having demonstrated the de- 
sirability of load-spreading wheel groups, then loca- 
tion relative to the aeroplane structure is of some im- 
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portance, since in large aeroplanes with relatively 
flexible wings the dynamic effects of landing are con- 
siderably modified by the conditions of support. It is 
desirable to place the wheel group on each side at 
approximately the radius of gyration in roll. The nett 
result of this arrangement has proved lighter than alter- 
native arrangements of three and four groups of wheels 
without placing great demands on runway width and 


strength. 
CONCLUSIONS 


The results of this analysis of transport aeroplanes 
designed specifically for long-range operation and fitted 
with turboprop engines show that the wings are the 
only component of the structure where the proportion 
of weight inevitably increases with increase of gross 
weight. The actual amount of increase will depend on 
the values of wing loading and aspect ratio selected. 
If these values are maintained constant, the order of in 
crease is approximately 1'/»2 per cent’of the gross weight 
for each increase of 100,000 Ibs. in gross weight. The 
minimum basic value for reference should be that 
appropriate to an aeroplane of 100,000-Ib. gross weight. 

This means that for an aeroplane of, say, 300,000-Ib. 
gross weight it is possible to design the wings having a 
weight of 12 per cent and the total structure including 
the fuselage, tail unit, landing gear, flying control 
system, engine mountings, and cowlings having a weight 
equal to 28 per cent of the gross weight. 

It has been my purpose to discuss fuel, electric, 
hydraulic, and other equipment necessary for the 
functioning of the aeroplane and its engines, not opera- 
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tional and passenger services. This information can 
best be analysed by gathering evidence from estab- 
lished aeroplanes, and this procedure is normal prac- 
tice to all investigators interested in design trends. 
The relative weights of the various aeroplane services 
will vary from one installation to another according to 
local preference, but their sum total does show a rela- 
tive reduction with increase of gross weight. Also, it is 
fairly evident that the route to be flown will fix some of 
the operational equipment and that much of the re- 
mainder will not vary at the same rate as the gross 
weight. 

Again, as the proportion of pay load diminishes with 
increasing range, the weight of air conditioning and 
other passenger services, being largely dependent on the 
number of passengers carried, is also relatively smaller. 

Thus it is not surprising to find that in practice the 
ratio of the weight of a fully equipped and furnished 
structure decreases slowly with increasing gross 
weight. 

To claim that the installed weight of turboprop en- 
gines is also less than that for piston engines is true, 
but a claim on this basis that the percentage basic 
operational weight is also reduced, thus giving a higher 
disposable load, is fictitious if used as a comparison 
with piston-engined aeroplanes. In two _ particular 
cases it has been found that the sum of engine plus 
fuel weight is approximately constant at long range for 
alternative piston and turboprop engines, thus leaving 
a constant pay load. However, for long-range opera- 
tion the turboprop engine gives an advantage of ap- 
proximately 100 m.p.h. at most economic cruising speed. 

In this lecture I have confined myself to certain 
physical aspects of large transport aeroplanes. While 
this side is the major reason for the employment of 
scientists and technicians, it is but the key that unlocks 
the door to the debating chamber. Within, other im- 
portant questions must be asked before such aeroplanes 
can be justified. Certainly, the results of an estimate 
of economy in operation will be an overriding factor 
in the determination of policy. An authoritative as- 
sessment of this factor is due to Masefield,”* 
presenting the Fourth British Commonwealth & 
Empire Lecture in London, made a classic contribution 


who, in 


towards establishing data for use in the evaluation of 
the economics of civil aviation. Often a careful anal- 
ysis of all the items in the buildup of fixed annual 
costs, hourly cruising costs, and take-off and landing 
costs and making the necessary allowances for fuel 
reserves appropriate to normal civil operation, he 
derived the optimum stage length for a series of aero- 
planes covering a wide range of gross weight. Con- 
versely, he inferred the most economic size of aeroplane 
to suit specified stage distances. 

On the assumption that in each instance it was pos- 
sible to achieve 300 hours of annual utilisation with a 
65 per cent passenger load factor, he found that for 
3,000- and 3,500-mile stage length the optimum gross 
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weight was 400,000 and 400,000 Ibs., respectively 
Moreover, he showed that the lowest possible operating 
cost per revenue ton-mile does not start to rise appre. 
ciably until the stage length exceeded 3,500 miles, 

On a route such as the North Atlantic, where mediym. 
range aircraft are subject to frequent deviations jy 
adverse weather to off-route airports, an aeroplane with 
longer range will both fly a shorter overall distance ang 
take less time for the journey. 
distances the economic factor is favourable to larg 


Thus for long-stage 


aeroplanes. 

But in the end the customer is always right, for the 
psychology of the travelling public does much to deter. 
mine the passenger load factor when a choice is avai. 
able. 

With increasing size of body, more scope is possible 
in the layout of seating, lounges, and toilet arrange. 
ments, including possibly two or more standards oj 
comfort. One has a simple faith that passengers will 
prefer to have a little more space provided for their 
accommodation, and, while all will appreciate the 
lower noise and vibration levels of turbine engines, the 
few with an aeronautical upbringing will, in time, 
cease to speculate on the hazards of explosive decom- 
pression. 

Summing up the overall picture, when long range isa 
necessary or desirable requirement for civil air trans 
port, the role can only be filled by large aeroplanes 
The point has not yet been reached where the maximum 
operating efficiency is adversely affected by fundamental 
technical barriers. 
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Discussions of the Lecture 


R. M. Hazen (East Coast), Director of Engineering, 
Allison Division, General Motors Corporation: The 
stature of the Wright Brothers Lecture Series has been 
increased by Mr. Russell's fine contribution. His 
modest lecture gives an inkling of not only the great 
ability and effort put on the large transport problem 
by him and his associates but of the vision and courage 
required to see it through to completion. Having 
both the airplane and the turboprop engine under one 
roof, so to speak, has been unusual and perhaps has 
contributed to the well-rounded picture he has pre- 
sented. 

From the turboprop engine point of view one can 
search for important factors that have been missed 
but can conclude that any discussion must be confined 
to observations on points already mentioned. It 
might be observed in a large transport with four dual 
unit turboprop engines that there are really five alti- 
tudes at which the power plants can be running at or 
near their maximum efficiency by using 4, 5, 6, 7, or 8 
operating power sections and that these can be used 
without disquieting effect on the passengers since all 
propellers would be turning. While specific fuel con- 
sumptions would go up somewhat with reduced ex- 
pansion ratio at lower altitudes, it would appear that 
a number of actual operating conditions exist, such as 
head or tail winds, storm or icing conditions, cabin 
supercharging limits, standby and descent, which 
would find this advantageous in maintaining range 
ability on individual flights. 

Our experience indicates an extremely rapid rate of 
development of power output in turbine power plants 


of a given size as eompared to the reciprocating engine, 
due primarily to the unexplored possibilities of the 
There is continuous pressure for improved 
These 


former. 
specific fuel consumption on turbine engines. 
can be obtained by either or both component efficiency 
and higher temperature operation in a given size 
engine, but they result in concurrent power increases 
of a fairly large order. This may result in the power 
plant outgrowing an existing transport airplane. Does 
Mr. Russell consider that an existing airplane can be 
adapted to such a rate of change within present air- 
worthiness limitations at reasonable costs in develop- 
ment and tooling by increasing operating altitude 
and/or cruising speeds? Would he consider “‘shed- 
ding’’ or step-by-step elimination of power sections as 
a feasible alternate? It is reasonable to assume that 
a relatively small number of turboprop engine sizes 
will be developed, tooled for production and available 
to commercial transport operators at reasonable cost 
in view of the variety and kinds of propulsive units 
which the military services can support. The prob- 
lem of keeping transport aircraft and available engines 
matched may involve a lot of thought and seriously 
affect the ‘‘optimum”’ power plant availability too fre- 
quently assumed. It is with this consideration that it 
is my belief that new turboprop transports not now in 
existence are unlikely to be designed and that the ex- 
pectation should be that turboprop conversions of de- 
veloped and tooled current reciprocating powered 
transport and cargo aircraft will predominate. It 
would also seem advisable that such airplanes cover 
military transport and cargo requirements, as well as 
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commercial, in order to maintain reasonable costs for 
all concerned. 

One other point is desirably clarified with reference 
to power plant plus fuel weights of turboprop and 
reciprocating engines. Frequently, best economy 
fuel-air ratios of 0.058 to 0.062 are used on reciprocat- 
ing engines when compared to those of current allow- 
able operating temperatures on turboprop engines. It 
appears to me that most of the reciprocating engines 
are about as far away from actual satisfactory opera- 
tion at their optimum cruise condition because of the 
lead fouling problem on high-octane fuels as turboprop 
engines are from operation at 150°—200°F. higher tur- 
bine operating temperatures. If this is so, the initial 
weight advantage of the turboprop is available for 
any range the airplane may fly, and the cost of fuel 
used maintains approximately a one-third advantage 
in the turboprop engine. Other factors affecting this 
question are the oil tank weights and oil consumption 
weight and cost in favor of the turboprop and the cur- 
rent availability of these engines at well under !/2 Ib. 
per e.h.p. Would these considerations change the 
conclusions regarding relative range and pay-load 
capacity of the turboprop and reciprocating-engine 
installations? 


Dr. N. J. Hoff (East Coast), Professor of Aero- 
nautical Engineering, Polytechnic Institute of Brook- 
lyn: Of the 13 Wright Brothers Lecturers, eleven have 
been research scientists, professors, and research ad- 
ministrators, and one a test pilot. Mr. Russell is the 
first designer to be asked to deliver the Wright Brothers 
Lecture. This is a special reason for us to give hima 
cordial welcome, because the achievements of Wilbur 
and Orville Wright in the field of airplane design cer- 
tainly did not lag behind their accomplishments as re- 
search men and pilots. Asis well known, Mr. Russell 
is the designer of the Brabazon I, which flew for the first 
time last September. In the design of this 300,000- 
Ib., 20,000-hp. transatlantic transport, many difficul- 
ties had to be overcome. Those familiar with the work 
going on at the Bristol Aeroplane Company know that 
the problems were attacked with the best available 
methods of experimental and theoretical research and 
that no effort was spared to get the optimum result. 
Mr. Russell's paper gives a picture of the most up-to- 
date approach to the problem of the design of very 
large transport aircraft. It will undoubtedly be used 
as source material by those interested in today’s 
state of development of the art. 

Perusal of the manuscript shows the tremendous 
change airplane design has undergone in the last 10 
Before World War II, airplane stress analysis 
In the design 


years. 
was based almost entirely upon statics. 
of the Brabazon, static conditions were replaced to a 
large extent by the analysis of dynamic loads. Almost 
half of Mr. Russell’s paper deals with problems of 


aeroelasticity. To save weight, Mr. Russell found 
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himself compelled to install in the Brabazon a gust 
alleviator, whose purpose is to cut off the peak loads 
due to dynamic effects. The alleviator permitted g 
saving in structural weight amounting to approxi- 
mately 2 per cent of the gross weight of the plane, 
which meant an increase in pay load of about 5,000 
Ibs. From the economic standpoint, this increase may 
make the difference between failure and success jn 
commercial operations, and it is most encouraging 
that both an airplane and a transport company are 
willing to take the bold step of trying out such a de- 


vice. The step is undoubtedly a bold and imaginative 
one. The installation is complex mechanically and 


electronically, and many a transport company will be 
reluctant to use an airplane whose structural integrity 
depends upon the proper operation of such a complex 
installation until it has been proved over a number of 
years. The problem is one of added pay load versus 
the possibility of structural failure under adverse 
weather conditions due to improper operation of a 
complex servo installation. But progress cannot be 
made without courageous decisions, and a pilot from 
the time before World War I would consider the me- 
chanical and electronic equipment of our smallest and 
simplest personal planes as unnecessarily complicated 
devices compromising the safety of flying. The ques- 
tion arises now whether it would not be more suitable 
to try out the operation of a gust alleviator on a heavy 
bomber. The weight saved by the alleviator could be 
used to provide more protective armament and armor. 
In wartime, the accruing increase in safety would 
more than offset the danger of a mechanical or electric 
Of course in 


breakdown of the alleviator. 


time, safety considerations are similar for bomber and 


peace- 


transport. 

As the diameter of the Brabazon fuselage is about 
200 in. and the stringer spacing is small, one might ex: 
pect the fuselage to fail in general instability when 
subjected to a static bending test. An approximate 
analysis, by the method of \.A.C.A. Technical Note 
No. 1499, of data received from Mr. Russell yielded a 
value of 214 for the buckling index. On first sight 
one would conclude, therefore, that the frame spacing 
is too small and that weight was wasted by incorporat- 
ing too many frames in the fuselage. It turns out, 
however, that the general instability stress by far 
exceeds the proportional limit of the material and that 
the reduced modulus of the stringer decreases con- 
siderably the value of the buckling index. The calcu- 
lations could not be carried out completely because the 
Brabazon fuselage is outside the range covered by 
Technical Note No. 1499, and the material properties 
were not known exactly to this discusser. Never- 
theless, there is little doubt that the danger of general 
instability is only slightly greater than that of panel 
failure. 
initiated by the local instability of an unsupported 
edge of a stringer. 


It appears most likely that failure would be 
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Mr. Russell and the Bristol Aeroplane Company 
must be complimented on the amount and excellence 
of the work that went into designing the compression 
members of the wing structure. The arrangement 
having the maximum strength-to-weight ratio was 
determined analytically in England, while at the same 
time in America the N.A.C.A. carried out a long series 
of experiments to the same end. In this connection 
two remarks may be in place. First, it appears that 
scientific intelligence between our two countries was 
not very good, because the two groups did not know 
of each other’s work until both the analysis and the ex- 
periments were completed. Second, the difference in 
the approach to the problem seems to be typical for 
our countries, in spite of the rapid increase in theoreti- 
cal research in the United States in the last 20 years. 
Apparently, we have a little more money for expensive 
experimentation, while the British perhaps can get 
highly trained analysts at a somewhat lower cost. 

In conclusion, I should like to quote from a lecture 
delivered by Major J. D. Rennie in order to show how 
much aviation has developed since the early twenties. 
The paper was printed under the title ““Some Notes on 
the Design, Construction, and Operation of Flying 
Boats’ in the Journal of the Royal Aeronautical Soct- 
ety, Vol. 27, No. 148, p. 123, April, 1923: 

“Over long sea routes, say one hundred miles 
and upwards, the flying boat possesses over- 
whelming advantages over its only rival the 
steamship. .... Taking the maximum speed of 
a ship as 20 knots and the cruising speed of a fly- 
ing boat as SO knots, the voyage would be ac- 
complished in a quarter of the time.” 


Ray D. Kelly (East Coast), Superintendent, Tech- 
nical Development, United Air Lines: It ismostrefresh- 
ing to find so comprehensive an analysis of the charac- 
teristics and design problems of the turboprop transport 
aircraft as has been given us by Mr. Russell in this 
Thirteenth Wright Brothers Lecture. 

So many proponents of turbine-powered aircraft 
have jumped so quickly, and apparently so easily, to 
the turbojet configuration that some of us have been 
left gasping as to how the practical operating problems 
are to be mastered. Mr. Russell does not attempt to 
approach closely to these practical operational prob- 
lems, but at least he has established broad parameters 
for all to see and appreciate, so that the day-to-day 
experiences and requirements of the civil aircraft 
operators may be more effectively applied to turbo- 
prop-powered aircraft. 

In general, this lecture indicates that the turboprop- 
powered commercial transport has great potential 
possibilities for air-line operations, although a transi- 
tion to it from our present reciprocating-engine- 
powered craft represents a step of considerable magni- 
tude from an operations standpoint. The apparent 
lesser degree of flexibility will not be easily adapted to 
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and patterns. 


present operational methods The 
higher rate of fuel consumption, particularly at other 
than optimum altitudes, is significant. The probable 
need to design the aircraft to fit a specific route seg- 
ment and to build, generally, large aircraft (up to 
300,000 Ibs.) makes the venture into the new field 
somewhat alarming when viewed from the standpoint 
of privately financed air lines operating in stiff compe- 
tition. 

However, if the changeover to turboprops appears 
to be of considerable magnitude, that of going directly 
to turbojet has many more disturbing elements. To 
me at least, the turboprop appears to be a logical next 
step for the medium- to long-range civil air transports, 
and I am grateful to Mr. Russell for giving us so com- 
prehensive an analysis of the pertinent design criteria. 

It is believed that the following good operational 
characteristics can be claimed for a properly designed 
turboprop transport and that these are consistent with 
Mr. Russell's design criteria: 

(1) Improved take-off and climb ability, making it 
possible to operate still faster and heavier transports 
out of existing airports. (Enlarged landing areas 
close to population centers are not easily, cheaply, or 
quickly achieved. Many major airports are marginal 
for present equipment.) 

(2) A reasonable and usable increase in cruising 
speed of about 100 m.p.h. can probably be achieved 
without greatly deviating from the aerodynamic char- 
acteristics or safety features of present modern air- 
craft. 

(3) Ground handling characteristics will be essen- 
tially equivalent of present piston-engined craft. 
The important feature of propeller braking will be re 
tained. 

(4) The economic considerations for turboprops 
should not be too greatly different from those of 
present reciprocating engines. Even now, it appears 
that the fuel costs of the turboprops may not be too 
much greater on a ton-mile basis. Mr. Russell seems 
to think that this is true. If this is the case, then 
there seems to be every reason to believe that the later 
improvements in turbine efficiencies will possibly give 
specific fuel consumptions that are better than that for 
the current reciprocating engines. (This is extremely 
important to the air lines, because, even at the present 
time, some 10 to 12 cents out of every revenue dollar 
must be immediately paid out to the fuel producers 
for aviation gasoline.) 

It is encouraging to note Mr. Russell’s comments 
that passenger comfort with relatively large airplanes 
and at speeds in the range of 400 to 500 m.p.h. is 
likely to be improved. The question of passenger 
comfort at high speeds is one that is asked by all 
riders, inasmuch as most of them have had some un 
happy rough-air experiences in present equipment. 
Even some of the information which has been received 
from military sources has surmised that it may be un 
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wise to permit passengers to move about the cabin in 


very high-speed aircraft. The implication has been 
that they must of necessity keep the seat belt fastened. 
Certainly such a procedure would be a considerable 


deterrent to increased passenger traffic. 


Few people recognize the inability of high-speed air- 
craft to maintain block speeds that approximate the 
cruising speeds. Recently a study was completed by 
United Air Lines wherein comparisons between block 
speeds and cruising speeds were made for a number of 
contemporary aircraft, two of which were turbojet, 
one of which was a turboprop, and the other two were 
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relatively transports. Fig. | 
illustrates the results of this study. It will be seen 
that the two fastest airplanes (turbojets) quickly lose 
much of their inherent speed for the shorter nonstop 
In fact, they show little advantage for the 
first 300 miles. Airplane A is a turboprop, and it call 
be noticed that it quickly benefits from its increased 
speed over the reciprocating-engined aircraft. At 
1,000 miles its block-to-block speed is 287 m.p.h., as 
contrasted with its cruising speed of 320 miles. 


new  piston-engined 


distances. 


This study particularly emphasized the problem jet 
transports have in climbing to high altitudes and then 
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The turbo- 


coming down quickly at the destination. 
prop suffered somewhat but not to the same degree, and 
in this case it was possible to make comparisons with 
the turboprop craft operating at the same altitude as 
that normally utilized by reciprocating-engined craft. 

Fig. 2 represents another phase of this study wherein 
a pay load-range comparison was made for three of 
the subject airplanes. It will be noted that the pure 
jet airplane was not capable of going beyond about 
100 miles, when it required a fuel reserve equivalent to 
200 miles plus 45 min. plus a 10 per cent trip fuel 
safety allowance. Airplane A, using turboprops, also 
lost pay load rapidly for distances beyond 425 miles, 
because it also used up fuel at a fairly high rate. As 
would be expected, Airplane B, while carrying the 
smallest pay load for short distances was, because of 
the lower fuel consumptions of the reciprocating en- 
gine, less adversely affected by offloading of pay load 
on account of fuel, and therefore it ends up with the 
greatest amount of pay load for a nonstop distance of 
1,000 miles. 
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Fig. 3 is a comparison of three similar airplanes 
operating with normal fuel reserves on a 500-mile trip 
length. It will be seen that the airplanes start out 
with essentially the same pay load for this length trip. 
The gross weight of the piston-engined airplane for 
take-off is 43,200 Ibs. However, the take-off gross 
weight of the turboprop increases to 50,000 Ibs., pri- 
marily to accommodate the additional fuel required 
by this type of power plant for the trip. By the same 
token, the turbojet aircraft must take off with a gross 
weight of 60,000 Ibs. 

An inspection of Fig. 3 will now show that, if these 
same airplanes were to be operated on a 1,000-mile 
trip length instead of 500 miles, a good portion of the 
pay load for the turboprop would have to be offloaded 
in favor of fuel. In the case of the pure jet airplane, 
it probably could not be operated for a 1,000-mile trip 
length and still retain the necessary reserves, particu- 
larly if it must be assumed that the airplane will some- 
times descend from 30,000 ft., make an approach to its 
original destination, find that it cannot land and then 
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FIGURE 3 
WEIGHT BREAKDOWN | 
THREE 40 PASSENGER AIRPLANES WITH DIFFERENT ENGINE TYPES | 
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(INCLUDING FUEL REQUIRED FOR MISSED APPROACH.) 
A eee catind. PAYLOAD FUEL OPERATING WEIGHT EMPTY 
{10,400 lbs. [38004 29,000 1bs. 43,200 lbs. MAX. T.0. GROSS WT, 
ENGINES , 
BLOCK SPEED 225 MPH 
CRUISE " 260 MPH 
Bese + PAYLOAD FUEL OPERATING WEIGHT EMPTY 
jee ti 9,500 lbs} 76008 32,900 1bs. 50,000# MAX. T.0. GROSS vr, 
BLOCK SPEED 272 MPH 
CRUISE " 320 MPH 
a PAYLOAD FUEL OPERATING WEIGHT EMPTY 60,000 MAX. 1,0 
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must climb to the most economical altitude and pro- 
ceed to its alternate 200 miles away, and still have the 
legal 45-min. reserve. 

This analysis has made it apparent to us that, until 
the time when specific fuel consumptions for turbines 
can be improved greatly, the turbojet transport air- 
plane becomes inflexible and its success will depend 
upon its ability to dispense with reserves and alter- 
nates. 
critical in these respects and, for that reason, seems 


The turboprop appears to us to be much less 


to be a more logical next step. 

One other point, the variable high-velocity wind 
conditions that are associated with operation at high 
A study by United 
Air Lines’ meteorologists has resulted in two charts 


altitudes should be emphasized. 


analyzing the average wind velocities that exist at 
30,000 ft. in January and in June over certain seg- 
ments of United Air Lines’ route. The high-velocity 
westerly wind condition for January between Chicago 
and New York, shown on Fig. 4, is particularly per- 
tinent. 
be expected at this altitude and that they will be be- 


It will be noted that westerly head winds can 


tween 80 and 100 m.p.h. for 19 per cent of the time, 
between 100 and 120 m.p.h. for 20 per cent of the time, 
and between 120 and 140 m.p.h. for 18 per cent of the 
Still higher velocities exist for smaller per- 
centages to give a total of 71 per cent above SO m.p.h. 


time. 


Other route segments suffer from high wind velocities 
at 30,000 ft. in January but to a lesser extent. The 
relative frequency of high-velocity westerlies dimin- 
ishes during the spring to a minimum in summer, but 
the next chart, Fig. 5, for June, still shows some in 
tense westerlies. It is believed that these conditions 
essentially demand that any turbine-powered air 
transports will need to be operated at lower altitudesor 
make fuel stops under these conditions. Much more 
favorable wind conditions will often be found for west: 
bound trips during the winter months at low altitudes, 
as has been indicated in recent reports. (One of these 
by Namias and Klapp of the U.S. Weather Bureau, 
was published in the Journal of Meteorology, October, 
1949, and carries the title of ‘‘Confluence Theory of 
the High Tropospheric Jet Stream.’’) 

The air transport operators have found the newer 
postwar four-engined transports—specifically, — the 
DC-6, the Constellation, and the Stratocruiser—to 
have greatly improved schedule reliability, efficiency, 
passenger comfort, and independence of weather. 
The ability of these aircraft to fly long nonstop dis 
tances over, around, or under unfavorable weather 
It is hoped that these 
recent advantages, achieved as the result of improved 


has been of greatest importance. 


flexibility of operation, will not be overlooked by the 


turbine-engined aircraft manufacturers in _ thei 
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Sample distribution of wind components at 30,000 ft. during January, 1947. 
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Fic. 5. Sample distribution of wind components at 30,009 ft. during June, 1946. 
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LARGE TRANSPORTS 
emphasis on advertised cruising speeds achieved under 
optimum conditions. 

Mr. Russell’s lecture gives us assurance that he and 
his company are approaching the design problems of 
the turboprop aircraft from a realistic point of view. 
As a result, we have increasing confidence that the 
turboprop has an important place in the future of civil 
air transports and that it offers us a lesser degree of loss 
of recently acquired ‘‘flexibility’’ than the turbojet. 

I consider it a distinct honor to have been asked to 
discuss this comprehensive lecture given by A. E. 
Russell. We feel that both he and the Bristol Aero- 
plane Company, Ltd., are due a vote of thanks for the 
disclosures they have made concerning design criteria 
of the turboprop-powered aircraft. We hope that 
they will be encouraged to continue their work and 
that, because of this effort, the aircraft industry and 
the air transport operators will ultimately be able to 
offer to the public safer, more reliable, and more 
economical air transportation to every part of the 
world, having timesaving potentialities and passenger 
comfort features exceeding the standards of today, to 
which we point with considerable justifiable pride. 


Abe Silverstein (East Coast), Lewis Flight Pro- 
pulsion Laboratory, N.A.C.A.: The lecturer is to be con- 
gratulated for his excellent presentation of a timely and 
comprehensive treatment of the problems facing the de- 
signer of the large airplane and for his exposition of how 
these problems are resolved. The versatility and scope 
of knowledge displayed by Mr. Russell in preparing this 
most interesting and informative paper serve not only to 
reflect credit on the lecturer but to illustrate 
dramatically the many diverse and complex fields of 
study which the designer must master if he would seek 
toemulate Mr. Russell's successes. 

To the researcher in the aeronautical sciences the 
paper provides both comfort and despair—comfort in 
the recognition of the many problems that have been 
exposed for study and despair in the feeling that these 
problems accumulate in some multiple ratio of the rate 


also 


at which they are solved. 

The conclusions reached from the analyses and dis- 
cussions presented in the lecture are believed to be 
adequately tempered with engineering judgment, and 
the discusser found no large areas of disagreement. 

The author has chosen a most difficult task for the 
turbine engine and the one that it is not now best quali- 
fied to accomplish—namely, that of long range at high 
altitude. If he had considered it for medium range at 
medium altitude, his task would have been simplified. 
Mr. Russell states that airplanes powered by recipro- 
cating engines and by turbine engines could be de- 
signed for about equal range at high altitude, with the 
turbine-powered airplane about 100 m.p.h. faster. 
If he had considered the modern reciprocating engine, 
compounded with exhaust turbines feeding power into 
the crankshaft, it is believed that his analysis would 


WITH 
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have revealed a different answer, for it can be shown 
that, at altitudes of 35,000 and 40,000 ft., the range of 
an airplane equipped with a compound engine is 
greater than that of the airplane powered with turbine 
engines that operate under the modest values of pres- 
sure ratio, peak temperature and efficiencies assumed 
in Mr. Russell’s paper. Further, the cruising speeds 
of the airplanes equipped with compound and turbine 
engines would not be greatly different. 

Despite the author's failure to make a comparison 
of the present turbine engine with the compound en- 
gine, which would have revealed the inferior position 
of present-day turbine engines with respect to range at 
high altitude, the conclusion he has drawn that large 
prototype turbine-powered airplanes should be con- 
structed is believed to be a sound one. The turbine 
engine must be accepted now as a long-range engine 
on the basis of its future potentialities, which appear 
brighter than those of the compounded reciprocating 
engine. The gains to be made in the next few years 
from increasing the turbine cycle pressure, peak cycle 
temperature, and component efficiencies will place the 
turbine in a sound competitive position on the basis of 
performance at altitude with the attendant advantages 
of reduced vibration, larger power units, and simplifica- 
tion—swinging the pendulum in favor of the turbine 
engine. 

Since structural considerations, safety, and comfort 
in rough air plague the engineer as speeds increase, any 
reduction in the experience of rough air is worth-while. 
Perhaps the value of radar for this purpose has been 
underestimated. The use of radar has been shown to 
reduce gust experience through the avoidance of some 
turbulent clouds, and on the average, therefore, the use 
of radar would be expected to be beneficial. As in- 
formation is collected, it is not inconceivable that 
there may be weight allowances for the use of such de- 
vices where convective clouds are a major operational 
problem. While clear air turbulence is a factor, the 
more severe gusts appear to be associated with convec- 
tive cloud formations. 

The introduction of turbine-propeller power plants 
with higher operating altitudes may result in reduced 
gust load experience, with consequent benefits as ex- 
perience is gained. It would be expected that, by 
virtue of avoidance of clouds that cause difficulty at the 
lower levels, the intensity of gusts experienced at 
cruising speed might be appreciably lower for the higher 
The current ICAO requirements 
recognize this possibility. Future research and data 
collection for operations visualized by the author may 
result in radical changes in the statistics of gust load ex- 


level operations. 


perience. 

Gust alleviation still poses many problems of flying 
qualities, flutter, structural loads, and comfort, as so 
well recognized by Mr. Russell. With regard to ulti- 
mate strength, maneuver load factors impose a mini- 
mum strength that may not permit much gain to be 
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realized from gust alleviation. From the point of view 

of comfort, control, and fatigue life, it offers the greatest 
benefits, since it reduces all gust effects. If fatigue is 
a serious factor in the wing design, alleviation can be 
the difference between a good and a bad airplane. 

In passing, some problems of alleviators not men- 
tioned by Mr. Russell are the fact that a gust detector 
on the fuselage detects the local value of the gust and 
that, when the myriad gust shapes are considered, the 
possibility of a spurious signal due to a gust at the 
detector opposite to that acting on the wings is a matter 
of concern. The alleviation and the 
time lag between detection and action with speed are 
other questions that will need careful consideration. 
It is believed from Mr. Russell’s remarks that these 
items are under consideration, together with many 
others. 

The important effect of air damping on elastic re- 
sponse in gusts discussed by Mr. Russell has been 
recognized. It is felt, however, that generalization at 
this time is not warranted, since, as speeds are increased 
to sufficiently high values, the air damping may reach a 
maximum and decrease drastically at any higher speeds. 

Although not mentioned, the added conditions to be 
investigated by the designer are probably of serious 
concern as to engineering costs. The hope of simple 
solutions through analogs or models bears directly on 
this problem. The determination of the proper aero- 
dynamic and structural coefficients remains a serious 
problem. As aircraft get larger and carry more mass 
in the wing, the precision of determining dynamic air- 
plane response must become greater. 

It is of interest to note that in the discussion of flutter 
and vibration the lecture brings to bear on the problem 
several combined modes of attack, including flight 
testing, normal mode calculations, use of models, ana- 
logs, and high-speed digital computing machines. 

The discussion on gust detection and alleviation and 
the airplane structural problems is based on informa- 
tion obtained from Philip Donely, of the N.A.C.A. 
Langley Aeronautical Laboratory. 

Let me close by assuring Mr. Russell that all of us 
have learned a great deal from his work and its dis- 
cussion, and I offer him my best wishes for continued 
success in his pioneering work on large turbine-powered 
airplanes. 


variation of 


Frank W. Fink (West Coast), Chief Engineer, Con- 
solidated Vultee Aircraft Corporation: It isa pleasure 
to be asked to comment on Mr. Russell’s very fine 
paper. Just 2 months ago, I had the pleasure of visiting 
the Bristol Aeroplane Company in England where Mr. 
Russell showed me around. 

The subject of Mr. Russell’s paper is of extreme in- 
terest at the present time in the United States, and I 
should like to congratulate him on his splendid presen- 
tation in this field. Mr. Russell states that his paper 
embraces a smaller field than that of the scientists who 
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have made papers for earlier lectures in this series 
however, I believe he has covered rather a broad field, 
The paper touches on many points, any one of which 
could be the subject of an entire presentation—gych 
as turbine engine material, turbine blade Cooling, 
turbine using heat-exchanger cycle, optimum Wing 
loading, power loading, span loading for a turboprop 
airplane, etc. 

One point, however, that he did not bring up but 
which is of great interest in this country is the contro] 
of the turboprop engine. I mention it here merely 
because the control for the type of turboprop Mr. Rus. 
sell describes in this paper can be extremely simple. 
Since the turbine stage that drives the propeller is iy 
no way mechanically interconnected with the turbine 
stages that drive the compressor, there are no problems 
of high torsional shaft stresses caused by an unbalance 
between the prop load and the compressor load. 

Mr. Russell confines his paper to the long-range 
turboprop transport airplane. I believe something 
should be said for use of turboprop in a short-range 
transport airplane. 

It is true that the fuel consumption will be some- 
what higher than a reciprocating engine for some time 
to come, but it will be much lower than for a jet engine 
of the size that will give equivalent take-off thrust 
and when operated at medium aititudes. 

It is pointed out in the paper that for endurance 
below 5 hours, a comparative low compression ratio 
engine may be used. A lower compression ratio 
turboprop appears also to fit in better with the re- 
quirements of a short-range transport, in that it will 
operate efficiently at a somewhat lower altitude, say 
20,000 to 25,000 ft., and is less susceptible to high am- 
bient temperatures. A short-haul airplane is likely to 
have to operate continuously during some seasons of 
the year under tropical or hot-day take-off conditions. 

It is believed that a turboprop engine can be oper- 
ated just as economically on short-haul airplanes asa 
reciprocating engine, since the engines have much 
lower maintenance costs and should eventually be just 
as economical to build and, further, since the ex- 
tremely low vibration should greatly decrease air- 
plane maintenance costs—to say nothing of increased 
passenger travel generated by the increased comfort. 
A turboprop also has a distinct advantage over turbo- 
jet short-haul transports, in that short-haul trans- 
ports must, by necessity, operate from many airports 
that will not see improvements necessary for jet trans- 
ports for many, many years to come. Propeller- 
driven airplanes will be able to take off on shorter 
fields, land on shorter fields because of reversing, and 
cope better with traffic conditions at fields that are 
not completely equipped for blind landings. 

One interesting point that Mr. Russell made was 
that of the present range of propeller-driven sizes now 
being manufactured or developed; the specific weight 
of similar designs appears to be roughly constant and, 
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LARGE TRANSPORTS 
therefore, in certain cases it is beneficial to couple two 
small engines together in order to make a more com- 
pact and more suitable installation for burying in a 
wing. This happy condition could also be used to 
geat advantage in short-haul transport airplanes, 
since an airplane could be designed with four engines 
paired together so that only two nacelles and two pro- 
pellers would be required. This would give the com- 
pactness and aerodynamic efficiencies of a two-engined 
airplane with the reliability and one-engine-out per- 
jormance of a four-engined airplane. 


E. G. Stout (West Coast), Assistant to the Chief En- 
gineer, Consolidated Vultee Aircraft Corporation: Mr. 
Russell pointed out that there isnoadvantage in cutting 
out turbine units on multiunit airplanes. His conclu- 
sions are based, however, on the ideal condition of 
cruising at optimum altitude, where the full power of 
all units can be used efficiently from a fuel economy 
standpoint. For the condition where a specified alti- 
tude must be held, as in the case of radar search with 
a military airplane, sizable improvements in range 
are obtainable by appropriate selection of units. For 
axample, the Convair XP5Y-1 T-40 propeller-turbine- 
powered long-range flying boat, which for certain 
radar search missions is operated continuously at 


10,000 ft., shows increases in range as follows: 

All unit operation Basis 

75 per cent units operation 30 per cent increase 

50 per cent units operation 65 per cent increase 

105 per cent increase (or over 
double the range 
able without cutting units) 


50 to 25 per cent units operation 
obtain- 


Of course, the economical cruising speed drops as 
the range is increased (following the usual speed-range 
relation). 

E. C. Wells (West Coast), Vice-President 
neering, Seattle Division, Boeing Airplane Company: 
be commended for his 


-Engi- 


Mr. Russell should indeed 
courage in selecting the broad field of aeroelastic 
phenomena and dynamic loads for a portion of his 
paper. Considering the scope of his selected subject, 
he has come forth with a thorough and straightfor- 
ward discussion of the problems and the general sig- 
nificance of design trends. 

As he mentioned, the dynamic and aeroelastic as- 
pects of the structural loads problems have only re- 
cently become of importance in the design of basic 
structure. Thus, as in any new and rapidly expand- 
ing phase of design, we find ourselves proceeding into 
somewhat unfamiliar territory with little more than a 
teasonable understanding of the physical relation- 
ships involved together with analytical and experi- 
mental tools through which we can increase this under- 
standing. 

Today we find ourselves introducing the additional 
parameter, time, into our structural loads concepts. 
Unfortunately, a large percentage of the experimental 
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data that formed the foundation for the pseudostatic 
structural design criteria was taken without adequate 
recognition of the future importance of time as a vari- 
able in the loads problem. However, this foundation 
material should neither be arbitrarily discarded nor 
be blindly accepted as applicable. Rather, it should be 
re-evaluated, where possible, in the light of the new 
problem. Since time is now of significance, we are 
interested not only in the probability of occurrence of 
a load but also in the probability that it will occur 
with a given time history. As Mr. Russell has shown, 
the time history of the load application may be of 
equal or greater importance than the close definition of 
the maximum magnitude. Study of available gust 
data will reveal investigations of several airplanes 
from which graphs similar to Fig. | may be obtained. 
It will be noted on this plot of true gust velocity vs. 
gradient distance, that high velocity gusts are not con- 
fined to a narrow band of gradient length. X gusts 
of approximately Y velocity occur. There is there- 
fore a probability of X divided by the total number 
that this velocity will be attained. However, the 
probability of any given gradient length at this veloc- 
ity is considerably reduced. This difference of proba- 
bility should logically be reflected in the design. 

Another important feature of the gust problem 
which warrants some consideration is the velocity 
variation along the span and fore and aft during the 
passage of the airplane. The generally accepted 
linear chordwise gradient with the maximum velocity 
sustained and a constant velocity spanwise at any in- 
stant is open to question. The importance of the 
chordwise variation of velocity and therefore the time 
history of the externally applied loads can be seen on 
Fig. 2. Here is shown the dynamic response factor 
for the linear gradient with sustained maximum and 
for comparison the response to a force that increases 
linearly to a given value and thereafter decreases to 
zero at the same rate. The curves represent the re- 
sponse of a simple single degree of freedom lumped 
system without damping. The differences are pro- 
nounced for the shorter gradients, particularly those 
for which the sustained type loading indicates the 
greatest dynamic effect. Considerable evidence exists 
to support the use of forcing functions of approxi- 
mately the triangular time history rather than the 
sustained type for the usual gust. Again, probabili- 
ties become important. 

It is interesting to carry this triangular velocity dis- 
tribution concept one step further to examine its rela- 
tionship to the spanwise load distribution. As Mr. 
Russell has pointed out, statistically an airplane will 
enter a gust asymmetrically much more often than 
not. If the velocity variations in the atmosphere are 
indeed somewhat triangular, this would certainly be 
expected. The few measurements that have been 
attempted indicate a triangular or sinusoidal type 
variation to be equally as applicable to the spanwise 
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yelocity variation as to that in the direction of flight. 
It would, in fact, be somewhat naive to assume a dif- 
ference between the distributions in the two perpen- 
dicular directions, since it is unlikely that gusts would 
have any particular geometric alignment with the 
direction of flight. The distribution will have impor- 
tant influences on both the pseudostatic loading and 
the response of the elastic airplane. 

As with any problem, the best solution is the 
elimination of the problem. Certainly this is true of 
gust loads. Thus, the development of gust alleviators 
is a logical step. The development of aircraft for 
flight in the relatively high subsonic speed range has 
resulted in the use of sweptback lifting surfaces to 
obtain their favorable effects on compressibility phe- 
nomena. Fortunately, sweepback also produces a 
significant reduction in gust loads. 

The major alleviating effects of swept surfaces are 
twofold. As is well known, the slope of the lift curve 
is lower than that of the equivalent straight surface. 
Hence, the effective angle of attack change due to a 
given gust will produce less load change or accelera- 
tion on a Sweptwing airplane. Second, bending of the 
wing produces alleviating effects similar to those of a 
skew hinge. The trend can be conveniently illustrated 
by considering conditions of stabilized flight. Fig. 3 
shows typical lift curves for a wind-tunnel model with 
the wings deflected to the several positions correspond- 
ing to flight at different load factors and rigidly re- 
strained in these positions as the angle of attack is 
varied. The curves therefore individually represent 
lift vs. angle of attack for a rigid wing with a particular 
warpage of the plane of zero lift lines, and their slopes 
would include the loss due to sweepback mentioned 
previously. On each curve, a circle has been drawn 
at the approximate lift coefficient appropriate to the 
load factor, assuming a given constant air speed. The 
lift change would tend to progress along the dashed 
curve between these points in a gust. Hence, the 
elastic deflections produce an additional effective re- 
duction in the slope of the lift curve. 

Examining these elastic effects further, it is seen 
that the major changes in lift occur near the tips of the 
wing for the usual wing deflection shape—a favorable 
arrangement structurally. Not only do the elastic 
deflections of a sweptwing reduce the gust loading, but 
they also serve to shift the spanwise center of pressure 
of the remaining load inward. 

In some respects, sweepback is unfavorable as a 
gust alleviator. For example, for typical gusts en- 
countered symmetrically, the forward areas of the wing 
in the vicinity of the root enter the gust first with a 
resulting tendency to pitch the airplane in a direction 
to increase the effectiveness of the gust. Similarly, 
the inward shift of the aerodynamic load due to elastic 
deflection is accompanied by a forward shift. This 
latter phenomenon is closely related to the short- 
period stability of the airplane, however, and therefore 
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is somewhat compensated for in the selection of per- 
missible center-of-gravity positions. Thus, it is 
found that the usual swept surface is also conveniently 
a powerful gust alleviator without the ever-present 
disadvantages of a mechanical system. 

It is felt that Mr. Russell’s rather unfavorable views 
on the value of testing scaled flutter models are not 
entirely justified. At any rate, our own experience in 
that field has led us to believe that an extensive model 
testing program is essential as a supplement to ana- 
lytical flutter studies. Of course, it must be admitted 
that full-scale Reynolds Numbers cannot be ap- 
proached with models of reasonable size. However, 
the limited experimental data available in this coun- 
try tend to suggest that Reynolds Number is not of 
primary importance in flutter phenomena, except pos- 
sibly at transonic speeds where separation and com- 
pressibility effects may be fundamentally interrelated. 
Certainly, the existing theory of oscillating airfoils is 
not sufficiently general to encompass scale effects, and 
it may be argued that model data are no more suspect 
than analytical results in that respect. Also the criti- 
cal flutter modes of a large, flexible airplane generally 
involve simple distortions of the primary components; 
for that reason it is felt that a simplified model struc- 
ture may be employed without consideration of the 
effects of shear deflection. Actually the number of 
normal modes of the complete airplane which involve 
only simple deformations of the primary components 
is relatively large; hence, it is possible with a simpli- 
fied model to obtain a more comprehensive representa- 
tion of the dynamic characteristics of the airplane 
structure than would be possible at the present time 
in an analytical investigation of reasonable scope. 
Our usual practice is to conduct both analytical and 
experimental investigations concurrently during the 
development of a new model. Fig. 4 illustrates some 
of the essential features of one of our flutter models. 
It will be noted that the primary wing structure con- 
sists of a single metal spar of solid section, suitably 
proportioned to give appropriate distributions of 
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A similar arrange- 
The contour is 


flexural and torsional stiffnesses. 
merit can be adopted for the fuselage. 
accomplished by the use of cellular cellulose acetate of 
low density or balsa wood so arranged as to provide 
no additional stiffness over the basic spar. 

Unquestionably, one of the greatest advantages to 
be derived from a flutter-model test program is the 
ability to determine relatively quickly the consequence 
of design modifications and changes in flight condi- 
tions or airplane loading. Because of the simplified 
construction, it is not difficult to manufacture new 
components incorporating changes of stiffness or of 
geometric form. Springs of different stiffnesses are 
usually provided in order to investigate the effect of 
variations in control surface natural frequencies, and 
variations in airplane altitude are simulated by adding 
or removing weight from the model. With present 
computational facilities we have found that the paral- 
lel analytical program of flutter studies simply cannot 
keep pace with the experimental program. In spite 
of its limitations the flutter model does provide an- 
swers in time to be of use during the design of an air- 
plane, thus materially reducing the likelihood of costly 
changes during the later stages of development. 

Also, it has been our experience that flutter model 
tests generally reveal the existence of critical modes 
that might be overlooked otherwise because of the 
understandable tendency to oversimplify the analyti- 
cal representation of the airplane. These tests have 
served to reinforce substantially our conviction that a 
large, flexible airplane structure should be regarded as 
a single, interconnected dynamical system. In one 
particular instance, a moderate change in elastic re- 
straint of the elevator served to introduce a very large 
v.ing vibration component into a flutter mode, in 
addition to the tail and fuselage motions that were 
already present. As a matter of fact, four distinct 
flutter modes were found within the practical range of 
elevator frequencies in this case. The experimentally 
determined flutter speed is shown as a function of ele- 
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vator frequency in Fig. 5. Experiences of this king 
serve to lessen further our confidence in the validity 
of any sort of simplified flutter criteria. 

In the design of our recent models, we have als 
maintained dimensional similarity of elastic deforma. 
tions due to weight and steady air loads in order ty 
investigate other aeroelastic phenomena as wel] as 
flutter. Recently, consideration has been given ty 
the possibility of testing aeroelastic models at hig) 
subsonic Mach Numbers in order to determine com. 
pressibility effects. Unfortunately, it appears to be 
impractical to maintain dimensional similarity of 
flutter and the other aeroelastic characteristics, as wel] 
as compressibility effects, with a model of reasonable 
size. 

I certainly do not mean to imply that the analytical 
investigations should be abandoned or reduced jn 
scope. On the contrary, we are in agreement with 
Mr. Russell with regard to the desirability of exploit- 
ing the modern developments in electronic computers 
to the fullest extent for the analysis of aeroelastic 
phenomena. It is anticipated that these machines 
will ultimately permit us to conduct a sufficiently 
complete analysis of the entire airplane as a single, 
interconnected, dynamical system. Our own studies 
in this field tend to favor the digital computer. In 
order to employ an analog computer for flutter 
analysis, it is necessary to accomplish a considerable 
amount of preliminary computational work before the 
problem can be put on the machine. Also, the prob- 
lem of introducing such refinements as compressibility 
and aspect ratio corrections into the analog com- 
puter is difficult and possibly impractical. These 
modifications, as well as any others that may be con- 
sidered necessary in the light of future research, could 
be handled without particular difficulty on a digital 
machine. 

Mr. Russell’s comments on ground vibration testing 
of large aircraft have suggested an ingenious solution 
to a difficult and important problem—namely, the 
provision of a supporting system that will permit 
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Fic. 5. Flutter speed vs. elevator frequency. 
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Application of Simultaneous Equation 
Machines to Aircraft Structure 


and Flutter Problems 


P. A. DENNIS* ann D. G. DILLT 
Douglas Aircraft Company, Inc. 


SUMMARY 


Flutter, divergence, and control surface reversal are generally 
Not only aerodynamic but also 
Considerable mathe- 


known as the aeroelastic problem. 
elastic and inertia forces are involved. 
matical effort is required along with the development of a proper 
physical conception to deal with this problem. It is used as an 
example in this paper, since the time expended on this problem in 
modern aircraft design led the authors to the study of the linear 
algebraic simultaneous equation machine. It also represents a 
typical aircraft engineering problem that may be solved by such a 


machine. 


TYPES OF PROBLEMS 


—— problems may be classified as deter- 
minate or indeterminate. The simultaneous equa- 
tion machine bears the same relationship to algebraic 
indeterminate problems as the slide rule to determinate 
ones. Indeterminate problems may be further sub- 
divided into simple indeterminate, in which only one 
condition of equilibrium exists, and stability problems, 
where a multiplicity of solutions exists. 


SIMPLE INDETERMINATE PROBLEMS 


These may be expressed by ordinary simultaneous 
equations of the form: 


AyX; + ApXo + 


Ay3X3 occ = B, 
AnX,; + AwX2 + Ao3X3 


= B 


Some of the numerous aircraft engineering problems 
that may be so expressed are: (1) statically indeter- 
minate structural problems; (2) forced vibration prob- 
(3) span load distribution; and (4) control sur- 
face effectiveness. 

Some of these problems illustrate an additional as- 
pect of equation machines, in that integrations may be 
performed by breaking them down into increments. A 
single integration appears explicitly in the control 
surface effectiveness problem, and they are implicitly 
contained in some of the others. 

The equation computing machine solves these prob- 
lems by generating, multiplying, and adding electrical 
voltages in exactly the same form as the equations. 
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lems; 


STABILITY PROBLEMS 


The second class of problems is much more difficult 
than the first because of the multiple conditions of 


equilibrium possible, and they may be expressed 
mathematically in the form: 
(Ay N)X1 + AywX2 oe Aj32 3 coos = 0 


Ajo a+ (Ave — r)Xe + Ao3. fen os @ @ 


Here the additional variable \ introduces n possible 
solutions to the equations, representing m possible 
equilibrium conditions, where m equations are involved. 

Aircraft engineering problems in this class are, for 
example: (1) buckling; (2) free vibration problems; 
(3) divergence; (4) flutter; and (5) aero stability, 
static and dynamic. 

In order to handle these problems, an additional set 
of multipliers is supplied which are ganged together and 
vary \ as required by the equation. 

This \ term is also useful in some of the first class of 
indeterminate problems, although it is not essential to 
their solution as in the second class. 


MISCELLANEOUS PROBLEMS 


Other problems of a more general nature may also be 
solved, such as: 

(1) Algebraic polynomial equations (a machine with 
n equations will solve an uth order equation. The A 
term is required.) 

(2) Multiplication and inversion of matrices. (These 
require the solution of m ordinary simultaneous equa- 
tions, in general. A complete matrix multiplying and 
inverting machine, that could perform this in one in- 
stead of m operations could be built on the same 
principles. ) 

(3) Step-by-step solution of variable coefficient and 
certain nonlinear problems. 


COMPLEX NUMBERS 


Many of the above problems require complex numbers 
for their representation. The machine used as a basis 
for this paper is based on real numbers and handles 
complex equations by expanding them to twice the 


number of real equations. Thus, 
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(a + ib)(x + ty)*c + id 


is equivalent to 
ax — by =c 


bx + ay =d 


THE AEROELASTIC PROBLEM 


The aeroelastic problem consists of the combined 
aerodynamic, elastic, and inertia properties of an aero- 
dynamic surface which manifest themselves in the 
phenomena of torsional divergence, control surface 
effectiveness and reversal, and flutter. A corollary to 
these are the vibration modes and frequencies of the 
structure. The difficulties in dealing with these prob- 
lems may be seen from the fact that we must not only 
represent mathematically a complicated set of struc- 
tural and aerodynamic forces but also their performance 
as a time history because of the oscillatory nature of 
vibration phenomena. 

An aircraft structure has continuously distributed 
properties. For analysis purposes these are sectional- 
ized and represented as lumped parameters strung to- 
gether to approximate the true structure. The aerody- 
namic forces are dealt with in a similar fashion. Varia- 
tion with time is dealt with by assuming uniform sinu- 
soidal motion, and the forces may then be represented 
by rotating vectors with given magnitudes and phase 
relationships. If it is further assumed that all of the 
forces involved are linear with their related displace- 
ments, the problem is then formulated so that the solu- 
tion may be obtained. This composition of the problem 
is not ultimate, and improvements in it are being made; 
however, it is the present working theoretical basis. 

The above assumptions have made the problem 
solvable and even simple in the physical sense; how- 
ever, they have not left it easy to obtain answers. The 
expression of the above characteristics as force equa- 
tions leads to simultaneous equations of orders varying 
with the complexity of the problem from two complex 
equations in simple problems up to eight or more in 
cases requiring the additional refinement. The solution 
of these simultaneous equations, while in principle very 
simple, are, in fact, laborious, and some of the best 
mathematicians and engineers have applied their efforts 
to finding short cuts to reduce the labor involved. One 
approach is the use of an analogy computing machine, 
set up to represent the equation coefficients. 

It was believed that such machine solutions repre- 
sented a powerful and worth-while approach to the 
problem of saving labor and giving a physical feeling for 
the relative importance of various parameters in specific 
cases—that is to say, their variations are easily ex- 
plored. To this end a simultaneous equation machine 
was constructed. A four-element machine was chosen 
because it was the smallest machine that would show up 
the essential characteristics of the problem. The prob- 
lems solved on it are thus of a simple nature, but they 
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fully illustrate the principles involved in the solution of 
the more refined cases and thus form an interesting 
evaluation of the capabilities of this method. 


THE SIMULTANEOUS EQUATION MACHINE 


The equation machine solves the problem by multi. 
plying and adding voltages in exactly the same form as 
the mathematical expressions. A potentiometer circuit 
is used for multiplication, R2> > R; and a parallel circuit 
for addition, R >> Z,. 

The machine thus multiplies directly and divides by 
the inverse process by adjusting the multiplier (un. 
knowns) until the solution is obtained. 


Basic CLASSIFICATIONS 


There are four basic types of simultaneous equation 
analogy machines. This classification does not depend 
on whether or not the machine is electrical, mechanical, 
hydraulic, or based on probability multiplication. Any 
physical process with a multiplicative law that is also 
capable of addition may be used as a basis for such a 
machine. The classifications are: (1) reversible, 
manual; (2) reversible, automatic; (3) semireversible, 
manual; (4) semireversible, automatic. 

Whether a machine is manual or automatic is a 
matter of design refinement rather than the mathe- 
matical principles on which it is based. The reversi- 
bility or the degree of reversibility is, however, a part of 
the mathematical design of the machine. 

The title ‘‘reversible’’ is used because machines based 
on reversible physical processes—such as the lever, 
gear, or transformer—immediately indicate the com- 
plete solution. A synthetic mathematical reversibility 
may be supplied to irreversible processes to achieve the 
same results. This may be either automatic or manual. 
If automatic, the answer is immediately indicated; if 
manual, as is this machine, convergence to it is guaran- 
teed. In this particular calculator, this is accomplished 
by minimizing the sum of the absolute values of the 
errors of each equation. This can only be zero when all 
equations are satisfied. The absolute value is taken by 
a rectifier before adding the error signals. 

Semireversible machines are based in an approxima- 
tion to the reversible process. The most common of 
these approximations is mathematically equivalent to 
the Gauss Sidel iteration procedure. These machines 
will converge for some problems and diverge for others, 
whether manual or automatic. These machines may be 
given the ‘‘synthetic’’ reversibility by considerable 
mathematical preparation of the problem. This also 
may be done in the machine at the expense of halving 


its capacity.! 


METHOD OF SOLUTION 


A problem is prepared for the machine by placing it 
in simultaneous equation form. The equations are then 





vol 





Fic. 


Fic. < 





ution of 
resting 


multi. 
form as 
circuit 
circuit 


ides by 
‘r (un- 


uation 
lepend 
anical, 

Any 
is also 
such a 
rsible, 
rsible, 


18 a 
1athe- 
versi- 
art of 


based 
lever, 
com- 
bility 
e the 
nual. 
i; @ 
iranl- 
shed 
the 
n all 
n by 


ma- 
n of 
t to 
ines 
ers, 
y be 
ible 
also 
ing 


x it 


SIMULTANEOUS EQUATION MACHINES 109 



























































































































































oa oe 1.5 | T T 
OCALCULATED WITH EQUATION MACHINE, @ *1,125 CPM 
x AEXACT CALCULATION, W =1I,II4¢ CPM 
R w 
UNIT ' 2 1.0 | 
A F 
a 
VOLTAGE Re a AT 
a 
AX ws ae 
- rs 
< 
O 2 
« 
R, »» R, 0 ; ; : 
Fic. 1. Multiplying circuit. € , 
WING SECTION 
Fic. 5. Comparison of torsion modes. 
R R | 4 LIFT (25 % CHORD) 
£1 + 5s ELASTIC 
. —?1 e — AXIS (37.5 % CHORD) 
E, E 2 a 6 
R > Ze Fic. 6. Wing section. 
Fic. 2. Adding circuit. 
¥ 240" an 
- " 
Fic. 7. Wing plan form. 
Fic. 3. Experimental simultaneous equation calculating 
machine. 
O CALCULATED WITH EQUATION MACHINE, Vo #664 MPH 
A EXACT CALCULATION, Vp * 858 MPH 
S 
raRRE anaemia 2 1.00 —— 
T 3% 
— | — Sz a 
_ a A 
| TORSIONAL! | — _| ws 
\ RIGIDITY INERTIA = 50 oe 
pt | ook 5 
| | | 3 it 
PET PTET 
| | : 
cone SPAN 1 2 3 4 
Fic. 4. Idealized four-section wing with concentrated torsion WING scones : 
springs. Fic. 8. Comparison of torsional divergence 











110 JOURNAL OF THE 


multiplied by such factors as necessary to bring the 
largest coefficients near 1, since the machine is cali- 
brated to handle numbers between +1 and —1. The 
coefficient dials are then set, and the unknowns are 
manipulated until the output signal is zero. 


In frequency problems, the ratio of the unknowns 
(amplitudes) is the desired information, while the over- 
all value of the unknowns (amplitudes) is arbitrary, so 
one of them is fixed and the diagonal term is manipu- 
lated in its place. In the flutter problem, the phase and 
amplitude of one component are fixed, and the diagonal 
and one of the constant terms are used in their place. 
The sign and magnitude of this constant term are indica- 
tions of the stability of the system at the conditions of 
the computed point. There are many possible varia- 
tions to the above technique. 


In general, because of the full reversibility of this 
machine, any 7 of the dials—chosen from the unknowns, 
coefficients, constants, or \ term 
as unknowns, depending on the information desired. 
Such freedom of procedure has been little discussed 


mathematically and is believed to have great possi- 


may be manipulated 


bilities. 


ILLUSTRATIVE PROBLEMS 


These problems, solved by the four-equation machine, 
are illustrative in character. For the sake of concrete- 
ness, an actual wing is invisaged, and typical properties 
are assigned to it, as shown in Figs. 4, 6, 7, and 9. 


For many practical cases, it would be desirable, of 
course, to divide the wing into more than four sections 
and to otherwise refine the treatment used in these 
sample problems. It is believed, however, that the 
basic principles of the use of the machine are adequately 
covered and that further detail refinements in tech- 
nique will be sufficiently ohvious of those familiar with 
these problems. 


The following problems are solved, and the solutions 
are compared to those obtained by a hand calculating 
machine: (1) first natural torsion vibration modes and 
frequency; (2) torsional divergence; (3) aileron effec- 
tiveness; and (4) wing bending—torsion flutter. 


WING PROPERTIES USED IN 
ILLUSTRATIVE TORSION PROBLEMS 
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IDEALIZED FOUR-SECTION WING WITH CONCENTRyp 
TORSION SPRINGS 


Spring Constants Mass Inertias 


; 1X 10° 
C=-< 23x ihe ; 
28.41 
zm, Oe 
G-—0wxe h= 
28.41 
ae ' 1 xX 10! 
C3 =) xX 106 Ts = 
28.41 
0.5 X 108 
C=25xX%10° = — : 
28.41 


DEFLECTIONS 

The twists along the wing for general applied torgy 

may be written: 
(1/Ci)(T1 + T2 + T3 + T,) 

a = a + (1/C2)(T1 + To + Ts) 

a = as + (1 C3) (T3 + Ts) 

a4 = a3 + (] C,)(T4) 
or 

a = 0.057; + 0.0572 + 0.0573 + 0.057; | 
0.057, + 0.157% 4. 0.1573 + 0.157; 
aga = 0.057; + 0.1575 + 0.3573 aa 0.357, | 
= 0.057; + 0.157, + 0.3573 + 0.757; 


Qa = 


II 


X 10 


(a) The Frequency Problem 


In this case, the 


forces: 
T; = Wha 
T. = Wleae 
sae Bae), 06 
T3 ai W *T30%3 x ' 
T; = W?2Is,a4 
Substituting, we have 
1/W2a, = 0.20aq, + 0.10a2 + 0.05a3 + 0.025q 
1/W2a. = 0.20q, + 0.3002 + 0.15a3 + 0.075a 
1/W2az; = 0.20a, + 0.3002. + 0.35a3 + 0.175q 
1/W2a, = 0.20a, + 0.30a2 + 0.35a3; + 0.375a 


or calling 1/W? = X, 


(0.830 — X)a; + 0.10a2 + 0.05a3 + 0.025; = 0 
0.20a, + (0.30 — X)az + 0.15a3 + 0.075a, = 0 
0.20a, + 0.30a2 + (0.35 — X)az3 + 0.175a; = 
0.20a, + 0.30a2 + 0.35a3 + (0.375 — X)ay = 0 


| 


which equations are solved for a’s and X’s on tht 


machine 
W = 1/¥V x in thousands of cycles per min. 
For the results of this calculation see Fig. 5. 
(b) The Torsional Divergence Problem 
In this case 


T, = gsiéidia,, T2 = 
7 = 


GS2€202 A2 


Ts = 530303a3, GS4040404 


torques are given by the inerty 


(qd 


where 
thous 
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SIMULTANEOUS 


= dynamic pressure 

= area of section 

= distance from elastic axis to 
q = slope of the lift curve 

(a, = a2 = G3 = A, = 6, Say.) 


'/, chord point 


where g in Ibs. per sq.in. = 17.8V? and where V is in 


thousands of m.p.h.). 


S, = 4,800 sq.in. é: = 10 in. 
S: = 3,830 sq.in. é = S8in. 
S; = 2,880 sq.in. és = Gin. = fs 
Ss = 1,920 sq.in. é = 4in. = fi 


symbol f refers to Fig. 9), or 


T, = 5.12 X 10° Via, 
T; = 3.26 X 10% Via 
Ts; = 1.84 X 10® V2a; 
T, = 0.82 X 10° V2a, 


Substituting this and setting XY = 1/V?(V in thousands 


of m.p.h.), 

(0.256 — X)a; + 0.163a2 + 0.092a3 + 0.04la, = O 
0.256a, + (0.489 — X)ae + 0.276a3 + 0.123a, = 0 
0.2560, + 0.489a. + (0.644 — X)a3 + 0.287a,; = 0 
0.256a, + 0.489a2 + 0.644a3 + (0.615 — X)ay = 0 


The results of the solution are shown on Fig. 8. 


(c) The Aileron Effectiveness Problem 


T’s under 
and 


In this case we have, in addition to the 
problem (b), torques due to the aileron in section 3 
4, 

Now the additional lifts are: 


L3 Gsx3kB, k == 
Ly = GSsQ4kB 
‘eo 3 = gs33kf36, 


7 ‘4 = gssaskfs B 








4 LIFT OUE 


E.A. TO AILERON 


C tt ~ p 


TAKEN FROM E.A. (37.5%) TO 50 % CHORD 




















Wing with aileron. 
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Fic. 11. Comparison of flutter speeds. 


Thus, 7’s are: 


T, = 5.12 X 10® V?ay 
T2 = 3.26 X 10° V2ae 
T; = 1.84 X 10® V?a3 + 0.92 X 10° V8 
T; = 0.82 X 10° Vay + 0.41 X 10° V*8 


and the lifts: 


L, = 5.12 & 10° V?ay 

L,; = 4.17 X 10° V*ae 

L; = 3.07 X 10° V?as e 1.53 & 105 V? 
Li = 2.05 X 10° Vay + 1.02 X 10° _ 


and the rolling moment: 


M = 30 Li a 90 Le aa 150 Ls + 210 Ly 
Neglecting wing twist, we would have a moment 


Mo = 150L3 + 210L;, 


Thus, 

a 1+ 0.347 + 0.835 @ + 1.04 = + 0.97 & 
Mo B 8 8 B 
and 

AM 


~ ia per cent increase of rolling moment due to wing 
- 0 
twist (actually a decrease) = 0.347 (a ‘B)+ 


0.835 (a2/B) + 1.04 (a3/8) + 0.97 (a4/8) 
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Substituting for the T’s, we obtain equations similar to 
that for torsional divergence: 


[0.256 — (1/V2)]ar + 0.16302 + 0.09203 + 0.041a, 


= 0.067 6B 
0.256a, + [0.489 — (1/V?)] ae + 0.276a3 + 0.12304 

= 0.200 B 
0.256a; + 0.489a2 + [0.644 — (1/V?)a3 + 0.287 a4 

= 0.465 B 


0.2560; + 0.489a2 + 0.644a3 + [0.615 — (1/V?)] a 

= 0.630 8 
(V in 1,000’s of m.p.h.). When these equations are 
solved for the a’s, the decrease in rolling moment is 
given by: 


in, ~ 0387 ($ )osss( Z)1.04($)o.7() 
= = 0.347 0.835 1.04 0.97 
My B 


To find the speed of aileron reversal, we put AM//Mo = 
—l and have B = —(0.347aq, + 0.835a2 + 1.04a3 + 
0.97 a4) 
and substitute in the above, obtaining an equation 
similar to that for torsional divergence. 

Results of the computations are given on Fig. 10. 


(d) The Flutter Problem 


The flutter problem is of a more complicated charac- 
ter than the problems treated previously, since ampli- 
tudes and phases of both bending and torsion of the 
wing are involved. 

We may, however, for purposes of illustration, em- 
ploy a ‘‘typical’’ section of the wing. In this case, we 
have only two degrees of freedom, bending and torsion; 
but, because the motions are not in phase, each degree 
of freedom is represented by a complex number instead 
of a real number. 

This situation can be handled on a four-equation 
machine by using one equation for the in-phase (real, 
say) forces and one for the out-of-phase (imaginary) 
forces in each degree of freedom. In substance, a two 
degree of freedom complex system is replaced by a four 
degree of freedom real system. 

Now the equations of equilibrium for a fluttering wing 
may be written: 


ju + ( how M )] — vat Ea (A pa/ Je _ 0 
\ (o,/ Wa)” d ated ( h We)? 


An A ah A aa r e 
Sp Oh tt + —™ — X(1 + fg.) la = 0 
‘ we uVa? 


where 
X = (w,/w?), the square of the ratios of the 
torsion frequency to the flutter fre- 
quency 
(w,/#,) = the ratios of the bending frequency to the 


torsion frequency 


im = ratio of the wing weight to that of, 
cylinder of air of diameter equal to th 
chord 

ae = distance from elastic axis to center 


gravity in terms of the half chord 


Ve = radius of gyration in terms of the haj 
chord 
h = dimensionless bending amplitude (coy 
plex) 
a = torsion amplitude (complex) 
j = V/-1 
La = torsion structural damping constant 


and the A’s represent aerodynamic force coefficien: 
that are complex functions of the reduced velocit 
(V/bw), where V is air velocity, 6 is the half-chord, an 
w is the flutter frequency. 

In the following example, we put » = 10, (a,/wa)=|+ 
4/10; 4. = 1/5, Va= 1/+/5, 980 that the equations 
become 


(10 + A, — X)k + (2 + A,,)a = 0 


po” Aa. ’ ; 
(1 + <= )n + [ ~~ =~ ou + it) Ja = 


Since the A’s depend on (V’/bw) in a rather complicated 
fashion, the equations will be solved for various as- 
sumed values of this quantity.* A table of the A’sis 
shown here as Table 1. 

Now these two complex equations may be put on the 
machine as four real equations. For example, with 
(V/bw) = 1.00, 


(10.80 — x)h® + 1.08h’ + 1.00a? + 1.5807 =0 
(10.80 — x)h’ — 1.08h? + 1.000’ — 1.58a® =0 
1.14h® — 0.14h/ + (1.25 — X)a® + (0.28 + Xg,)a’ 


1.14h’ + 0.14h* + (1.25 — X)a’ — (0.28 + Xqg,)a' 


where h = h® + jh’ = a® + ja’. 

Some further modifications are desirable to adapt the 
problem to the machine. For (V/dw) = 1 and 2, the 
equations are divided by 6 and a new “X”’ term, X = 


(X — 5)/6, is used. 


Thus, for example, the equations for (I’/bw) = 1 are 
(0.97 — Y)h® + 0.18h/ + 0.17 = 0 
—0.18h* *. (0.97 — X)h’ — 0.26 = 0 


0.19h* — 0.0247 — 0.63 — X = 0 
0.02h® + 0.19h’ — 0.05 (—Xg,/6) = 0 


In the above a has been taken as 1 + 0, since one of 
the amplitudes is arbitrary. 

For the cases (V/bw) = 3 and 4, the first equation 
was divided by 50, the second by 10, and the first 
column was multiplied by 5. X was taken as(X —1)/10: 


* These are for an assumed elastic axis at 38.6 per cent chord 
to fit conveniently available aerodynamic terms. 
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These equations are then solved in the machine for h, 
X, and a quantity set in the last constant dial, propor- 
tional to —Xg,, and called Y in Table 2. 


ACCURACY IMPROVEMENT 


The accuracy of a machine solution may be extended 
to any degree desired by calculating the errors exactly 
and solving for the increments in the unknowns to an 
enlarged scale. Suppose, for example, that the equa- 
tions are written in matrix form AX = B, where A isa 
square matrix of the coefficients A7j and X and B are 
column matrices of the unknowns and constants. If 


X = X,+ AX 


where X, is an approximate solution, then the column 
matrix of the errors e will be given by 


AX,-—B=e 
and since, by definition, 


A(X, + AX) -B=0, AAX =-e- 


Now the e’s and AX’’s will be small so that the scale 
may be changed to use the full machine capacity. It 
should be noted that the setting of the coefficients in the 
machine need not be altered in this process; all that is 
required is an accurate computation of the errors. 

The above process may be applied also to matrices 
involving a \ term. Thus, if \; and X, are approxima- 
tions to a solution of (A — AJ)X = 0, where J is the unit 
matrix, then 


[A — (4. + AA)](X, + AX) = 0 
or 


(A — v,J)X, + (A — AJ) AX — AAX, — AAX = 0 


Now (A — ,/)X, = e, the matrix of the errors, and 


if these are computed exactly and if, furthermore, the 
second order error AX AX is neglected, 


(A — AJ)AX — AX, —<€ 
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i TABLE 1 
_— (V/bw) Abn Ana Aa, Aaa 
1 0. 80-1. 08j —1.00-1.58j 0.28 + 0.28) 0.50 — 0.56 
2 0.40-2.397 —4.99-2.537 0.20 + 0.647 1.60 — 1.307 
3 —0.05-3 .907 —12.23-2.69) 0.50 + 1.067 3.56 — 2.26) 
4 —0.48-5.547 —23.01-2.11j 0.64 + 1.507 6.50 — 3.427 
-_ TABLE 2 
Results for the Flutter Problem 
(V/bw) h = Y 4 La V, 
M.P.H. 
1 —0.13 + 0.13) —0.65 +0.02 3 —0.11 192 
2 +0.23 + 0.337 —0.52 0.02 1.88 —0.064 285 
3 0.95 + 0 75] +0.265 —0.06 3.65 +0.164 322 
4 1.65 + 2.25) +0.38 —().27 4.8 +0.562 375 
o5V 
Where b = 30 in.; wa = 1,125 cycles per min.; V = 5.95 X 30 1.1251 bw (See also Fig. 11.) 


VX 


Now as one of the X’s is arbitrary, the corresponding 
error may be taken to be exactly zero, so that the  un- 
known increments will be 


AX,, AX2.... AXg-1, AA 


Thus, if the last column of the A matrix is replaced by 
the column matrix X, and the matrix of constants 
(originally set zero) is replaced by —e, a solution for the 
AX’s and the Ad is obtained by solving the above 
ordinary simultaneous equations to an enlarged scale. 
If an extra column of coefficients were provided, the 
AXX, could be set in without changing the (A — XJ) 
matrix. 


In summary, by means of the above methods, a 
solution of any desired degree of accuracy may 
be obtained. For example, if the original solution were 
accurate to two figures, the first improvement would be 
good to four figures, the next to six, etc. 


This method of accuracy extension also provides an in- 
dependent check of the calculations, including machine 
errors, dial setting, and reading errors. 


It is possible to devise systems of equations which are 
nearly singular, and any method of computation will 
have difficulty with these. Experience to date has indi- 
cated that engineering problems do not, in general, fall 
in this class. If a problem is near a singularity, it is 
sufficient to know this, since any solution that de- 
pends on a critical combination of parameters is usually 


not satisfactory from an engineering standpoint. 


TIMESAVING 


The sample problems worked above were done both 
on the simultaneous equation machine and on desk 
calculators with experienced operators. A timesaving 
of 3 to 1 was shown by the equation machine. It is be- 
lieved that this will be even more favorable on larger 


order problems. The labor involved in solution on the 


(Continued on page 127) 








A Note on the Approximate Plane Motion 
During the Burning Period of Rocket- 
Propelled Missiles Launched at Small 

Angles of Yaw from Aircraft’ 


RAY E. BOLZt 
Rensselaer Polytechnic Institute 


SUMMARY 


The differential equations of plane motion during the burning 
period of a rocket-propelled missile launched from a moving air- 
plane are written, simplified for approximate application to small 
angles of yaw (20° or less), and solved. In writing the equations 
of motion, the gravitational force was neglected, as well as the 
influence of any pitching disturbances on the yawing forces or 
moments. 

The solution determines the approximate motion ot the center 
of gravity and the spatial orientation of the missile axis as a func 
tion of missile design, the initial launching conditions, and the 
distance along the trajectory. The results are plotted in general 
ized coordinates, and the curves may then be used for the solution 


of specific cases. 


INTRODUCTION 


Ws A ROCKET-PROPELLED MISSILE is launched 
from a moving airplane in any direction other 
than in the direction of flight, the missile is subjected 
to the aerodynamic forces arising from the resultant 
angle of yaw. This angle of yaw, which is the angle 
between the missile axis and the relative wind vector, 
has in the general case a magnitude within the range 
of 0° to 180° depending upon the angle of launching, 
the aircraft velocity, and the missile velocity at launch 
ing. 

In this paper, the differential equations of motion 
during the burning period of a rocket-propelled missile 
launched from a moving airplane are written, simpli 
fied for the specific approximate application to small 
angles of yaw (say 20° or less), and solved. In writing 
the equations of motion the gravitational force, the 
influence of any pitching disturbances on the yawing, 
and the change in mass of the missile during the burning 
period are all neglected, and all motion is assumed to 
take place in the plane determined by the vectors of air 
plane velocity and missile velocity. 
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others—-have attacked the problem in various ways 
This paper merely presents a somewhat different and 
rather simple solution of an approximate form of the 
equations of motion. 

The solution derived herein determines the motion oj 
the center of gravity and the spatial orientation of the 
missile axis as a function of missile design, the initial 
launching conditions, and the distance (or time) along 
the trajectory, and it is approximately valid for the 
range of yaw angles up to about 20 The results are 
plotted in generalized coordinates and the curves may 
then be used for the solution of specific cases. 


SYMBOLS 
a = constant, (A /2m) + (C,/2/) 
b = constant, i—(C/I) — [(K/2m) — (C,/2] 
ri = wing chord 
“8 = constant, (dC,,/da) (pSc/2) (negative for stable 
missiles) 
C = constant, [dC,,/d(Bc/2V)|(pSc?/4) (always posi 
tive 
Cr, Cy, = aerodynamic coefficients of lift and moment 
D = drag 
F = force 
K = constant, (dC,/da) (p.S/2 
l = bz 
'g = lift 
M = moment 
m = mass 
S = characteristic area, usually wing area 
1 = thrust 
t = time 
V = scalar velocity of missile 
V = launching velocity of missile with respect to 
ground 
V = vector velocity of missile 
= distance along trajectory 
a = angle of yaw, between missile axis and _ relative 
wind 
8 = angle between missile axis and y-axis 
+ = constant of integration, y = Vo?m/27 
\ = by 
rr angle between relative velocity vector and y-axts 


EQUATIONS OF MOTION 
Consider a coordinate system in which the orthogonal! 
axes x and y are located in space in the plane determined 


by the vectors of the airplane and missile velocities such 


that 1 
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that the y-axis is parallel to the airplane velocity vector 
at the instant of firing. As shown in Fig. 1, let V be 
the velocity of the center of gravity of the missile and 
the vector V is then tangent to the trajectory at every 
instant. Let @ be the angle between V and the y-axis, 
3the angle between the missile axis and the y-axis, and 
athe angle of yaw. All angles are taken positive in 
the clockwise direction. 

Now since the gravity forces are neglected the motion 
of the missile will take place in the x-y plane. The 
equations of motion are then written with respect to 
axes fixed in space and the linear motion is resolved 
along axes parallel and normal to the relative velocity 
vector at every instant. 

The equations of motion are: 


mV = >F (la) 
Ig = >M (1b) 


The boldface italicized symbols signify a vector. 
Since V = vl’, then 
Vv = Vo,+ Vo, = Vov, + V», (2) 
where v, and v, are unit vectors parallel and normal to 
the vector V. 


Eqs. (1) now become: 


mV = Tcosa —D (3a) 
mV =Tsina+L (3b) 
Ig =M (3c) 


If the equations are restricted to small angles of 
yaw, the two following assumptions may be made: 
(1) Sina=a; cosa=1. 
(2) The drag is negligible compared with 7. 
The lift, L, and moment, J/, are written as: 
L = KV’a 
M = CV’a — GVB 


where the coefficients 


ll 


K (dC,/da) (pS/2) 
C = (4C,/da) (95c/2) 
C, —|{ [dCn d(Bc/2V )}} (p.Sc*/4) 


are considered as constants. It will be noted that in 
the equations above, K and C, are always positive 
numbers whereas C is negative for a statically stable 
missile and positive for an unstable one. 

The approximate equations of motion finally may be 
written : 


mV=T (4a) 
mVo = (T+ KV*)a (4b) 
I(¢ + &) = CV2a — CV(o + &) (4c) 


where 8 = @ + a and a dot over the symbol signifies 
differentiation with respect to time 
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Coorpinate System 
Fic. 1. Schematic diagram illustrating the geometry used for the 
analysis. 


SOLUTION OF EQUATIONS OF MOTION 


Eq. (4b) divided by V and differentiated gives: 


mo = (T+ KV2)[(a/V) — (aV/V2)] + 2KVa (5 


Substitution of Eqs. (5) and (4b) into Eq. (4c) gives: 


m | V? 
sod atc 
fa = Cla — CiI E + (* + . ) a | (6) 
m\ 


From Eq. (4a) 


I (=~ of aV7 + KVa + Kal’) + 


V = T/m (7) 
and 
d (>) _ T — T dz 
dt \2 m m dt 
V2? = (27/m) (2 + y) (8) 
where 
y=Lotu= Vo2m/2T (Sa) 


Here Ly is the length of the launcher and uy is a con- 
stant to account for the part of the launching velocity 
of the missile contributed by the airplane velocity, 
That is, atz = 0, 

Vo? = (2T/m)y (9) 

The independent variable, /, is now changed to the 
distance, z, along the trajectory of the missile. 

Then in the usual manner 


d/dt = V (d/dz) (10) 


and 


d?/dt? = V(d?/dz*) + V (dV/dz) (d/dz) (11) 
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) 4 12 16 20 
R= bz 
Fic. 2. The dimensionless quantity ¢ — go/ ao ( proportional to 


the change in the angle between the missile velocity vector and 
the y-axis) plotted as a function of the dimensionless number 
1 = bz (proportional to the distance along the trajectory) for the 
value a/b = O and various values of \. Here a/b is a number de- 
pendent upon missile design and \ is a number dependent upon 
the initial launching conditions. 


The use of Eqs. (7), (8), and (10) with Eq. (6) re- 
sults in 


a” a! ( aa + . +> “) + 
sty m I 
a( K _ I an C 4 
2m(z + vy) 4(2 + vy)? I 
C <*) — 


- (12) 
2I(2 + vy) 


Im 


where a’ = da/dz and a” = d’a/dz*. From Egs. 


(4b) and (8) then, 
¢’ = {[1/2(¢ + y)] + (K/m)}a 
Now let z + y = Xin Eq. (12). 
a ta’(s + : + ot) + 
Xx, m I 
( K l t C 
oo = = 
2mX, 1X? I 21X 


(13) 


Then, 


ar =( (14) 
Im 


If 


a ee . 
a = e **X,—* cos bX,* 


(15) 


* The general torm of this solution was suggested by Dr. D. P. 
Ling, of the Bell Telephone Laboratories. 


AERONAUTICAL 


SCIENCES FEBRUARY, 1950 

is the assumed form of the solution, an insertion 9 

Eq. (15) into Eq. (14) simplifies to the following gy. 

pression : 

a?X,~ '/* cos bX, + 2abX,~"”* sin bX, — b2X,7' x 
— eae — 

cosbhxX, — Ky “eostxz,. — —Z, 

m m 


Cy mite 


‘sin bX, — 


Cia 

I 
C X,~ '” cos bX; + GK Xi~"” cos bX; = 0 (Ig 
I Im 


X,~” cos bX; — sin bX, — 


If the two terms involving the sin bX, function ar 
eliminated by proper choice of a, the remaining equa. 


tion can be solved for 6. Obviously a must equal 


a = K/2m + G/21 (17 
Eq. (16) becomes 


b= V —-(C/I) — [(K/2m) — (C,/21)]? (18 

The constant C is negative for a statically stable 
missile and an inspection of 6 shows that C/J > 
[(K/2m) — (C,/27)]* for usual designs. 
however, breaks down for a missile designed for nearly 


The solution, 


neutral static stability or instability. . 





2. aS es eee ee | | 
oO 4 8 12 4 20 
&= bz 
Fic. 3. The dimensionless quantity @ — qo/ao as a function 
of the dimensionless number / = bz for the value a/b = 0.15 and 


various values of X. 
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0 AIRPLANE MOTION DURI 
Sertion of Thus Eq. (15) is a solution of Eq. (12) for values of a 
Wing ey. | and b defined by Eqs. (17) and (18). Similarly 
a =e ™X,~'” sin bX, 
x isalso a solution and therefore the general solution may 
be written : 
- 
p= ett”) (5 $+ y)7"* {A sin [be + y)] + 
B cos [b(z + y)]} (19) 
From Eq. (13) 
0 (1 ¢’ = {[1/2(s + y)] + (K/m)}a 
(if 
fyvaluation of the Constants A and B from Boundary 
‘tion are 9} Conditions 
1g equa: At the exit of the launcher, z = 0, the missile has no 
lal yawing angular velocity. 
(li) Fats =0 
p’ = By’ = 0 
a = ao 
(18 & = do 
Since 
’ stable 
"ts I » B ! = ron) ss 
ution, a= 
nearly § From Eq. (13) 
go’ = [(1/27) + (K/m)] ap = — ay’ (20) 
> And from Eq. (19) 
a! = (—ae~*%y~"/? — (1/2)e~*7 y~*”) (A sin by + 
—+ Boos by) + ey” (Ab cos by — Bb sin by) (21) 
And also 
ia a = e*%y~'(A sin by + B cos by) (22) 
z Combination of Eqs. (20), (21), and (22) results in the 
evaluation of A and Bas: 
| A = (ay/b) e*’y'” (b sin by — a cos by) l (23) 
* B = (a/b) ey” (b cos by + a sin by) f a 
- SOLUTION FOR ANGLE @ 
+ Returning to Eqs. (13) and (19) 
+] ’ = {[1/2(¢ + y)] + (K/m)}e~*@t” (¢ + y)7”" X 
~ 1A sin [b(z + y)] + Boos [b(2 + y)]} (24) 
| Integrating from zero to z, 
F t= om toala f eM tVig 4 vy)? &X 
0 
[ sin [big + y)] dz +B f e~ tV(g +)" X 
0 
t 
cos [b(z + y)]dz} + (1/2) 4 ff g Ser 
0 
a: Zz 
/ ty) “sin [b(s + y)]dz + Bf grr" x 
ction 0 
) and 8/. . a. 
(s + vy)” “*cos [b(z + y)]dz} (25) 
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If the last two terms are integrated by parts once and 
combined with the first terms, the equation for ¢ be- 
comes 


—@ =m — et 4 4)" x 
{A sin [b(s + y)] + B cos [b(z + y)]} + 
I(aA — bB)+1{aB + bA) (26) 


where A and B are defined by Eq. (23) and 


I, -f{ ety (g + ¥)~'”* sin [o(s +») ds} 
0 


‘ (27) 
I, = f e “Ty (g + y) 2 cos [b(z + y)|as| 
0 
Now let 

bz = 1 

and 
by =A (28) 

then 
dz = dl/b 


If Eqs. (28) are inserted into Eqs. (26) and (27) and 
ao factored out: 
(6 — ¢o/ao) = 1 —e "FM + )~* XK 
[A’ sin (J + A)+ B’ cos (1 + d)] + 
I,{(a/b) A’ — B’]) + I,[(a/b) B' + A’] (29) 


where 


A’ = e'/*)"* [sin X — (a/b) cos d] t 
B’ = e@/* 2 [cos \ + (a/b) sin A] f 


> 


and 


l 
I. = f e @/EF™ f+ y)—" sin (I + rat) 
0 
(31) 


l > 
I, = f e~ e/b)+a) (1 + )~'” cos (i+ ryal| 
0 

Thus Eq. (29) shows that the change in the angle of 
the missile velocity vector A¢/ap is a function of three 
dimensionless variables—namely: / = 62, propor- 
tional to distance z along the trajectory; A = by, a 
function of the initial launching conditions; and a/b, 
a function of missile design. 


Here 
y=Lotpu = Vo2m/2T (Sa) 
a = (K/2m) + (C,/2/J) (17) 
b = V -(C/D) — [(K/2m) — (C,/21)]? (18) 


Since 8B = @+a 


EF wk F A’ - B'| +I. ¢ B’+A | (32) 
ao b b 


Thus A/a is similarly dependent upon /, 4, and a/d. 
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=30 Fic. 5. The dimensionless quantity 8 — Bo/ao (proportiona 
| — | | Thee to the change in the angle between the missile axis and the y-axis 
8 12 16 fe) plotted as a function of / = bs for the values of a/b = t 
Q: bz various values of X. 
Fic. 4. The dimensionless quantity @¢ — @o/ao as a function 
of the dimensionless number / = bz for the value a/b = 0.30 and 
various values of X. 
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DISCUSSION 


Presentation of the Calculated Results 


The dimensionless quantity (¢ — ¢9)/ao which is 
proportional to the change in the angle of the missile 
velocity vector from its initial value, is plotted in Figs. 
9 3, and 4 against the dimensionless number bz, pro 
portional to the distance along the trajectory, for vari 
ous values of the dimensionless parameter a/b. Figs. 
5, 6, and 7 show the change in the orientation of the 
missile axis and Figs. 8, 9, and 10 show the respective 
changes in angle of yaw all plotted in an identical 
manner. The parameter » [see Eq. (28)] is a function 
of the absolute launching velocity of the missile and 
therefore of the launcher length, airplane velocity, 
launching direction, and the thrust/mass ratio of the 
missile. The parameter a/b and the proportionality 
factor 6 [Eqs. (17) and (18)] are functions of the missile 
design only. 

The design of a specific missile chosen for study then 
fixes the values of a 6 and 6; and the launching velocity 
or launching conditions) fixes the value of the param 
eter \. The particular (¢ 
3, and 4+ and the particular (8 — 6»)/ap curve in Figs. 
5, 6, and 7 corresponding to the a/b and X previously 
determined show, respectively, the angle of the missile 
velocity vector (and therefore the slope of the trajectory 
curve) and the orientation of the missile axis at any 
The curves 


goo) /ao curve in Fig. 2, 


point along the missile flight trajectory. 
presented then give the resultant plane motion during 
the burning period of a// missiles included within the 
scope of the parametric values chosen and for launching 
conditions which satisfy the approximations made in 
the equations of motion. 

It may be noted from the equations that the wave 


length of yaw ois 


o = 2r/b 


which, for usual missiles where C// > [(K/2m 
C,/27)]*, agrees with the usual definition of o given 
i reference 1. The abscissas for the curves, / bz, 


can then be written 
l 22/0 


in terms of this wave-length factor. 


Approximations 


The approximations made during the analysis are (1) 
the coefficients of force and moment are constant; (2) 
gravity is neglected; (3) rocket thrust is constant and 
large compared with the drag forces; (4) the change in 
lass of the missile during the burning period is neg 
lected; and (5) the angle of vaw is restricted to values 
Where sin a a@and cos a i. 

The first approximation cannot be fulfilled when the 
. 


velocity of the missile attains a very high subsonic « 


Supersonic value except by use of a rough mean value 
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Fic. 8. The dimensionless quantity a — ay/ay (proportional 
to the change in the angle of yaw) as a function of / = bs for the 


value a/b (and various values of \ 
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Fic. 10. The dimensionless quantity @ — ay/ay as a function 
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for the coeflicients or by a step solution. However, 
as noted in reference |, the variations in the coefficients 
will occur in the latter part of the trajectory (where 
the high missile velocities are attained) and therefore 
will have small effect on the trajectory itself. The 
first approximation also breaks down when the angle of 
yaw over the initial part of the trajectory is large enough 
to cause This latter 


condition was investigated in the Rensselaer Poly 


variations in the coefficients. 
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technic Institute subsonic wind tunnel for a general 
supersonic missile configuration and the results are 
shown in Fig. 11. The lift and center of pressure 
(c.p.) curves show an approximately linear relation for 
the lift and moment coefficients with angle of yaw up 
to about 30° or 40°. 


Brief Description of the Resulting Motion 


The motion of the missile is, of course, a result of 
the missile design, the wind velocity, and the initial 
angle of yaw. Upon launching, the moment of force 
acting on a statically stable missile rotates the missile 
very quickly toward the direction of the realtive wind 
vector. At the same time the acceleration of the missile 
along its own axis resulting from the rocket thrust and 
the velocity induced by the lift force cause the relative 
wind vector to swing around toward the missile axis. 
The relative magnitude of these two effects, as well as 
of the damping moment, then determines the missile 
trajectory as well as the orientation of the missile 
axis throughout the rocket burning period. 

If the missile, however, were very nearly neutrally 
stable, then almost no rotation of the missile axis 
would take place and the wind vector would rapidly 
align itself with the missile axis resulting in little angu- 
lar deviation from the launching condition. 


Discussion of Calculated Results 


The values of a/b chosen for the calculations are ex- 
treme and were chosen to illustrate the general char- 
acteristics of the solution. The value a/b = 0 has no 
physical significance but is of interest as a limiting 
value only. 

From a comparison of the curves and an inspection 
of Eq. (29) it is seen that the lift coefficient term (A /2m) 
acts as a damping factor in almost identically the same 
manner as the damping moment term (C,/2/) for stat- 
ically stable missiles where C/J > [(K/2m) — 
(C,/27)]? and thus 6 & V C/I. However, for missiles 
with very small static stability, the value of b is appre- 
ciably affected by the value of [(AK/2m) — (C,/2/)] 
while a = K/2m + C,/2/. Since both a and 6b are 
damping factors, the a in the exponential term and 
the db in the (J + \)~”* term, the damping action of the 
lift force differs from that of the damping moment for 
this particular case. 

The physical action of the lift force in producing 
damping results from the fact that lift induces a ve- 
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center of pressure (c.p.) measured in the R.P.I. tunnel as a fune. 
tion of angle of yaw from 0° to 180°. The missile consisted of a 
cone-cylinder fuselage 25 in. in length and with cruciform wings 
and tail surfaces attached and was tested at 150 ft. per sec. air 
velocity. 


locity component normal to the trajectory in a diree- 
tion which always aids the relative wind vector to 
swing around into the path of the missile, thereby 
decreasing the angle that the missile will rotate in its 
simultaneous attempt to swing around into the wind 
vector. Any factor causing a decrease in the angular 
displacement of the missile appears in the resultant 
motion as a damping factor. 

The relative movement of the wind vector and the 
missile axis is given by a comparison of the A@/ ap and 
the AB/ao curves for any missile within the range of 
variables chosen for calculation. The damping effec 
tiveness is observed in the A8/a curves. 
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Seventy-Fifth Anniversary 


of 


Ludwig Prandtl 


HEN THE ART OF FLYING began its modern growth 
W: the early days of the century, little was known 
of what today we call aerodynamics. Thoughtful ob- 
servation of birds and model experiments with wings 
played an important part in the minds of daring in- 
Today everyone is aware that the production 
It is well known that, 


ventors. 
of airplanes is a major industry. 
in the rapid development of this new branch of engi- 
neering, scientific research played an important part. 
One of the outstanding figures in this field is Ludwig 
Prandtl in G6ttingen, who on February 4, 1950, will 
complete his 75th year. In recognition of his world- 
wide influence and prominent leadership, a survey of 
his scientific work is here presented to our readers. 
Long before governments recognized the importance 
of airplanes or expended any funds on their develop- 
ment, Prandtl] began his studies of airflow and threw 
light on one of the darkest corners of applied mechanics. 
At the International Mathematics Congress in Heidel- 
berg in 1904, he presented the famous paper on the 
motion of fluids with extremely small viscosity, which 
laid the foundation of boundary-layer theory. It gave 
the explanation of the great discrepancy between the 
theory of potential flow and observed facts and showed 
how and why this discrepancy could be avoided in cases 
of aeronautical importance. Much later, when air- 
craft had developed so far that extreme reduction of 
drag became a question of vital importance, Prandtl’s 
boundary-layer theory became a subject of world-wide 
mathematical investigation. This investigation still 
continues 
Soon after his publication on boundary layers, 
Prandtl developed an active interest in airships, .which 
led to the creation of Germany’s first wind tunnel. 
The design of this tunnel is visibly influenced by 
Prandtl’s conception of the boundary layer: The free 
jet avoids all wall friction in the test section, and the 
convergent collector reduces the inhomogeneity of the 
air stream because of the boundary layer developed in 
the duct. The laboratory that developed around this 
tunnel ranked among the first of its kind and became one 
of the leading centers of aerodynamic research. Not 
only did it yield many valuable experimental and 
theoretical results, but it also turned out a considerable 
number of excellent research workers who spread all 
over the world and brought with them that Gottingen 
spirit of applied research which so happily connects 
pure science with the practical problems of the engineer. 
The second of the great impulses which Prandtl gave 
to aeronautical science developed in the secrecy of a 
war period. Published first in restricted form in 1918, 
then in 1919 in the Nachrichten der Géttinger Gesellschaft 


121 





Ludwig Prandtl 


der Wissenschaften, his papers on airfoil theory soon 
became the base of research activities in all countries 
where airplanes are built. It is well known how far 
the shape of airplanes was changed by this research. 

Prandtl’s publications on both subjects, boundary 
layer and wing theory, demonstrate his specific way of 
scientific work. He shows surprising relations, pre- 
viously overlooked or incompletely understood, gives 
the basic ideas for their analysis, and so inspires many 
other workers to continue the mathematical develop- 
ment and work out the technical applications. This 
has been characteristic of him from his Doctor's thesis 
to the present day. 

In the nineteen-twenties, public interest in aviation 
increased enormously. The strange flying machines 
of the prewar period were developing into the new 
streamlined airplanes. The war having shown beyond 
doubt their value for destructive purposes, the ques- 
tion of peacetime use presented itself. The days of the 
pioneers had ended, and commercial aviation began. 
In these years, Prandtl’s position in aeronautical re- 
search also changed from that of the scientific trail- 
blazer to that of one of the most prominent men in a 
world-wide family of scientists. His further publica- 
tions are no longer explorer’s trips into new lands 
but are leading contributions to questions of great 
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interest which find an immediate response beyond the 
circle of his personal followers. 

Thus he participates in the investigation of the turbu- 
lence problem, seeing its full breadth from the velocity 
distribution in pipes to the turbulence of the free 
atmosphere. His particular interest is devoted to the 
transition from laminar to turbulent flow in pipes and 
along walls, and he initiates theoretical investigations 
on the stability of the flow in the boundary layer. The 
mixing zone on the surface of free jets also attracts his 
attention, and the notion of the mixing length, formu- 
lated by analogy with the free-path length in the 
kinetic theory of gases, proved to be a new and useful 
tool in this field. 

As early as 1904, Prandtl had published papers on 
the dynamics of compressible fluids. The topic was 
resumed in more recent times, when airplanes reached 
such speed that the compressibility of air became of 
substantial importance for their performance. 

Less known to aeronautical engineers is Prandtl’s 
work outside aerodynamics, and many among the 
younger generation may not know that he began his 
scientific life with a Doctor’s thesis on a problem of 
elastic stability. This part of the theory of elasticity, 
today so vital in aircraft design, was rather undeveloped 
at the end of the last century. In his thesis, published 
in 1900, Prandtl presented and solved a new problem 
that of the lateral instability of beams in bending. 
More than that, in his experiments that accompanied 
the theory, he conceived and developed the method of 
adjusting the test specimen so that the influence of its 
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imperfections on the critical load is compensates 
now a standard practice in all careful buckling ey 
periments. 

The second of Prandtl’s early papers on ¢lag 
problems is known to a larger public. In 1904 }, 
discovered the soap-film analogy for the Saint-Vengy, 
torsion problem. Because of its great suggestivenes 
it became one of the standard topics of textbooks 
dealing with torsion, and its practical value for streg 
analysis has been recognized by many engineers, 

In 1920 his first paper on the theory of plasticity 
peared, and, subsequently, we find the subject appearing 
at intervals in his publication list in contributions ; 
knowledge which have stimulated interest in the subjeq 
in G6ttingen and elsewhere. 

This sketch of Prandtl’s scientific activities, neces. 
sarily far from exhaustive, is not enough to deserih 
his influence on present-day applied mechanics. Pe 
haps an even more important aspect of his activity js 
that of an academic teacher of the first rank, gathering 
through almost half a century an ever-changing grou 
of pupils about him. These young men (and not onl 
young ones) acquired from him both scientific know. 
edge and a scientific attitude. They continue the work 
so begun in many parts of the world—not a few of them 
in this country—contributing to the training of a new 
generation of research workers and engineers. Their 
thoughts will be with Prandtl on his 75th anniversary, 
and they will wish to join with fellow scientists the 
world over in saluting him as one of the greatest 
figures of engineering science. 





Discussion of Wright Brothers Lecture 


(Continued from page 106) 


accurate measurement of the free-body vibration 
modes. In our own thinking on this problem, we 
have come to regard the ground vibration test merely 
as a means for verification of the mass and stiffness 
distributions and of the method employed to calculate 
vibration modes. If this philosophy is adopted, then 
there is no real necessity for simulating the free-body 
condition, since comparisons may be made between 
measured and calculated modes for a rigidly supported 
structure. 

Finally, with regard to flight vibration tests, it is 
contended that the use of a vibration exciter does not 


in any way increase the hazard of the tests. Of course, 
it is essential to differentiate between vibration testing 
and flutter testing in flight. If, as implied by Mr. 
Russell, the test is conducted solely for the determina- 
tion of vibration modes and natural frequencies of the 
airplane in free flight, then it must be assumed that the 
airplane will not be operated deliberately near 4 
critical flutter speed. If such a deplorable situation 
should occur, then the vibration exciter at least gives 
some assurance that the critical speed will not be ex 
ceeded by a large margin, thereby reducing the likelt- 
hood of catastrophic flutter. 
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Note on Rotational Gas Flow 
james E. Broadwell 
Ann Arbor, Mich. 


October 10, 1949 


SYMBOLS 


H angular momentum about center of mass 

‘ unit vector perpendicular to plane of flow 

L = torque about center of mass 

m = mass 

p pressure 

P vector to center of mass, with rectangular components x and yi: 
, = velocity, with rectangular components wu and 7 

# = distance from origin with rectangular components x and y 

5 area 

p density 


_ LAWS OF conservation of mass and momentum for the 
two-dimensional, rotational, diabatic flow of a perfect com 


pressible fluid may be put in the form: 


d 
(Vp!) X (Vp) = —p (p~! w) (1)! 

dt 

where 
Ov Ou 
Ps om 

Ox oy 
md the symbol d/dt indicates differentiation following the 
motion of the fluid. In the following, Eq. (1) is derived by 


equating the torque which acts on a fluid particle, assumed to be 
square, to the rate of change of angular momentum of that par 
ticle. Such 
fuid mechanics where the concept of rotation causes some diffi 


a derivation would perhaps be useful in teaching 


culty. 
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Consider a square fluid element of unit thickness at the origin 
For any element of the fluid, 


dl 
dt 


Al ff le - B x doas 


For the above element, 


— a /2 ea /2 
H=«x J G-ale- G- y)ulpdxdy (3) 
—a/2J —a/2 z 


Let the particle be of such size that the density and velocities 


as shown in Fig. 1. 


(2) 


where 


inside are given by the equations 


Op Op 
p = po + x y 
Ox/o Oyo 
(4) 
Ou Ou 
u= uw +t x + y 
ovo oy/i 


and likewise for v, where the subscript 0 denotes values at the 


The center of mass P is defined by the equation 


in 1 
P =a ff 7 dm 
WU JS 


where ./ is the total mass of the element. 


origin 


With Eqs. (4) the 


components of P become 


l (°”) a’ 
x, = 

Po Ox /o 12 

l 2) a’ (9) 
"= 
is Po Oy o 12 


Substitution of Eqs. (4) and (5) in Eq. (3) and integration result 
in 
” , wo a” . 
H = kxpoa? : (6) 
2 6 


if terms of order a‘ and higher are dropped. Note that this is the 
expression for the angular momentum, measured about the center, 
of a square of uniform density and angular velocity wo/2 

The torque, which results only from pressure forces, is given by 


9 


ea/2 a/2 
K (y — w)pdy 
e/2 ~ a/2 
Pa /2 
/ 
J a 


a/2 
(x — %:)p ax | (7) 
2 a/2 


L 


Setting 


| *) ‘ ) 
= ” + x y 
P P Ox Jo Oy /o 
substituting from Eq. (5) and integrating, changes Eq. (7) to 
- Op Op Op Op a* 
L=-—k = . (8) 
Ov) i\Ox/i Ox] )\Oy/o Jl 2p, 


Eqs. (6) and (8) combine through Eq. (2) to 


[EN@)-Ce-se] 


Since Eq. (9) is to be applied to a particle as it moves along the 


d 
dt 


Op a’ 


ox 


op 
Oy 


re) p 
Oy 


op 


Ox L2p; 


streamline, the mass, pov? to the order of Eq. (9), must be held 


He) - CD) -sle 


or, dropping the subscript, to Eq. (1), 


constant 


1 
dt 


Op 
Ox 


Op 
oy 


Op 
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op 
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A ‘ 
-/ 
v/ 
oP 
Fic. 2. 


(Vp!) X (Vp) = ut (pw) 

dt 

From the above derivation it is seen that the quantity (Vp~!) 
(V p) is proportional to the torque acting on the element and that 
wp! is proportional to the angular momentum. In the simple 
case of a circular disc of compressible matter held at uniform 
density, p, it is easily shown that the angular momentum is pro- 
portional to Qo~! where { is the angular velocity. 
disc of radius R, 


For a circular 


H = (1/2)rR%pQ 


If the body expands radially and if the mass, prR?, remains 
constant, 
1 (prR?)? Q Q 


H = - = constant 
2 Tv p p 


The physical reason for the angular momentum remaining 
constant when p is a function only of p (so that (Vp~') X Vp 
is zero) is pointed out by F. H. Clauser in a Douglas Company 


report.2, Consider again a circular element in a flow field (see 
Fig. 2). If pis a function p only, then the center of mass lies 


somewhere on the dotted line parallel to Ap and the torque 
about this point is zero. 
independent of p as in diabatic flow. 


This is not the case, of course, if p is 


REFERENCES 
1 Vazsonyi, A., On Rotational Gas Flow, Quart. of Applied Math., Vol. ITI, 
No. 1, April, 1945 


2 Clauser, F. H., Theory of Compressible Fluids, Douglas Aircraft Company 
Report, April 8, 1944. 


Errata—The Aerodynamic Forces on an 
Oscillating Flap at Supersonic Speeds 


John W. Miles 
Department of Engineering, University of California at Los Angeles 
September 26, 1949 


eS TAYLOR of Douglas Aircraft, Inc., has pointed out to 
us that Eq. (22) and Figs. 8 and 9 in this paper were cor- 
rect as originally stated and not in error, as stated in a previous 
errata (JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 16, No. 
7, p. 442, July, 1949) 
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The Mechanism of Unstable Combustion 


W. C. Randels 
Northrop Aircraft, Inc., Hawthorne, Calif. 
November 4, 1949 


HE PHENOMENA OF unsteady burning in forced-draft op 
bustion systems has been frequently observed. The hi 
pose of this note is to show that the possibility of such OCCUrrences 
might be expected on the basis of a conventional stability analy,, 
Consider a system consisting of a compressor, combustig, 
chamber and nozzle as in Fig. 1. The performance of the eg, 
pressor can be expressed by the ‘‘characteristic’’ curve of presgyp 
as a function of mass flow as Curve A in Fig. 2. Similarly th 
total pressure required by the combustion chamber can } 
graphed, and the intersection with the compressor characteris} 
gives the operating point. Two types of intersection are jl}, 
trated by Curves B and C in Fig. 2. 
these types can be explained as follows: Consider Curves A ayj 
B in the vicinity of the point of intersection in Fig. 3. If for any 
reason an increase in flow changes the operating point to Pain 


The difference betwee 


1 at some instant, then at that time the pressure supplied by th 
compressor will be less than that required by the combustigy 
chamber and nozzle, and this pressure deficiency produces , 
decelerating force on the air column so that the mass flow tend 
to return to the equilibrium position. Such an intersection i 
called stable. On the other hand, for an intersection of the typ 
A,C (Fig. 4) a displacement to Point 1 produces an excess of avail 
able pressure over that required which applies an accelerating 
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Fic. 6. 


force to the air stream which tends to move the operating point 
to 2, and so on. Such an intersection is called unstable. If 
Curve C recrosses A, as in Fig. 5, vibration may occur. If it 
does not recross, the solution diverges steadily, causing blowout. 
Since most systems operate on the region of the compressor 
characteristic with negative slope, it must be shown tiat it is 
possible to have pressure-required curves with negative slope. 
In order to make a simpler problem computationally, pressure 
losses in the combustion chamber are neglected, and all dynamic 
pressures (1/29V*) are assumed to be small compared to am- 
bient pressure (Pg). Incompressible formulas can then be used. 
The pressure required by the nozzle is 
P, — Pq = (1/2)p2V2? 
and, since 
p2V2A» = Q 
and 
po = (11/72) pr 
this can be written 


P, — Pq = (1/2) (Q?/p:A2) (T2/Th) 


The slope dP, /dQ, is 


dP, " l (2 T, + 0? =) " 1 (20° ‘ a) 
dQ =. 2p AoT; secithalinde dQ} 2wA:T\QO § dQ 


The quantity in brackets will be negative if the combustion 
efliciency drops off with velocity enough that 
—dT»/dQ > 2(T2/Q) 
The data of reference 1 have been analyzed from this point of 


These data were taken from tests of an annular turbojet 
This is not too close to the one-dimensional model, 


view. 
combustor. 
since it is possible to have local instability at isolated points 
without general overall instability. As a matter of fact, three 
kinds of instability are reported as having occurred separately 
or in combination. One was observed visually as a flicker in the 
burning zone. One was observed as fluctuations in the exit tem- 
perature. The third was observed as noise and duct vibration 
Since the above analysis involves pressure stability, it should 
apply only to the third case, and so only such points have been 
studied. In Fig. 6 (data from reference 1, Figs. 8c and 9c) have 
been cross-plotted for constant fuel flow, and the points of un- 
stable combustion are indicated. Since the compressor char- 
acteristics are not given, it can only be said that the agreement 
with analysis is plausible. 

On the basis of this analysis, certain things can be pointed out 
which can help stabilize combustion. For one thing, operating 
on the steep (constant flow) portion of the characteristic curve 
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FORUM 


should help stability. For another, a metering system designed 
for constant heat release would maintain stable combustion if 
the metering system were fast enough. 


REFERENCE 


1 Childs, J. Howard, McCafferty, Richard J., and Surine, Oakley W., 
Effect of Combustor-Inlet Conditions on Performance of an Annular Turbojet 
Combustor, N.A.C.A. Report No. 881, 1947 


Errata—‘‘ Finite Deflections of Sandwich 
Plates’’* 


E. Reissner 

Department of Mathematics, Massachusetts Institute of Technology, 
Cambridge, Mass. 

November 3, 1949 

be THIS PAPER three of the equations should be corrected as 


follows: 


M, = (Ni+ — Nz-)(h/2), My = (Ny Ny) (h/2) | 
} (28) 
H = (S+ — S-)(h/2) } 
OV: OV, _ Ow 42 O*w +N O*u ‘ 
p Ox oy * ox? Oxoy a oy? 
07e 07e 07e 
M. — + 2H + M,— =0 (35) 
Ox? Oxoy oy? 
4 h y 07e 495 0*e +N o% 4 
q- iV a 5s yi . 
. 7 «4 ox? oxdy ” dy? 
1 ow o*w o*w 
M; + 2H + M, = 0 (36) 
( Ox? oxody : =) ” 


The error committed was in the nature of an oversight in 


algebra. I am very much indebted to Prof. C. T. Wang, of 
New York University, for discovering the necessity of the above 
corrections. As the terms that have been changed in Eqs. (35) 
and (36) were not used in the remainder of the paper, none of the 
subsequent results is affected by these changes. 
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Detonation and Static Electricity in the 
Tropical Atmosphere 


G. Broersma 

General Secretary, Organisation for Scientific Research, Konings- 
plein, Z. 11, pav., Batavia 

October 31, 1949 


ib 1941, THE AUTHOR submitted an explanation of the experi- 
ments of Corbitt and Clark! on aero-engines detonating when 
electrostatically cleaned air was admitted to the intake, respec- 
tively,? of experiences of aero-engines detonating during flight in 
or after thunderstorms in the tropics. This explanation was not 
published until 1946 and 1947.%* It probable that 
ozone had been formed in both cases, which induced and/or 
accelerated chain reactions so that detonation occurred 
Recently, when taking stock of literature that had appeared in 
Indonesia before the war and which had survived, the author 
found some figures on measurements of electrostatic charges 
In these measurements, de 


seemed 


occurring during tropical rains.‘ 
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Kruyff found positive charges up to 2,362 c.g.s. and negative 
charges up to 2,324 c.g.s. The largest specific charge found was 
4.7 c.g.s. per cu.cm. 

Small showers had higher specific charges than heavy showers. 


REFERENCES 

1 Corbitt, H. E., and Clark, N. J. Electrostatic Precipitation for Aircraft, 
Aero Digest, Aviation Engineering, p. 132, December, 1940. 

2 Bioersma, G., Detonation, Natuurwetenschappelijk Tijdschrift, Vol. 
102, p. 80, 1946 

8 Broersma, G., Detonatie, Polytechnisch Tijdschrift, p. 519a, 1947. 

‘de Kruyff, J. P. H., Metingen van electrische lading van regens te Soer 
abaia, Handelingen van het Achtste Nederlandsch-Indisch Natuurweten 
schappelijk Congres te Soerabaia, 140, July 1938, 20-23. 


Second-Order Two-Dimensional Supersonic 
Airfoil Theory 


Bruno W. Augenstein 
The RAND Corporation, Santa Monica, Calif. 
October 20, 1949 


N DETERMINING the higher order correction terms to the linear 

supersonic airfoil theory of Ackeret, the usual procedure has 
been‘ to expand either the two-dimensional oblique shock equa 
tions or the Prandtl-Meyer equation into a power series in 0, the 
flow deflection angle.!. A third method to be presented here 
consists in replacing the airfoil or wedge by sources producing 
either isentropic waves of finite amplitude or weak shock waves, 
thus extending Taylor’s result.2 The flow can then be thought 
of as being essentially one-dimensional if only outgoing waves 
are present. The computations will be carried out to terms in 
62 


First, by utilizing the isentropic relations 


p (*) EP vot ( ; (y= 
Po po . To i a 


where a = speed of sound, the one-dimensional equation of con- 
tinuity and momentum are: 


l(y — 1)/2] auz + a, + uaz = 0 


(1) 
u, + uuz + [2a/(y — 1)Jaz = of 


where ¥ is the ratio of specific heats and un is the particle velocity. 
Substituting @ = a(u) (which is possible since the motion is sup- 
posed isentropic), and multiplying the second equation by da/du 
and subtracting from the first, we get: 
da/du = +(y —1)/2 or a =a + {(y — 1)/2]u 

Considering only compression flows for convenience, the positive 
sign isappropriate. Hence: 

P/Po = il + [(y 1) /2] (u/ao)} 7/7 — YD (2) 
The particle velocity must be normal to the wave front and of 
such magnitude as to equal the component of free-stream velocity 
normal to the airfoil, @ sin 6. Therefore, x = a sin @/cos @ = 
u9/cos@,@small. (See Fig. 1.) 

The wave velocity itself is given by: 


T +1 +1 
iV => 1 a + ay + € )s = a + ¢ )s 
2L 2 4 


so that: 


( 42 + ] ) 1 ea +1f 0 ) ; 
sin 7 = dy ] “= (< 
. . 7 Vu 4 \cos¢ 


By expanding Eq. (2), we find, to terms in 6?: 
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P-—Po._ AP 26 v++1 
yPoM? gq Mecos¢ 2 
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cos? @ ‘ 


Putting 0, = sin~! (1/M) + 6, where 6 is small, into Eq. (; 
gives an equation for 5, which to the necessary order of approxi 


mation is solved b 


y: 


— y+ ‘)o M ) 
} M?—-1 


0 oy i £. y+ ‘ya = 
, = sin 
ileal M 4 M? — 


Also, by keeping only the lowest order terms in @, we find, fro, 


o = Oy — 0: 


— cos? ¢ 


This equation is solved by iterating, the first approximatig; 


= 2V 1 — cos? @ cos @ X 


4 (y — 3)M 
<0 
4(M?2 — 


) 


aa | 
1) Jf 


being cos @ = V (MM? — 1)/M; the next approximation is: 








jelta wi 
V M2 — 1S, _9V M1 [(y — 3)M? + 4)) plan for 
cos ¢ = Vv ) 2 2(M? — 1)? ( side of t 
verse t 
Substituting into Eq. (4) gives: assumpt 
AP 2 | A (vy + 1) — 4244 
= 0 oO 9 2 2 5 
q VM? —-1 2(M? — 1)* 
which is the Busemann relation 
If, on the other hand, we consider shock waves, we have the 
well-known relations: 
: : ¥+1 e 2 
AP = pWu, W = e+ [te sags 
} \ } 
To the second order we therefore get, as above: 
AP 20 + | (i ; 
= 4. , equati 
¢ AM cos $ cos? 
d ° ? metho 
and for the shock angle: No 
sin Oy = (1/M) + [(y + 1)/4] (@/cos ¢) the pr 
and hence Eq. (8) follows again. would 
The above method can be extended to unsteady airfoil prob quent] 
lems, but the work becomes virtually intractable unless we as howev 
sume that @y is constant during the motion. tion of 
may | 
REFERENCES ae 
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On Simple Plan Forms in Supersonic Flow 


John W. Miles 
Department of Engineering, University of Calitornia at Los Angeles 


September 26, 1949 


[° A RECENT ITEM! in this column, we pointed out that the prob- 


lem of nonsteady motion of a “wide’’ delta wing could be 


reduced to an equivalent two-dimensional problem. For the 
calculation of the force or moment coefficients, this reduction 
establishes the validity of strip theory (i.e., the calculation of the 
flow at any spanwise station as if it were two-dimensional). 

The purpose of this note is to point out that strip theory is 
yalid for the calculation of the forces and moments on any simple 
plan form with a straight trailing edge in either steady or un- 
steady flow. (A simple plan form is one for which the compo- 
nent of free-stream velocity normal to the boundary of the plan 
form in its own plane is everywhere supersonic, so that the zone 
of action of the plan form has no region in common with its zone 
of influence.) This follows directly from the result for the wide 
jelta wing, for we may circumscribe a wide delta about any simple 
plan form and regard the prescribed normal velocity as null out- 
side of the simple plan form. The restriction toa straight (trans- 
verse to the flow) trailing edge follows from the fact that the 


assumption of zero normal velocity aft of the trailing edge does 
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not meet the physical requirement of zero pressure jump across 
the wake. Our result may be stated mathematically as follows: 
If cp,(«,y,t) is the true three-dimensional pressure coefficient for 
a simple plan form with a straight trailing edge, if f(x) is an 
arbitrary (but finite) weighting function, and if S is the plan 


form, we have: 


(a + by)f(x) cp;(x,y,t)dx dy = 


(a + by)f(x) cp.(x,y,t)dx dy (1) 


Ss 


where ¢»,(x,y,/) is calculated on the assumption that the flow at 
any section is two-dimensional. In the case of steady flow, this 
result may also be proved with the aid of the reversed flow 
theorem for simple plan forms.” It is important to note that Eq- 
(1) is valid only for a (at most) linear, spanwise (y direction) 


weighting function. 
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Application of Simultaneous Equation Machines to 
Aircraft Structure and Flutter Problems 


(Continued from page 113) 


CoNcCLU 


equation machine varies as n*, while mathematical 
methods vary as n*. 

No check was made against I.B.M. solutions, since 
the problems were so small that undue preparation time 
would have been involved and the comparison, conse- 
quently, unfair to the digital machine. It is believed, 
however, that, because of the flexibility of the explora- 
tion of the problems and the intimate way the machine 
may be integrated into an engineering department, 
simultaneous equation analogy computers have a place 
alongside of the digital machines. Digital machines are 
useful in preparing the data for insertion in the equation 


machine. 


SIONS 


The problems presented here are only an example of 
what may be done with equation machines, and many 
problems, some undoubtedly not yet envisaged, may be 
solved with useful results. 

The accuracy is adequate for engineering practice, 
and the accuracy improvement methods described not 
only provide a means of improving it if required but also 
an independent check of the calculations. 


A machine with a capacity of between 10 and 20 real 
equations, equipped with provisions for complex X's, is 
recommended as being the proper size for an aircraft en- 


gineering department. 
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Enratum 


The Readers’ Forum item “Supersonic Airfoils Simplified’? (JouURNAL, Vol. 17, No. 1, p. 61, 
January, 1950) was incorrectly attributed to James E. Broadwell. 
H. Upson, Professor of Aeronautical Engineering, University of Minnesota, Minneapolis, Minn. 


The correct author is Ralph 
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